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Abstract 

We sketch some of the different roles played by Whitham times in connection with averaging, 
adiabatic invariants, soliton theory, Hamiltonian structures, topological field theory (TFT), 
Seiberg-Witten (SW) theory, isomonodromy problems, Hitchin systems, WDVV and Picard- 
Fuchs equations, renormalization, soft supersymmetry breaking, etc. 



Contents 



1 INTRODUCTION 



2 BACKGROUND 



2.1 Riemann surfaces and BA functions 



2.2 General remarks on averaging 

2.3 Averaging with ip*if) 

2.4 Dispcrsionlcss theory 10 



3 ISOMONODROMY PROBLEMS 



4 JMMS EQUATIONS 



5 GAUDIN MODEL AND KZ EQUATIONS 



6 ISOMONODROMY AND HITCHIN SYSTEMS 



17 

22 
27 
30 



7 WHITHAM AND SEIBERG-WITTEN| 40 

7T Background from \8$ 41 



7.2 Connections to [168fl 47 



8 SOFT SUSY BREAKING AND WHITHAM 



S.l Remarks on susy 



49 

49 

S.2 Soft susy breaking and spurion fields 53 



|9 RENORMALIZATIO N 
|9.1 Contact terms . . . 



56 

61 



1 



10 WHITHAM, WDVV, and PICARD-FUCHS 






10.1 ADE and LG approach 










10.2 Frobcnius algebras and manifolds 


10.3 Wittcn-Dijkgraaf-Vcrlindc-Vcrlindc (WDVV) equations 



1 INTRODUCTION 



Whitham theory arises in various contexts and the corresponding Whitham times play roles such as 
deformation parameters of moduli, slow modulation times, coupling constants, etc. When trying to 
relate these roles some unifying perspective regarding Whitham theory is needed and we try here 
to provide some steps in this spirit. The idea is to display Whitham theory in many (not perhaps 
all) of its various aspects (with some sort of compatible notation) and thus to exhibit (if not always 
explain) the connections between various manifestations of Whitham times. Some background work 
in_this dkectionap^ears in RMMM H 11 jlMMjl 115 JUL 11 11 PL , H 
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165], {Teg, [168[ [L6|, flTg ^175| |176|, |188|, |gg |190|, |191|, [Tgg |193|, |194|, |195|, |198[ [L9g foff Symbols 
such as (•), (Jit), (A), (XIII), etc. may be repeated in different sections but references to them are 
intrascctional. 



2 BACKGROUND 



Perhaps the proper beginning would be Whitham's book [199 or the paper |£4j where one deals 
with modulated wavetrains, adiabatic invariants, etc. followed by multiphase averaging, Hamiltonian 
systems, and weakly deformed soliton lattices as in 56, pft 126]. The next stages bring one into 



contact with some current topics in theoretical physics such as topological field theory (TFT), strings, 
and SW theory to which some references are already given above. In Section 2 we begin with the 
averaging method following |2^, |78|, |126| plus dispersionless theory as in @ H |3l], [mj |l8| (which 
i s re l ated to TFT) ; then we develop relations to isomonodromy and to Hitchin systems following 
|P1 pEj, ITM Il7i| Il76l EM |92l [193, EM Eif] in Sections 3-6, while connections to Seiberg-Witten 



(SW) theory, the renormalization group (RG) and soft supersymmetry (susy) breaking appear in 
Sections 7-9. In Section 10 we add a few remarks about t he st r ucture of contact t erms in twisted 
Af = 2 susy gauge theory on a 4-manifold following Q |l37| , |l45| , |l47| , |l62[ |l89 . 



The averaging method for soliton equations goes back to [|74j and early Russian work (cf. ]56] , 
p7{ for references) but the most powerful techniques appear first in Jl26| (with refinements and 
clarification in |1 M). For background on KdV, KP, and RS see e.g. |TC|23l |35l |1 |23l p7|. 



2.1 Riemann surfaces and BA functions 

We take an arbitrary Riemann surface S of genus g, pick a point Q and a local variable 1/k near 
Q such that k(Q) — oo, and, for illustration, take q(k) — kx + k 2 y + k 3 t. Let D = P\ + ■ ■ ■ + P g 
be a non-special divisor of degree g and write ip for the (unique up to a multiplier by virtue of 
the Riemann-Roch theorem) Baker- Akhiezer (BA) function characterized by the properties (A) ip 
is meromorphic on £ except for Q where ^i(P)exp(—q(k)) is analytic and (*) after normalization 
?/> ~ exp(q(k))[l + near Q- (^) On S/Q, V nas on ly a finite number of poles (at the 
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Pi). In fact ip can be taken in the form (P € E, Pq 7^ Q) 



ip(x,y,t,P)=exp[ (xdfl 1 + ydfl 2 + tdfl 3 )} ■ rv/um, ^ t 2 ' 1 ) 

where dft 1 = dfc + • • • , dft 2 = d(fc 2 ) + • • • , dtt 3 = d(k 3 ) + .-.,U j = J B . dil\ Vj = J B . dfl 2 , Wj = 
j B dVt 3 (j = l,---,g), zq = — A(D) — K, and is the Riemann theta function. The symbol 
~ will be used generally to mean "corresponds to" or "is associated with"; occasionally it also 
denotes asymptotic behavior; this should be clear from the context. Here the dilj are meromorphic 
differentials of second kind normalized via J A dQj — (A,, Bj are canonical homology cycles) and 

we note that xdtt 1 +yd£l 2 +tdQ 3 ~ dq(k) normalized; A ~ Abel-Jacobi map A(P) — {Jp o dwk) where 
the dujk are normalized holomorphic differentials, fc = l,---,g, J. dujk = Sjk, and K — (Kj) ~ 
Riemann constants (2K = —A(Ks) where Ks is the canonical class of E ~ equivalence class of 
meromorphic differentials). Thus Q(A(P) + zq) has exactly g zeros (or vanishes identically. The 



paths of integration are to be the same in computing Jp dfl 1 or A(P) and it is shown in [^5, 23 



Po 

that ip is well defined (i.e. path independent). Then the £j in (*) can be computed formally 
and one determines Lax operators L and A such that d y ip = Lip with dtip = Aip. Indeed, given 
the £j write u = —2d x ^\ with w — 3^i9^^i — 3<9 2 £i — 39 x ^2- Then formally, near Q, one has 
(-d y + dl + u)tp = 0(l/k)exp(q) and (-d t +d 3 + (3/2)ud x + w)4> = 0(l/k)exp(q) (i.e. this choice 
of u, w makes the coefficients of k n exp(q) vanish for n = 0,1,2,3). Now define L = d 2 + u and 
A = d 3 + (3/2)ud x + w so d y ip — Lip and dtip = Aip. This follows from the uniqueness of BA 
functions with the same essential singularity and pole divisors (Riemann- Roch). Then we have, via 
compatibility L t — A y = [A,L], a KP equation (3/A)u yy = d x [ut — (l/4)(6uuj, + u xxx )] and therefore 
such KP equations are parametrized by nonspecial divisors or equivalently by points in general 
position on the Jacobian variety J(E). The flow variables x,y,t are put in by hand in (A) via q(k) 
and then miraculously reappear in the theta function via xU + yV + tW; thus the Riemann surface 
itself contributes to establish these as linear flow variables on the Jacobian. The pole positions Pi 
do not vary with x, y, t and (f ) u — 2d 2 log<d(xU + yV + tW + z ) + c exhibits Q as a tau function. 

We recall also that a divisor D* of degree g is dual to D (relative to Q) if D+D* is the null divisor 
of a meromorphic differential dCl = dk + ((3/k 2 )dk + ■ ■ ■ with a double pole at Q (look at ( = 1/k to 
recognize the double pole). Thus D + D* - 2Q - Jf s so A(D*) - A{Q) +K = - [A{D) - A(Q) + K] . 
One can define then a function ip* (x,y,t, P) — exp(—kx — k 2 y — k 3 t)[l + £*/k) + ■ ■ ■} based on 
D* (dual BA function) and a differential dCl with zero divisor D + D* , such that (p = ipip*dti is 
meromorphic, having for poles only a double pole at Q (the zeros of dtl cancel the poles of ipip*). 
Thus ipip*dCl ~ tpip*(\ + (f3/k 2 + • • -)dk is meromorphic with a second order pole at oo, and no other 
poles. For L* = L and A* = —A+ 2w - (3/2)u x one has then (d y + L*)ip* = and {d t + A*)ip* = 0. 
Note that the prescription above seems to specify for ip* (U — xU + yV + tW, = —A{D*) — K) 

- (xd^+ydtf+tdtl 3 ) Q(A(P) - U+z%) 

* ~ e ° ' ecA(F) + *s) [2 - 2} 

In any event the message here is that for any Riemann surface E one can produce a BA function 
ip with assigned flow variables x,y,t,--- and this ip gives rise to a (nonlinear) KP e quation with 



solution u lin eariz ed on the Jacobian J(E). For averaging with KP (cf. J23| |78|, 126 1) we can use 



formulas (cf. (gjj) and fl2.2|)) 

iP = e px+Ev+Qt ■ (P{Ux + Vy + Wt, P) (2.3) 
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^* = e - P x-Ey-nt . ( f ) *^_ Ux -Vy- Wt, P) (2.4) 

with </>, (/)* periodic, to isolate the quantities of interest in averaging (here p = p{P) 1 E = E(P), il = 
51(P), etc.) We think here of a general Riemann surface S s with holomorphic differentials du>k 
and quasi- moment a and quasi-energies of the form dp = dil 1 , dE = dfl 2 , dfl = dfl 3 ,--- (p — 
Jp g dfl 1 etc.) where the cKP = dflj — d(X^ + 0(X^ 1 )) are meromorphic differentials of the second 



kind. Following [126, 127 one could normalize now via 3? J A . dVt k = 3? J B ,dfl k = so that e.g. 



Uk = (l/27ri) § A dp and Uk+ g = — (l/2m) § B ^ dp (k = 1, • • • , g) with similar stipulations for Vk ~ 
§ dil 2 , Wk ~ § dil 3 , etc. This leads to real 2g period vectors and evidently one could also normalize 
via § A dft k = (which we generally adopt in later sections) or 3 § A dft k = Q 1 § B dfl k = (further 
we set Bjk = § B ^ dojj). 

2.2 General remarks on averaging 

Averaging can be rather mysterious at first due partly to some hasty treatments and bad choices of 
notation - plus many inherent difficulties. Some of the clearest exposition seems to be in Eq, U% |7q, 



1 57 1 for which we go to hyperelliptic curves (cf. also jUj |2^, ^5[). For hyperelliptic Riemann surfaces 
one can pick any 2g + 2 points Xj € P 1 and there will be a unique hyperelliptic curve E s with a 
2-fold map / : E s — > P 1 having branch locus B — {Xj}. Since any 3 points A^, Xj, Xk can be sent to 
0, 1, oo by an automorphism of P 1 the general hyperelliptic surface of genus g can be described by 
(2g +2) — 3 = 2g — 1 points on P 1 . Since / is unique up to an automorphism of P 1 any hyperelliptic 
Tig corresponds to only finitely many such collections of 2g — 1 points so locally there are 2g — 1 
(moduli) parameters. Since the moduli space of algebraic curves has dimension 3g — 3 one sees that 
for g > 3 the generic Riemann surface is nonhyperelliptic whereas for g = 2 all Riemann surfaces 
are hyperelliptic (with 3 moduli). For g = 1 we have tori or elliptic curves with one modulus r and 
g = corresponds to P 1 . In many papers on soliton mathematics and integrable systems one takes 
real distinct branch points Xj, 1 < j < 2g + 1, and oo, with Ai < A2 < • • • < X2 g +i < 00 and 

23+1 

^ 2 = H(X-X j ) = P 2g+1 (X,X j ) (2.5) 
1 

as the defining equation for S a . Evidently one could choose Ai = 0, A2 = 1 in addition so for g = 1 
we could use < 1 < u < 00 for a familiar parametrization with elliptic integrals, etc. One can take 
dX/ ix, XdX/fi, ■ ■ ■ , X 9 ^ 1 dX/^ as a basis of holomorphic differentials on S g but usually one takes linear 
combinations of these denoted by du)j, 1 < j ' < g, normalized via § A dujj = 8ij, with period matrix 
defined via § g dujj = Ily. The matrix If = (11^) is symmetric with 311 > and it determines 
the curve. Frequently in situations arising from KdV (Korteweg-deVries equation) for example one 
regards the intervals [Ai, A2], • • • , [A2 g +i, 00) as spectral bands and intervals (A2, A3), • • ■ , (X 2g , A2 9 +i) 
as gaps with the a, cycles around (X 2 i, A2i+i) (i — 1, ■ ■ ■ , g). One will also want to consider another 
representation of hyperelliptic curves of genus g via 

2g+l 

f x 2 =l[(X-X j ) = P 2g+2 (X,Xj) (2.6) 


where 00 is now not a branch point and in fact there are two points /x± corresponding to A = 00. 
We recall next some of the results and techniques of |74|, |157| w here one can see explicitly the 



nature of things. The presentation here follows |E3[. First from [157], in a slightly different notation 
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write q t = 6qq x - q xxx with Lax pair L = —d% +q, B = -Adl + 3(qd x + d x q), L t = [B, L], Lip = Xtp, 
and ipt = Bip. Let ip and <f> be two solutions of the Lax pair equations and set '5 = ^0; these are 
the very important "square eigenfunctions" which arise in many ways with interesting and varied 
meanings (cf. 0, |25|, EH], ||^, ||j|). Evidently \I> satisfies 

[-dl + 2(qd x + d x q)]$ = 4A5 X *; d t V = -2q x V + 2(q + 2A)9 X * (2.7) 



From (2.7) one finds immediately the conservation law (C): dt[^] + d x [6(q — 2\)^> — 29^] = 0. If 



one looks for solutions of (2.7) of the form ^(x,t, A) = 1 + J2T\^ j( x ^)]^ 3 as ^ ~~ * 00 then one 



obtains a recursion relation for polynomial densities 

a^+x = [-1^ + ( g a x + cU)]^- u = 1, 2, • • •); *o = 1 (2.8) 

Now consider the operator L = —d\ + q in L 2 (— 00,00) with spectrum consisting of closed intervals 
separated by exactly N gaps in the spectrum. The 2N + 1 endpoints Afc of these spectral bands 
are denoted by —00 < Ao < Ai < • • • < X2N < 00 and called the simple spectrum of L. They 
can be viewed as constants of motion for KdV when L has this form. We are dealing here with 
the hyperelliptic Riemann surface determined via -R 2 (A) = Ho^t^ — ^fc) (°° i s a branch point) and 
one can think of a manifold Ai of A^-phase waves with fixed simple spectrum as an ./V-torus based 
on 0j £ [0,27r) . Hamiltonians in the KdV hierarchy generate flows on this torus and one writes 
Q = Qn(0i, • • • , #iv) (the 9 variables arise if we use theta functions for the integration - cf. also 
below). Now there is no y variable so let us write 0j = xKj + twj (we will continue to use dujj 
for normalized hol omor phic differentials). For details concerning the Riemann surface we refer to 
H ||, H, II H, |m| and will summarize here as follows. For any q^ as indicated one can find 



functions fJ.j(x, t) via ^(x, t, A) = J^f^A — fij(x, t)) where fJ,j(x, t) £ [Aaj— i, \2j] and satisfies 

d xH = -2i{R{^)/ - Mi)); (2-9) 

2N 

dtnj = -2i[2(J2 Afc - 2^Mi)] ' (i*0*i)/II^ ~ ^ 

In fact the /ij live on the Riemann surface of i?(A) in the spectral gaps and as x increases \ij travels 
from \2j-1 to \2j on one sheet and then returns to \2j-1 on the other sheet; this path will be called 
the j //-cycle (~ Aj). In the present context we will write the theta function used for integration 
purposes as 6(z, t) = J2 m ez N exp[iri(2(m, z) + (m, Tin))] where z S and r denotes the N x N 
period matrix (r is symmetric with 3r > 0). We take canonical cuts Ai, Bi (i = 1, • • • , N) and let 
du>j be holomorphic diffentials normalized via J A diok = Sjk (the cycle Aj corresponds to a loop 
around the cut Aj). Then q^ can be represented in the form 

2N 

q N (x,t)=A + T-2d*loge(z(x,t);T); A = ^Aj-; (2.10) 



N 

( T i 3 ) = (j> du 3 ); r*j = -Tij\ r = -2^ <t Xdujj 



1 Bi 

andz(a;,t) = —2i[c JV (x—a:o)+2(Ac Ar +2c Ar_1 )t]+d where (c N )i = CiM arises from the representation 
duji — (^2 1 (kj\3~ )[dX/ R(X)] (d is a constant whose value is not important here). Then the wave 
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number and frequency vectors can be defined via H = —Ainr 1 c N and w — —8iirr ^Ac^ + 2c N x ] 
with 9j(x, t) = KjX + Wjt + 6 j (where the 6® represent initial phases). 

To model the modulated wave now one writes now q = qN(&i, • • • , 6n] A) where Xj ~ Xj(X, T) and 
A ~ (Aj). Then consider the first 2N+1 polynomial conservation laws arising from (2.7) - ( |2.8| ) and C 



for example (cf. below for KP) and write these as d t Tj(q) + d x Xj(q) = (explicit formulas are given 
in [jFI] roughly as follows). We note that the adjoint linear KdV equation (governing the evolution 
of conserved densities) is d t jj + d x jj — 6qd x jj = (jj ~ ViJj) and ( [j~§l ) has the form djj + i = 
(— (l/2)<9 3 + qd + dq)-fj. One then rewrites this to show that 6qd x jj = d x [6jj + i — 6q-fj + 3d 2 jj] so 
that the adjoint equation becomes 

d tlj + d[-23 2 l3 + 6q 7j - 67 i+x ] = (2.11) 



which leads to ( [2.12] ) and ( 2.16 ) below (after simplification of ( 2.11 )). Then for the averaging step 



write dt = £(?t, etc., and average over the fast variable x to obtain 

8 T < T 3 (q N ) > +d x < Xjitot) >= (2.12) 
(( [2.12 ) makes the first order term in e vanish). The procedure involves averages 



<Tj(q N ) >= UmL^oc-^jT J Tj(q N )dx (2-13) 

for example (with a similar expression for < Aj(qjv) >) and an argument based on ergodicity is 
used. Thus if the wave numbers Kj are incommensurate the trajectory {q^(x,t); x € (— oo,oo)} 
will densely cover the torus M. Hence we can replace x averages with 

< TMn) >= -^v j o ' ■ ■ J o T AMe)) J] Mj (2.14) 

For computational purposes one can change the 9 integrals to fj, integrals and obtain simpler calcu- 
lations. By this procedure one obtains a system of 2N + 1 first order partial differential equations 
for the 2N + 1 points Xj(X,T), or equivalently for the physical characteristics (k(X, T), w(X, T)) 
(plus < qw >)■ One can think of freezing the slow variables in the averaging and it is assumed that 
( [2.12 ) is the correct first order description of the modulated wave. 



The above argument may or may not have sounded convincing but it was in any case rather 
loose. Let us be more precise following |7^| . One looks at the KdV Hamiltonians beginning with 
H = H{q) = liniL^oo^^L) J_ L (q 2 + (l/2)ql)dx (this form is appropriate for quasi-periodic sit- 
uations). Then {/, g} = Um^^^iX^V) J^ L (Sf/5q)d x (Sg/5q)dx (averaged Gardner bracket) and 
q t = {q,H}. The other Hamiltonians are found via 

d 5 J^±l = {qd + dq - Itf)**^ (m > 0); ^° = 1 (2.15) 
oq 2 oq oq 

where jj ~ Vi/j ~ (dHj/Sq) (cf. here |7^|). It is a general situation in the study of symmetries 
and conserved gradients (cf. J32|) that symmetries will satisfy the linearized KdV equation (dt — 
6d x q + d x )Q = and conserved gradients will satisfy the adjoint linearized KdV equation (dt — 
Qqd x + d x )Q* = 0; the important thing to notice here is that one is linearizing about a solution q of 



G 



KdV. Thus in our averaging processes the function q, presumed known, is inserted in the integrals. 
This leads then to 



(2.16) 



with ( [2.12| ) holding, where < d 2 <p >= implies 



1 

< Xj >= LimL-taz — 



SH j+1 5H 
(-6—^ h QqN T — )dx 



Sq. 



N 



Sq 



N 



(2.17) 



In [frij the integrals are then simplified in terms of fi integrals and expressed in terms of abelian 
differentials. This is a beautiful and important procedure linking the averaging process to the 



Riemann surface and is summarized in [ 157 1 as follows. One defines differentials 

2 l ^ 3 l R{X) 



(2.18) 



Ci 2 = \-ix 



N+l 



dX 



where the Cj, Ej are determined via <f b 
that 

< * >~< T >r 



Oi = 



OO 

E 





<X>r 



1, 2, • ■ • , N). Then it can be shown 



OO 

E 





< Xj > 
(2/x)j 



(2.19) 



oo (/I 



with Six ~< T > (d£/£ 2 ) and < X > (d£/£ 2 ) - 12[(d£/£ 4 ) - fl 2 ] where [i = <T 2 
(VV?) 1/2 ) so d/i = -2^~ a de (df/£ 2 ) - -(£/2)dp ~ -{dn/ZJfl. Since A x = 0(fi N /fi N +^)dn = 
0(fi- (1/2) d(j,, Cl 2 = 0(n 1/2 )dfi (with lead term -(1/2)) we obtain < <5 >~< T >= 0(1) and 
< X >= 0(1). Thus (2.7), ( |2.8| ), (C) generate all conservation laws simultaneously with < Tj > 
(resp. < Xj >) giving rise to fl± (resp. fi 2 ). It is then proved that all of the modulational equations 
are determined via the equation 

d T Cli = I2d x &2 (2.20) 

where the Riemann surface is thought of as depending on X, T through the points Xj(X,T). In 
particular if the first 2N + 1 averaged conservation laws are satisfied then so are all higher averaged 
conservation laws. These equations can also be written directly in terms of the Xj as Riemann 
invariants via drpXj = SjdxXj for j = 0, 1, • • • , 2N where Sj is a computable characteristic speed (cf. 



(2.29)-( p~30 )). Thus we have displayed the prototypical model for the Whitham or modulational 
equations. 



Another way of looking at some of this goes as follows. We will consider surfaces defined via 
i?(A) = rii ff+1 (A — Aj). For convenience take the branch points A, real with Ai < •• • < A 2g+ i < 00. 
This corresponds to spectral bands [Ai, A 2 ], • ■ ■ , [A 2g +i, 00) and gaps (A 2 , A 3 ), • • • , (A 2g , A 2g +i) with 
the Ai cycles around the gaps (i.e. a« ~ (A 2 ^, A 2 ^ + i), i = 1, • • • , g). The notation is equivalent to 
what preceeds with a shift of index. For this kind of situation one usually defines the period matrix 
via iBjk = § B ojj and sets 



A k 



diOj = 2nSjk (j, k 



,3); duj 



<A Cj q K q - l dK 



1 



(2.21) 
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(cf. [[l5[ ^||). The Bi cycles can be drawn from a common vertex (Pq say) passing through the 
gaps (A2j, A2i+i). One chooses now e.g. p — J dp and Q = J dQ in the form 

P (A) = /,KA> = /»P = A- + X>A*- j ; <"» 

fi(A) - y dfi(A) = y ^Z- [ > dA; O =^bjA^; b = -3^ A * 

with normalizations J^ a,+1 dp(A) = f A ^ ,+1 dfl(A) = 0; i = 1, ■ ■ ■ , g. We note that one is thinking 

here of (*) V = exp[ipx + iQt] ■ theta functions (cf. ( |2.ip with e.g. p(A) = —i(log%p) x ; f2(A) = 
—i{log^) t . Recall that the notation < , > K simply means re-averaging (or ergodic averaging) and 
(log^) x ^< (logij))x >x z / z here since e.g. (logvp) x is not bounded. Observe that (*) applies 
to any finite zone quasi-periodic situation. The KdV equation here arises from I/ip = Atp, L = 
—d 2 + q, dtip = Aip, A — 4<9 3 — 6qd — 3q x and there is no need to put this in a more canonical form 
since this material is only illustrative. 

In this context one has also the Kruskal integrals Iq , • • • , li g which arise via a generating function 

X j 

p(A) = -i < (logi>) x > X =VA + Y / (2 ^ )2s+1 (2-23) 
where I a =< P s > x = ~P S (s = 0, 1, ■ • ■) with -iQogip)* = \/A + EH^A 2 ^) 2 ^ 1 ]. Similarly 

A J OO Q 

- i{logi>) t = -i— = 4(VA) 3 + V (2.24) 

V J ' V> o (2VA) 2s+1 

and one knows dtP s = d x fl s since (♦) [{logip) x ] t — [(logip)t] x . The expansions are standard (cf. 
p3| , ^3, ^5|, [n| ) . Now consider a "weakly deformed" soliton lattice of the form 

6{t\B)= exp(-^J2 B 3kn 3 n k +tY,n ] T J ) (2.25) 

— oo<ni<-<n m <oo j,k j 

with the Aj (i = 1, • • • , 2g+ 1) (or equivalently the parameters u l = Li; i = 0, • • • , 2g) slowly varying 
functions o f cc, f (e.g. u l = u l (X, T), X = ex, T = et, i = 0, 1, • • • , 2m). Now one wants to obtain a 
version of ( [2.20 ) directly via (4). Thus insert the slow variables in (4k) and average, using ec?x or 
ec?T in the external derivatives, to obtain 



d T {logi>) x = d x {logip) t (2.26) 
or 8tp(A) = dxE(A). Then from ( |2.22| ) differentiating in A one gets 

d T dp = d x dQ (2.27) 



Now (recall dtP s — d x Q s ) expanding (2.27) in powers of (vA) 1 one obtains the slow modulation 
equations in the form 

d T u s =d x n s ( a = 0,->-,2g) (2.28) 
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where £l s is a function of the u l (0 < i < 2g). This leads to equations 

d T A t =v l (A 1 ,---,A 2g+1 )d x A. l (i = l,---,2g + l) (2.29) 
for the branch points as Riemann invariants. The characteristic "velocities" have the form 



dn 

dp 



A=Ai 



(2.30) 



To see this simply multiply (2.27) by (A — A^) 3 / 2 and pass to limits as A — > A; 



2.3 Averaging with t/;*ifi 



Let us look now at [126] but in the spirit of [E3l [781]. We will only sketch this here and refer to 



]23| for more detail. Thus consider KP in a standard form 3a 2 u yy + d x (iut — 6uu x + u xxx ) = 
via compatibility [d y — L, dt — A] =0 where L = a^ 1 (d 2 — u) and A — d 3 — (3/2)ud + w (a 2 = 1 
is used in [Q which we follow for convenience but the procedure should work in general with 
minor modifications - note d means d x and u — » — u in the development of (^J)). We have then 
(d y — L)ijj = with (dt — A)ip ~ and for the adjoint or dual wave function ip* one writes in [127] 
i p*L = —dy^p* with ip*A — d t ip* where ip*(fd j ) = (—d) j (ijj*f). We modify the formulas used in 
jl26| (and @) in taking (cf. J23| - 4> and 4>* are periodic) 

V> = e px+Ev+nt ■ <f>(Ux + Vy + Wt, P); if}* = e -P*-E V -at . ^*^_ Ux -Vy- Wt, P) (2.31) 



where p = p(P), E = E(P), — £l(P), etc. (cf. (2.3) - (2.4)), which isolate the quantities needed 
in averagi ng. The arguments to follow are essentially the same for this choice of notation, or that 
in |7q] or |l26f|. Now one sees immediately that 



(tp*L)ip = Tp*Lip + a^iV) + d 2 x (^*L 2 ^) + ■ ■ 



(2.32) 



where e.g. L r = ((-l) r /r!)(d r L/d(d) r ). In particular L 1 = -2d and L 2 = 1 while A 1 = -3d 2 + 
(3/2)u, A 2 — 3d, and A 3 = —1. We think of a general Riemann surface E g . Here one picks 
holomorphic differentials dto^ as before and quasi-momenta, quasi-energies, etc. via dp ~ effii, dE ~ 
dQ.2, dfl ~ (KI3, • • • where A ~ k, p = Jp o (Kli, etc. Normalize the dflk so that $t f A . dfik = = 
3ft f B dfi/c (cf. remarks after (2.4)); then U, V, W, ■ ■ ■ are real 2g period vectors, and one has BA 
functions tj)(x,y,t, P) as in (2.31). As before we look for approximations based on uq(xU + yV + 
tW|I) = Uo(dj, Ik)- For averaging, 6j ~ xUj +yVj +tWj, 1 < j < 2g, with period 27r in the 9j seems 



y ji "g+j 



natural (but note 
becomes < <p >= (1/(2tt) 2 s /■••/ 
bounded. In [Q one thinks of (f>(xU 



Uj, etc.). Then again by ergodicity < <p >x— Hm_L_ +00 (l/2L) J_ 1 



and one notes that < d x (f> >= automatically for q 
) with c/) x = J2 Ui(d4>/d0i) and J • • • / (d^/dO^d 2 ^ = 0. 



Now for averaging we think of u ~ ua(jS\I) with (S, I) ~ (S, I)(X, Y,T),d x S = U, d Y S = V, 
and drS = W. We think of expanding about uq with d x — > d x + edx- This step will cover 
both x and X dependence for subsequent averaging. Then look at the compatibility condition 
(ft) : d t L — d y A + [L, A] = 0. As before we will want the term of first order in e upon writing e.g. 
L = Lq + eL\ + • • • and A = Aq + eA\ + • • • where slow variables appear only in the Lo, A$ terms. 
Details are indicated in p3 and after some calculation (Jfr) becomes 



d t U - d v A + [L , A ] + eidtL, - d v A x + d T L - d Y A + [Lq, A x ) + [L u A ]} + 0(e 2 



(2.33) 
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Now the dtLo — d y Ao + [Lq,Aq] term vanishes and further calculation shows that to make the 
coefficient of e vanish one wants 

dtLi - d v A l + [L , Ai] + [Lx, A ] + F = 0; F = d T L-d Y A- [L l dxA - A 1 d x L) (2.34) 

However via ergodicity in x, y, or t flows, averaging of derivatives in x, y, or t gives zero, so from 



(2.34) and further computation we obtain the Whitham equations in the form < ip*Fip >= 0. In 
order to spell this out in Jr8j one imagines X, Y, T as a parameter £ and considers £(£), A(£), etc. 
and ip ~ (with / ~ (/&) being any additional factors dependent on the slow times) so that 

</-(£) = e P(i)*+E(i)y+m)t . cj>(U(Z)x + V(0y + W(0t\I(0) (2-35) 

but one keeps ip* = exp(—px—Ey—Clt)(fr*(—Ux—Vy—Wt\I)&xedm£ (i.e. assume p, E, fi, U, V, W, I 
fixed in ip*). We recall that one expects \k = Xk(X,Y,T) etc. so the Riemann surface varies with 
£. The procedure here is somewhat contrived but seems appropriate for heuristic purposes at least. 



Also recall that x, y, t and X, Y, T can be considered as independent variables. Note also from ( 2.35 ) 
for P fixed (0 ~ xU + yV + yW) 

W>(£)l«=o = (px + Ey + f2t)V>*V> + (2.36) 

+{Ux + Vy + Wt) ■ iP*d ip + I ■ ip*diil] 

where / ~ df /d£. In |Q one assumes that it is also permitted to vary £ and hold e.g. the Ik constant 
while allowing say the P to vary. Now one computes terms like d^[dt(ip*ip{O)]\^=0j averages, and 
subsequently inserts X, Y, T for £. Taking into account compatibility relations 

d Y U = d x V; d T U = d x W; d T V = d Y W (2.37) 

after considerable calculation one arrives at the Whitham equations in the form 

PT = n x ; py=E x ; E T = n Y (2.38) 



(cf. Section 7 for general forms). We feel that this derivation from J23[, based on |126| , is 
important since it again exhibits again the role of square eigenfunctions (now in the form ip*tp) in 



dealing with averaging processes (a derivation based on (2.26) also seems important). In view of the 
geometrical nature of such square eigenfunctions (cf. ]32] | for example) one might look for underlying 
geometrical objects related to the results of averaging. Another (new) direction involves the Cauchy 
kernels expressed via ij)*ip and their dispersionless limits (cf. J27], |2q|). It is also proved in p3j that 
dp =< ipip* > dO, dE ~ - < ^L 1 -^ > d(l, and dil — - < ^A 1 ^ > dti. This shows in particular 
how the quantity ipip* determines the Whitham differentials dp, dE, dQ, etc. 

REMARK 2.1. Since D + D* — 2oo <~ Ks, ipip* is determined by a section of Ky, (global point 
of view) but it is relations based on < ^L 1 ^ >, < ?/;*A 1 -0 >, etc. (based on the Krichever averaging 
process) which reveal the "guts" of ipip* needed for averaging and the expression of differentials. 



2.4 Dispersionless theory 

We give next a brief sketch of some ideas regarding dispersionless KP (dKP) following mainly 
188| to which we refer for philosophy. We will make various notational adjustments as 
we go along. One can think of fast and slow variables with ex = X and et n = T n so that d n — > ed/dT n 
and u(x,t n ) — > u(X,T n ) to obtain from the KP equation (l/4)u xxx + 3uu x + (3/4)9 _1 9|u = the 
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equation dru = 3udxu+ (3/4)9 _1 (9 2 m/9T|) when e — > — > (l/e)9 _1 ). In terms of hierarchies 
the theory can be built around the pair (L,M) in the spirit of |2^, |32| 188 1 . Thus writing (t n ) for 
(x,t n ) (i.e. a; ~ ti here) consider 

oc oo oo 

L e = ed + Y,u n+ i(e,T)(ed)- n ; M e = ^ nT n I£- 1 + £ u„+i(e, T)^"- 1 (2.39) 
l 11 

Here £, is the Lax operator L = d + Y^T w n+ i<9~™ and M is the Orlov-Schulman operator defined 
via ip\ = Mip. Now one assumes u n (e,T) — U n {T) + 0(e), etc. and set (recall Lip — Xip) 



1 + 



explf^— X l = exp(-S(T,X) + 0(l 



T = exp{ \F(T) + O 



(2.40) 



We recall that d n L = [B n ,L], B n = L\, d n M = [B n ,M], [L,M] = 1, Lip = Xip, d x ip = Mip, 
and ip = t(T- (l/nA"))ea;pE~r n A' i ]/T(r). Putting in the e and using <9„ for <9/<9T„ now, with 
P = Sx , one obtains 

oo oo 

A = P + ^C7„ +1 P-"; P = A-^PA" 1 ; (2.41) 

1 1 

OO oo 

M = Y, nT n \ n - 1 + V n+1 \- n - x ; d n S = B n (P) => d n P = 8B n (P) 
i i 

where d ~ dx + (dP/dX)dp and M — > M. Note that one assumes also v i+ i (e, T) = V*+i (T) + 0(e); 
further for B n = J2o bnmd" 1 one has B n = Y^q b nm P m (note also B n = L n + CT i ^~' J )- We ust a 
few additional formulas which are easily obtained (cf. |2q| ); thus, writing {A, P} = OpAdA — dAdpB 
one has 

9„A = A}; d n M = {B n , M}; {A, M} = 1 (2.42) 
Now we can write S = ^n^ n + Si° Sj'+i^ ~ J with d m Sj+i = aj 1 , V n +\ — —nS n +i, and c^S 1 = 

ov n+1 dd n F 



M (of -*af). Further 



B n = A" + dnS j+ i\- j ; dS n+1 ~ -P„ 



(2.43) 



We sketch next a few formulas from [121]. First it will be important to rescale the T n variables 
and write t' — nt n , T' n — nT n , d n — nd' n = n(d/dT n ). Then 



A" K 
d' n S=^; d' n X = {Q n ,X} (Q» = ^)5 

d' n P = dQ n = dQ n + dpQ n dP; d' n Q m - d' m Q n = {Q n , Q m } 



(2.44) 



Think of (P, X, 7£), n > 2, as basic Hamiltonian variables with P = P(X, T n ). Then -Q„(P, X, T n ) 
will serve as a Hamiltonian via 



(2.45) 
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(recall the classical theory for variables (q,p) involves q = dH/dp and p = —dH/dq). The function 
S(X,X, T n ) plays the role of part of a generating function S for the Hamilton-Jacobi theory with 
action angle variables (A, — £) where 

A" 

PdX + Q„dT' = -£,d\ - K n dT' + dS; K n = -R n = ; (2.46) 

n 

-^7 = Ki = 9(Rn = 0; -jjp =i' n = -d\Rn = -A" -1 

(note that X' n — ~ d' n A = {Q n , A}). To see how all this fits together we write 
jp ftp ap 

m = 3 » p + ox m ~ dQ " + ox** = aa - + apa - a " + dp *'« {2A1) 



This is compatible with (2.45) and Hamiltonians — Q n . Furthermore one wants 

Sx = £ S x = P; d'J = Q n - Rn (2.48) 



and from ( 2.46 ) one has 

PdX + Q n dT' n = -£dX + R n dT' n + S x dX + S x d\ + d' n SdT n (2.49) 

which checks. We note that d' n S = Q n = B n /n and Sx = P by constructions and definitions. 
Consider S = S- A " T r7«- Thcn Sx = S x = P and S' n =S' n - R n = Qn ~ Rn as desired with 
i = S x = S x - ET K^ n ~ l - lt follows that £ ~ M - Y% K^' 1 = X + tS+iA - * -1 . If W is 
the gauge operator such that L = \VdW~ 1 one sees easily that 

(CO \ / CO \ 

fca*0 fc_1 J = ( G + ^ tecfcA* -1 J ^ (2.50) 

from which follows that G = WxW^ 1 — > £. This shows that G is a very fundamental object and 
this is encountered in various places in the general theory (cf. |25|, 51, |202fl). 



REMARK 2.2. We refer here also to |24|, Q for a complete characterization of dKP and the 
solution of the dispersionless Hirota equations, and will sketch this here. Thus we follow ]24| (cf. 
also [121, 188|) and begin with two pseudodiffcrcntial operators (d = d/dx), 

OO OC 

L = d + J2 u n+id~ n ; W =l + ^2w n d- n (2.51) 
i i 

called the Lax operator and gauge operator respectively, where the generalized Leibnitz rule with 
d~ 1 d = dd^ 1 = 1 applies 

d i f = jr( t \&f)d i -i (2.52) 

3=0 

for any i 6 Z, and L = \VdW~ 1 . The KP hierarchy then is determined by the Lax equations 

{d n = d/dtn), 

d n L = [B n , L] = B n L - LB n (2.53) 

where B n = L™ is the differential part of L n = L\ + L n _ = £?d l + ElL^- 0ne can also 
express this via the Sato equation, 

d n WW~ x = -UL (2.54) 
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which is particularly well adapted to the dKP theory. Now define the wave function via 

oo oo 

il) = We !L = w{t,\)£\ £ = ^i„A"; w(t, A) = 1 + ^ w n (t)\~ n (2.55) 

1 i 

where t\ = x. There is also an adjoint wave function ip* — W*~ x exp(— £) = w*(t, A) exp(— £), 
w*(t, A) = 1 + w i(t)^~ l i an d one nas equations 

L^ = \^- d n ^ = B n ^; L*f=Af; = -B* n ^* (2.56) 



Note that the KP hierarchy ( 2.53 ) is then given by the compatibility conditions among these equa- 
tions, treating A as a constant. Next one has the fundamental tau function r(t) and vertex operators 
X, X* satisfying 

X(X)r(t) ^G.(X) T (t) ^-[A" 1 ]) 
A) = — ^_ = = W) ; (2.57) 

X*(A)r(t) e-tG + {X)T{t) e^rjt + [A' 1 ]) 

V {t, A = 777 = 777 = 777 

r(t) r(t) r(t) 

where G±(A) = exp(±£(<9, A -1 )) with 8 = (%, (1/2)%, (1/3)%, ■ • •) and i± [A" 1 ] = (ti ± A-\t 2 ± 
(1/2)A~ 2 , • ■ •). One writes also 

(oo \ oo 

X>" A j = J Exj(ti,t 2 ,---,t j )X j (2.58) 

where the %j are the elementary Schur polynomials, which arise in many important formulas (cf. 
below) . 

We mention now the famous bilinear identity which generates the entire KP hierarchy. This has 
the form 

i/j(t,X)ili*(t',X)dX = (2.59) 



where § (-)dX is the residue integral about oo, which we also denote Res\[(-)dX]. Using (2.57) this 
can also be written in terms of tau functions as 



T 



(i - [A- 1 ])^' + [A- 1 ])^*'^-^''^ = (2.60) 



This leads to the characterization of the tau function in bilinear form expressed via (t — > t — y, t' —* 

t + y) 

r£ X n(-2y)xn + i(d)e^T « 9 * j r • r = (2.61) 

where d™a ■ b = {d m /ds 7 p)a{t J + s s )b(U -gj)| s =q and d = (%, (1/2)%, (1/3)%, • • •)■ In particular, 
we have from the coefficients of y n in (2.61), 

%%t-t = 2 X „ +1 (%V-t (2.62) 

which are called the Hirota bilinear equations. One has also the Fay identity via (cf. ||, |3l], pjj 
- c.p. means cyclic permutations) 

^(so - «i)(«2 - s 3 )r(t + [s ] + [ Sl ])r(t + [s 2 ] + [s 3 ]) = (2.63) 

C . p ■ 
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which can be derived from the bilinear identity ( 2.60| ). Differentiating this in sq, then setting 
so = S3 = 0, then dividing by s\s 2 , and finally shifting t — > t — [S2], leads to the differential Fay 
identity, 



T(t)dr(t + [si] - [s 2 ]) - r(t + N - [s 2 ])dr(t) 

= (sr 1 - S2 1 ) Ht + M - [s 2 })r(t) - r{t + [«i])r(t - [s 2 ] 



(2.64) 



The Hiro ta eq uations (2.62) can be also derived from (2.64) by taking the limit si — > s 2 . The 



identity (2.64) will play an important role later. 



Now for the dispersionless theory (dKP) one can think of fast and slow variables, etc., or averaging 
procedures, but simply one takes t n — > et n = T n (pi = x — > ex = X) in the KP equation u t = 
(l/A)u xxx + 3uu x + (3/4)d~ 1 u yy , (y =t 2 , t = t 3 ), with d n —> ed/dT n and u{t n ) — > U(T n ) to obtain 
d T U = 3UUx + (3/4)9~ 1 C/yy when e — > (d — d/dX now). Thus the dispersion term u xxx is 
removed. In terms of hierarchies we write 



L c = ed + J2u n+1 (T/e)(ed)- n 



(2.65) 



and think of u n (T/e) = U n (T) + O(e), etc. One takes then a WKB form for the wave function with 
the action S 

V> = exp ~S{T,\) (2.66) 

Replacing now d n by ed n , where d n — d/dT n now, we define P = dS = Sx ■ Then e l d l ip — > P l i> as 
e — > and the equation Lip = \ip becomes 



\ = P + J2u n +iP- n ; P = \-J2 p >+i y 



(2.67) 



where the second equation is simply the inversion of the first. We also note from d n ijj = B n ip = 
Yso b n 7n(£d) m ip that one obtains d n S — B n (P) = A™ where the subscript (+) refers now to powers 
of P (note ed n i>/<ip -> d n S). Thus B n = L% -> B n (P) = A™ = Eo 5 «™ pm and thc KP hierarchy 
goes to 

d n P = dB n (2.68) 



from (2.57) in the limit e — » as 



which is the dKP hierarchy (note c^S" = i3 n =>■ 9 ra P = dB n ). The action 5 in (2.66) can be computed 

(2.69) 



s=yr„A™-y^A- 

Z ^ Z — ^ 



1 1 

where the function F = F(T) (free energy) is defined by 



r = exp 



>(T) 



(2.70) 



The formula ( 2.69 ) then solves the dKP hierarchy ( |2.41 ), i.e. P = B\ = dS and 



B n = d n S = A' 



Er nm 
711 



(2.71) 
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where F nm = d n d rn F which play an important role in the theory of dKP. 



Now following |188| one writes the differ ential Fay identity ( 2.64 ) with ed n replacing d n , looks 
at logarithms, and passes e — > (using (2.7C )). Then only the second order derivatives survive, and 
one gets the dispersionless differential Fay identity 



E /' 

m,n— 1 



-m \—n r mn 




(2.72) 



Although (2.72) only uses a subset of the Pliicker relations defining the KP hierarchy it was shown 
in 188] that this subset is sufficient to determine KP; hence (2.72) characterizes the function F 
for dKP. Following [^4| [h], we now derive a dispersionless limit of the Hirota bilinear equations 
( [2.62 ), which we call the dispersionless Hirota equations. We first note from ( 2.69| ) and (2.67) that 
F\ n = nP n+ i so 



Consequently the right side of (|2.72| ) becomes log[ P ^_^ A ^ ] and for \i 



(2.73) 



A with P = d\P we have 



logP(A) = J2 A 



-m—n r mn 

mn 



E E 

3=1 \n+m=j 



J? 

1 71 



Then using the elementary Schur polynomial defined in ( 2.5q ) and (2.67), we obtain 

oo oo 



(2.74) 



E 



(Zi = 0) 



m+n—i 

Thus we obtain the dispersionless Hirota equations, 



^ij — Xj+i(Zi — 0, ^2> • • • > Zj+i) 



(2.75) 



(2.76) 



These can be also derived directly from ( 2.62 ) w ith (2.70 ) in the limit e — > or by expanding ( 2.74 ) 
in powers of A~™ as in 124, 3l}|). The equations ( 2.76 ) then characterize dKP. 



It is also interesting to note that the dispersionless Hirota equations (2.76) can be regarded as 
algebraic equations for "symbols" F mn , which are defined via ( [2.71 ), i.e. 



B„ 



E 



F, 



(2.77) 



and in fact 



F nm — F mn — Resp [A m <iA™] 



(2.78) 



Thus for A, P given algebraically as in (2.67), with no a priori connection to dKP, and for B n defined 
as in ( [2.77 ) via a formal collection of symbols with two indices F mn , it follows that the dispersionless 
Hirota equations ( 2.76| ) are nothing but polynomial identities among F mn . In particular one has 
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from E4 



THEOREM 2.3. ( |2.78| ) with ( |2.76| ) completely characterizes and solves the dKP hierarchy. 



Now one very natural way of developing dKP begins with ( |2.67| ) and (2.41) since eventually 



the Pj+i can serve as universal coordinates (cf. here for a discussion of this in connection with 
topological field theory = TFT). This point of view is also natural in terms of developing a Hamilton- 
Jacobi theory involving ideas from the hodograph — Riemann invariant approach (cf. ^5|, |2^, |8C , 121 ] 
and in connecting NKdV ideas to TFT, strings, and quantum gravity. It is natural here to work with 
Q n := (l/n)B n and note that d n S = B n corresponds to d n P = dB n = ndQ n . In this connection one 
often uses different time variables, say T' n — nT n , so that d' n P = dQ n , and G mn — F mn /mn is used 
in place of F mn . Here however we will retain the T n notation with d n S = nQ n and d n P = n dQ n 
s ince one will be connecting a number of formulas to standard KP notation. Now given ( 2.67 ) and 
( [2.68 ) the equation d n P = ndQ n corresponds to Benney's moment eq uatio ns and is equivalent to 



a system of Hamiltonian equations defining the dKP hierarchy (cf. j80j, 121]); the Hamilton- Jacobi 
equations are d n S = nQ n with Hamiltonians nQ n (X, P = dS)). There is now an important formula 

(cf. ■ 



involving the functions Q 7 
by 



In particular one notes 



27, 121 1 ) , namely the generating function of dpQ n {\) is given 



1 



Pfo) - P(A) 



i 



d^i = d P Q n+ i{\) , 



(2.79) 



(2.80) 



which gives a key formula in the Hamilton-Jacobi method for the dKP [ 121 1 . Also note here that 
the function P(X) alone provides all the information necessary for the dKP theory. It is proved in 
§1 that 



THE ORE M 2.4. The kernel formula ( 2.79 ) is equivalent to the dispersionless differential Fay 
identity (|2.72p. 



The proof uses 

d P Qn=Xn-l{Ql,---,Qn-l) (2.81) 

where XniQii ■ ' • j Qn) can be expressed as a polynomial in Qi = P with the coefficients given by 
polynomials in the Pj+\. Indeed 



Xn = det 



p 


-1 











... o 


Pi 


P 


-1 








... o 


p 3 


Pi 


p 


-1 





... o 


Pn 


Pn-1 




Pi 


p. 


P2 P 



— dpQ n +i 



(2.82) 



and this leads to the observation that the F mn can be expressed as polynomials in Pj+\ = Fij/j- 
Thus the dHirota equations can be solved totally algebraically via F mn = & m n(P2, P3, • • • , Pm+n) 
where $ m n is a polynomial in the Pj+\ so the F\ n = nP n+ \ are generating elements for the F mn , 
and serve as universal coordinates. Indeed formulas such as (2.82) and ( [2.81 ) indicate that in fact 
dKP theory can be characterized using only elementary Schur polynomials since these provide all 
the information necessary for the kernel (2.79) or equivalently for the dispersionless differential 
Fay identity. This amounts also to observing that in the passage from KP to dKP only certain 
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Schur polynomials survive the limiting process e — * 0. Such terms involve second derivatives of 
F and these may be characterized in terms of Young diagrams with only vertical or horizontal 
boxes. This is also related to the explicit form of the hodograph transformation where one needs 
only dpQ n = Xn-i(Qi> ' ' ' ;Qn-i) an d the Pj+i in the expansion of P (cf. p4|). Given KP and 
dKP theory we can now discuss nKdV or dnKdV easily although many special aspects of nKdV for 
example are not visible in KP. In particular for the Fij one will have F n j — Fj n — for dnKdV. We 
note also (cf. j27, 121 ) that from ( 2.8l| ) one has a formula of Kodama 



P(jj) - P(A) 

3 ISOMONODROMY PROBLEMS 



(2.83) 



We begin with |192| where the goal is to exhibit the relations between SW theory and Whitham dy- 
namics via isomonodromy deformations. One considers algebraic curves C — > Co (spectral covering) 
over the quantum moduli space of z £ Cq ~ algebraic curve of genus or 1 via (3A) det(w — L(z)) = 
0. For the N = 2, SU(s) Yang-Mills (YM) theory without matter C = CP 1 and z = log(h) for 
h E CP 1 /{0,oo} with Lax matrix 



L{z) 



I h 


i • 





Cl 


b 2 ■ 








o • 








• 


• b a -i 


V h 


• 


■ c s _i 



Csh- 1 \ 



1 

Ik 



(3.1) 



Here bj = qj = dqj/dt with cj — exp(qj — qj+i) and one compares with 113 ] where q s +i = qi and 
qj — > — qj ~ Cj — * cj 1 . The eigenvalue equation (3A) becomes (A 2s = Y[l Cj = exp(qi — q s +i) 

1 s 
h + K 2s h- 1 = P(w); y 2 = P 2 (w) - A 2s ; h=-(y + P(w)); P(w) = w s + ^ UjW s ~ j (3.2) 



(where A 2s = 1 in [pL 13| ) . Thus one has hyperelliptic curves as in the case of a finite periodic Toda 
chain where such curves arise via a commutative set of isospectral flows in Lax representation 

d n L(z) = [P n {z),L{z)\; [d n - P n , 8 m - P m ] = (3.3) 

The P n are suitable matrix valued meromorphic functions on Co . The associated linear problem 

w^p = Lip; d n i) = P„V (3.4) 



then determines a vector (or matrix) Baker- Akhiezer (BA) function. Now in [ 192 1 one replaces such 
an isospectral problem by an isomonodromy problem 



oz 



(3.5) 



The idea is that the t flows leave the monodromy data of the z equation intact if and only if 

[d n - P n , ed z - Q] = 0; [d n - P n , d m - P m ] = (3.6) 
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An idea from [j75| is to write $ in a WKB form (e — * 0) 

'«2o\ -1/2 




exp[e 



"^(z)] 



(3.7) 



where d$ = wdz corresponds to the SW differential. Thus in the leading order term for ed z ^! = Q 1 ^ 
one finds (3B) (d z S)4> = Q4> and if we identify (3C) w = S z = dS = wdz then the algebraic formu- 
lation gives essentially the same eigenvalue problem as (3A) (where the relation between <fi and ip is 
clarified below). The idea from jr5| is that isomonodromic deformations in WKB approximation look 
like modulation of isospectral deformations. The passage isospectral — > isomonodromy is achieved 
via (3D) w —>■ ed z which is a kind of quantization with e ~ h and ^ ~ a quantum mechanical wave 
function. 

Now one separates the isomonodromy problem into a combination of fast (isospectral - t n ) and 
slow (Whitham -T n ) dynamics via multiscale analysis. The variables are connected via T n = et n 
and one assumes all fields u a — u a (t,T). Then writing d n ~ d/dt n we have (SE) d n u a (t,et) = 
d n u a (t,T) + e(du a (t,T) /dT n )\-T=et- Further take P n and Q as functions of (t,T,z) and look for a 
wave function as in (3.7) with <p n = <p> n (t,T,z) (say <po = <ft) and S = S(T,z). The leading order 
terms in (3.4) become, for (3F) w = d z S and ip — (p(z)exp(Y^t n (dS/dT n ), 



wtp = Qip\ d n ip = P n ip 



(3.8) 



which is an isospectral problem with associated curve C defined by det(w — Q(t, T, z)) = which is 
typically t independent but now depends on the T n as adabatic parameters. One can then produce 
a standard BA function tp — (f>exp(^2t n £l n ) (cf. |35|, [5^] and Section 2) where Q n — J z d£l n . The 
amplitude <fi is then composed of theta functions of the form 6(Y^ t n a n + ■ ■ •) where (3G) a n — 



(cr„ 



., li9 ) T with <r nj 



{l/2ni)§ B dfl n (we assume a suitable homology basis (Aj,Bj) has 



been chosen). The theta functions provide the quasi-periodic (fast) dynamics which is eventually 
averaged out and does not contribute to the Whitham (slow) dynamics. The main contribution to 
the Whitham dynamics arises by matching ip ~ ip and <j) ~ <f> with (3H) dS/dT n = Vl n (z) which 
in fact provides a definition of a Whitham system via (31) ddfl n /dT m = dd£L m jdT n representing a 
dynamical system on the moduli space of spectral curves. Thus one starts with monodromy and the 
WKB leading order term is made isospectral (via (3F)) which leads to a t independent algebraic 
curve. Then the BA fun ction corresponding to this curve is subject to averaging and matching to 
ip. We ignore here (as in [192Q all complications due to Stokes multipliers etc. The object now is to 
rel ate Q and L a nd t his is done with the N = 2, SU(s) YM example. The details are spelled out 
in HI (cf. also |||). 



We go next to [ 194 1 where isospectrality and isomonodromy are considered in the context of 
Schlesinger equations. Thus in the small e limit solutions of the isomonodromy problem are expected 
to behave as slowly modulated finite gap solutions of an isospectral problem. The modulation is 



caused by slow deformation of the spectral curve of the finite gap solution. Now from [ 194 1 let 
gl(r, C) N ~ ®fgl(r, C) - iV-tuples {A l7 ■ ■ ■ ,A N ) ofrxr matrices. There is a Gi(2, C) coadjoint 
action Ai — > gAig~ x and the Schlesinger equation (SE) is 



dti 



AiAl 



A, 



A, 



tk 



(3.9) 



and each coadjoint orbit Oi is invariant under the t flows. Thus (SE) is a collection of non- 



autonomous dynamical systems on Y\.\ ®i- We consider only semisimple (ss) orbits labeled by the 
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eigenvalues 6i a (a = 1, • • • , r) of A- Thus these eiqenvalues (and in general the Jordan canonical 
form) are invariants of (SE). There are also extra invariants in the form of the matrix elements of 



= — Ai, invariant via dA^jdti = 0. Assuming Aoo is ss it can be diagonalized in advance 
by a constant gauge transformation A — > CAiC~ Y and then only the eigenvalues 9oo,a (a = 1, • • • , r) 
of Aoo are nontrivial invariants. One can introduce a Poisson structure on gl(r, C) N via 

(3 J) {Ai a p, Aj pa } = Sij (—Sf3pAi i0ia + 5 aa Ai P ij) 

which in each component of the direct sum is the ordinary Kostant-Kirillov bracket. The (SE) can 
then be written as 



dAj 



1. 



{Aj,Hi}; H i =Res x = t .-TrM{X) 



(3.10) 



with {Hi,Hj} = 0. 



Now (SE) gives isomonodromy deformations of (3K) dY/dX = M(X)Y with M(A) = X)f [A/(A- 
ti)]. One can assume for simplicity that (3L) The A and A^ are diagonalizable and 9i a — 9i@ £ Z if 
a ^ (3. This means that local solutions at the singular points A = tx, ■ ■ ■ , tjv, oo do not develop log- 
arithmic terms. The isomonodromy deformations are generated by (3M) dY/dU = — [Ai/(X — ti))Y 
and the Frobenius integrability conditions are 



= 0; 



A 







A — U dti 



X-tj 



= 



(3-11) 



and these are equivalent to (SE). Next, since A = tx, 
there will be local solutions 



t/v,oo are regular singular points of (3K) 

oo 

Y i =Y i .(X-t i ) e *; % — ^2 Yin(X — U) n \ (3.12) 



-e 



where the Yi n are r x r matrices, l^o and Yoo,o are invertible, and ©i, are diagonal matrices of 
local monodromy exponents. This leads to (3N) A = YioQiY^ for i = 1, • • • , N, oo where 6^ = 
diag(9n, • • • , 6i r ). The tau function for (SE) can be defined in two equivalent ways (30) dlogr = 
'^2 1 Hidti or dlogr /d ti = Tr(QiY iQ 1 Yn and the integrability condition dHi/dtj — dHj/dti is 
ensured by (SE) (see | fl94| for details). The spectral curve i s no w (3P) det(M(X) — fil) = and 



this generally varies under isomonodromic deformations (sec [194 



Consider now Garnier's autonomous analogue of (SE). This is given by 



d_ 



A 



A 



-,M(A) 



= 0; 



d_ 



Ai 



A 



d_ 

dti 



Ai 



X 



= 



(3.13) 



and this is an isospectral problem with det(M (X) — fil) independent of t. An auxiliary linear problem 
is given by 

= M(A)V; |f - ~T~ — ^ ( 3 - 14 ) 

dU X — Ci 
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where ip is a column vector. This kind of problem can be mapped to linear flows on the Jacobian 
variety of the spectral curve and the case of r = 2 is particularly interesting since here the Painleve 
VI (and Garnier's multivariable version) emerge. 



For the geometry one goes now to |2j, |l4j, |95| |9j|. Let Co be the spectral curve (3Q) F(X, fi) = 
det(M (A) — fil) = which can be thought of as a ramified cover of the punctured Riemann sphere 7r : 
Co — > CP 1 /{ci, ■ ■ ■ , cjv, 00} where 7r(A, (i) = A. Then generically 7r -1 (A) ~ {(A, fj, a ), a = 1, • • • , r} 
and the /i Q are eigenvalues of M(X). Near A = c, one has /i Q = 9i a /(X — Cj)+ nonsingular terms, 
where the are the eigenvalues of A,-. Similarly near A = 00 one has /x a = —Ooo^X^ 1 + 0(X~ 2 ) 
where the #00, a are eigenvalues of A^,. One can compactify Co by adding points over the punctures. 
Thus near Cj replace (1 by p, = /(A)/i where /(A) = — c «)- Then the spectral curve equation 

becomes (3R) F(A,/2) = det(f(X)M(X) - jjtl) = and 7r _1 (A) = {(A,/2 a ), a = 1, ■•■,r} where 
(3S) /}„ = f'(ci)0i a + 0(X — Ci). Since the #j a , (a = 1, • • • , r) are pairwise distinct one is adding to 
Co at Ci, r extra points (A, p, a ) — (cj, f'(ci)6i a ) to fill the holes above Ci. At 00 one uses ft, — A/z with 
7r _1 (oo) = {(00, — #oc,a)}- Thus by adding rN + r points to Co one gets a compactification C of Co 
with a unique extension of the covering map 7r to a ramified 7r : C — > CP 1 . The genus of C will 
be ff = (l/2)(r - l)(rJV - r- 2). Taking (3T) M°(A) = Ei [^°/(A - <*)], A° = d*a 5 (0 u , ■ • • Ar) 
as a reference point on the coadjoint orbit of M(A) one compares the characteristic polynomials via 
(3U) F{X,(i)-F {X,fl) = /(A)EL 2 w(A)M r ^ where Pi (X) = Eo' h m iX m with & = {N-l)l-N, 
This leads to 



(3.15) 



1=2 



The spectral curve has therefore three sets of parameters (3V) (1) Pole positions Cj (i — 1, • • • , TV) (2) 
Coadjoint orbit invariants (i = 1, • • • , N, 00), and (3) Isospectral invariants h m g {I = 2, • • • , r; m = 
0,- • • - 1). 

Now one develops the Whitham ideas as follows. First reformulate the (SE) equation via (cf. 

d A, 



|9§) 



at* 



~&Ti dX OA' X-U A-Tj 
Note this corresponds to eij = Tj with <9i = ed/dTi and eA = A with d\ = ed\ and Ai/(X — t 
Ai/(A-Ti) so A - eA z . Then |Tf) becomes 



(3.16) 



(3.17) 



(the notation in [194| is confusing here since the step Aj — + Aj — eAj was not clarified). The 
auxiliary linear problem is (cf. (3K) and (3M)) 



dy ay Aj 



(3.18) 



(thus M — > M. and Y — > y). Then the e dependent (SE) equation ( J3 . 1 7| ) can be reproduced from 
the Frobenius integrability condition 



8 A i KA(k\ 8 



Ai d 



«9A 



d 



dTi A-Tj 9Tj A-Tj 



= 0; 



= 



(3.19) 
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Now start with (3. 17)- (3. 19) and introduce first variables ti = e 1 Ti so that (3.17) can be written 
as 

— (.1 — Oij)- 



Ti - Ti 



8 v ■ 



13 £-< T t - T k 



(3.20) 



In the scale of the fast variables the Tj can be regarded as approximately constant and hence (3.2C) 
looks approximately like Garnier's autonomous isospectral problem. Of course the T, vary slowly as 
do the h m £ (but not the 9ik)- 

For multiscale analysis following [[IlJ one writes now Ai — A4 (t, T, e) and ed/dTi — > d/dti+ed/dTi 
in ( 3.20| ) (note and Tj are considered to be independent now with ti = e -1 ?^ imposed only in the 

left side of the differential equation (|3.20D). Now assume (3W) Ai = A°{t, T) + eA}(t, T)-\ . The 

lowest order term in e gives 



dU ~ { %3 Ti — Ti 



Ti - T k 



(3.21) 



For Ti = Ci this is Garnier's autonomous system. The next order equation is 



dti 



(1 



[^4, , Aj] + [^4, , ^4^] 
Ti - Ti 



(3.22) 



_ §i ^[AlAl} + [Al,A° k ] 



Ti-T k 



where S refers to terms in A^ and their T derivatives. 

A standard procedure in multiscale analysis is to eliminate the t derivatives of A} by aver- 
aging over the t space but this involves some technical difficulties. The spectral curve for (SE) 
is hyperelliptic only if r = 2 and even in that situation there are problems. So another ap- 
proach is adopted. In addition to the multiscale expression for the Ai one assumes (3X) y = 

[<f>°(t, T, A) + <p 1 (t,T,A)e H }exp[e~ 1 S(T, A)] (y and the <p k are vector valued functions with S a 

scalar). Then one writes the auxiliary linear problem (3.18) in the multiscale form 



ay 



M(A)y; 



d_ 

dti 



_d_ 

l dTi 



y = - 



Ai 



A-T, 



y 



(3.23) 



The leading order term reproduces the auxiliar linear problem of the isospectral problem 



dS_ 

OA' 



M 



j.0. 



dti 



+ 



dS 

df~: 



-4 



A-T 



(3.24) 

t>°exp[Y,ti(dS/dTi)] so 
(3.25) 

Next one compares this ip with the ordinary BA function tp = cf>exp[^2 ijfii] where flj — j A '^ dfli 
and the vector function <f> involves 9 functions as before. For matching we want now (3Z) <fr = cj> 



where M° = E^U ?/(A - Ti). Define now (3Y) (j, = dS/dA and %p 
(3.24) becomes (3.14) in the form 
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(or more generally (j>° = h(T,A)<f>) and dS/dTi = f2j. Thus one proposes (•) fj, = dS/dA and 
dS/dTi = fli (and consequently dili/dTj — dftj/dTi) as the modulation equations governing the 
slow dynamics of the spectral curve (••) det(A4° (A) — /i/) = 0. Details in this direction are given 
in Section 6 of ]194| where it is determined that 

9 ' 

dQi — — —y^dh + nonsingular terms (3.26) 

near P a (A) 6 7r _1 (A) (a = 1,- • ■ ,r). Normalization is achieved via (• • •) § A dSli = for j = 
1, ' ' ' ,9- 

Next consider the hi for / = 1, • • • , g as the set of isospectral invariants h m l (cf. ( |3.15| )) and let 
S = S(T, h, A) with dS = fi(T, h, A)dA. One defines 

ddS dfi(T,h,A) JA 1 dfi 

dm = ^h~r dhj dA = W)dhi ( 7) 

Then (*) [dF{A,Jx)/dh m l] = f(A)A e fl r ~ m and there results 

(*) duw - -[A £ if- m /(dF/dp,)]dA 

These differentials dwj form a basis of holomorphic differentials on the spectral curve and one shows 
that 

;),, . /' ; i,/C r 

(3.28) 



daj i dds r 



is an invertible matrix (dQi = Aijdtuj where the dtuj are determined via § A dtuj = Su). Thus 
one has an invertible map h — ► a where the aj are the standard action variables of period integrals 
a j = § Aj dS. One can also express dS now as dS(T,a) with ddS/dTi = dfli and ddS/dai — dwi 
since 

ddS \ dhj ddS v-^, . . i ,„ „„. 

One can also introduce a prepotential as in SW theory with 



(cf. (3.3C) for Hi). 



4 JMMS EQUATIONS 



We go now to |193||, which is partly a rehash of 1 192 for the Jimbo-Miwa-Mori-Sato (JMMS) equa- 



tions, giving isomonodromic deformations of the matrix system (r x r) 
dY N A r 

— = M(X)Y; M(A) = u + £ u = £ u a E a (4.1) 



where E a ~ 5 a fj5 ai in the (/?, 7) position so u = diag{u\, ■ • ■ , u r ) where the u a ~ time variables of 
isomonodromic deformations. Evidently A = t j are regular singular points and A — 00 is an irregular 
singular point of Poincare rank 1. In order to avoid logarithmic terms one assumes the eigenvalues 
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6i a of Ai have no integer difference and that the Ui are pairwise distinct. Then isomonodromic 
deformations are generated by 



&Y 



A ■ BY 
T^tY; = (XE a + B a )Y- 

A — U OU a 



E 

p^a 



EaAooEp + EfjAozEa 



N 



A, 



with Frobenius integrability conditions 

d_ 



U a - U/3 



Ai 



X — tidX 



M(X) 



This leads to 



XE a — B a , — — M(X) 

OU a OX 



~q^~ = [tjE a + B a ,Aj]; 



= 



f^i - (i 5 \ l A i> A A I x 



+ Er^ 



t k 



(4.2) 



(4.3) 



(4.4) 



There are two sets of invariants (4A) Eigenvalues Gj a of Aj : (dOjp/dti) — (d8jp/du a ) = and 
(4B) Diagonal elements of : (dAoo^p/dti) = (cMoo,/3/3/ 'du a ) = 0. The JMMS equation is a 
nonautonomous dynamical system on the direct product FJ Q{ of coadjoint orbits and can be written 
in the Hamiltonian form 



9Aj -{A H V dA * 



= {Aj, Kot}] Hi = Tr\uAi + J2j^ 



N 



K a = Tr\E a ^2 UAi + £ 

1 P^a 



EaAocEpAoc 



Uf3 



(4.5) 



The dual isomonodromic problem can be formulated in a general setting where we suppose rank 
Ai = £i (which is also a coadjoint orbit invariant and constant under the JMMS equation in the 
generalized sense). Thus write first with r > li (cf. 95 96 , for more details and expanded 
frameworks) 

M(X) — u~ G T (XI — T)~ 1 F (4.6) 

where F = (F aa ) and G — (G aa ) are £ x r matrices, £ = Yl^i> an d T is a diagonal matrix of the 
form (4C) T = £f *iA with £>< = E ll+ ... +u _ 1+1 + ■ ■ ■ + E^ 
D% = Ee t + • ■ • + Ei, etc.) In particular Ai can be written as Ai 
understood as canonical coordinates with Poisson bracket 



k (so £»! = E x H h E t with 

-G T D l F and F aa , G aa may be 



{F aa , Ftp} = {G aa ,Gbp} = 0; {F aa ,Gbp} — S a b5 a f: 



(4.7) 



As an example consider li = 3 and £ 2 = 2 which implies £ — 5 and take r — 3 so F, G are 5x3 
matrices with D x = E x -\ h E 5 = I 5 and D 2 = E 4 + E 5 . Then A\ = ~G T I 5 F is 3 x 3 of rank 



23 



3 and A2 = —G T D2F is 3 x 3 of rank 2. The map (F, G) — + (A\, ■ ■ ■ , An) then becomes a Poisson 
map and the JMMS equation can be derived from a Hamiltonian system in the (F, G) space with 
the same Hamiltonians Hi and K a . The dual problem is formulated in terms of the rational I x i 
matrix (4D) L(y) = T — F(pl - u^G 7 where L(/x) = T + EI-Pa/G" - u a)} with P a = -FE a G T 
of rank 1 (but this restriction can also be relaxed). The map (F,G) — > (P±,P2, • • •) is Poisson for 
the Kostant-Kirillov bracket (4E) {Po,,a6, f/3,cd} = &up(~5bcPj3,ad + Sda,P/3,cb) and one can write the 
dual isomonodromy problem in the form 



— = -LU)Z- — 
dX ' du a 



-^—Z; ^- = -( f iD i + Q i )Z 

(1-U a Oti 



H tj 
3^£i J 



Now introduce a small parameter e via (cf. |l9q|) 



d_ 

8X 



d d 



d 



d 



dA' dU dT^ du a 
with u a — * U a E a and the JMMS equatins take the form 



dU a 



(4.8) 



(4.9) 



4£ = (i-<^)!^4U.v 



v ^ Ti — T 



a4j 



E 



U a -U p 



(4.10) 



(where .4, = eAj, = eij, U a = eu a , and Sq, = eB a ). Working as before we assume Aj 
Aj(t,u,T,U) and then, using ti — e~ x Ti and u a = t~ l Y ai one has e.g. 



dA, dA j dA, 
e — > + e - 

as in (^20|)-(^2il), leading to (Aj = A] + A)e + ■ ■ ■) 

<U; [TjE a + B a ,A°j] 



dU a 



(4.11) 



(4.12) 



% " (i j Ti — Tj 
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Ti - T fc ' 



where is given by the same formula as B a but with Aoa replaced by ^4^, = — -^H ■ The slow 
variables are parameters at the lowest order and to determine the slow dynamics one goes to the 
next order and one uses here the following approach. 

The lowest order equations are again an isospectral problem with Lax representation 



d 



dti A - T 



A 

3 ,M°(A) 



N 



0; M°(A) = U + J2 



A" 



A - T 
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^--AE a -B° a ,M° 
ou a 



(4.13) 



and det(fil — M°(A)) = E is constant under (t, u) flows. H depends on (T, U) however and the slow 
dynamics may be described as slow deformations of the characteristic polynomial of the spectral 
curve, namely H = 0. This spectral curve Co on the (A, fi) plane can be compactified to a nonsingular 
curve C; the projection n(A,fj,) — > A : Co — > CP 1 /{Ti, • • • , TV, oo} extends to C to give an r 
fold ramified covering of C over CP 1 . One adds points to the holes over A = (Tj,oo) as before. 
There will be r points (oo, U a ) (a = 1, • • • , r) over A = oo. For the 7$ write (4F) jl = / (A)/i with 
/(A) = nf (A-Ti) so S becomes (4G) F(A,/I) = det(p,I - f(A)M°(A)) = 0. Add then the r points 
(A, p) = (Ti, f'(Ti)6 ia ) (a = 1, • • • ,r) over T, to obtain a curve of genus 5 = (l/2)(r - l)(rAT - 2). 
To describe the slow dynamics one wants a suitable system of moduli in the space of permissible 
curves and the slow dynamics involves differential equations for such moduli. Now one can write 

r Ss-1 

F(A, fi) = F°(A, jl) + /(A) E h msA m fl s (4.14) 

s=2 m=0 

where (4H) S s = (N — l)s — N. F (A, jl) is a polynomial whose coefficients are determined by 
the isomonodromy invariants and (T,U). The coefficients h$ s _i tS (2 < s < r) are also determined 
by these quantities and the remaining coefficients /i TOiS for 2 < s < r, < m < S s — 2 give the 
appropriate moduli of number J2l=2^ s ~ 1) = 9- 

Thus the .4° are solutions of the isospectral problem with slowly varrying spectral invariants 
h m ,s(T, U) and to derive the modulation equations one writes 

Y = (<p°(t, u, T, U, A) + (j^e H ) x exp^SiT, U, A)] (4.15) 

The lowest order equations are 

9A^ V ' ati 9Ti A-Ti 

^ + || / = (A£ „ + B . W » (4 , 6) 

which determines 0° up to a multiplier — > (f>°h(T, U). These equations can be rewritten into the 
isospectral linear problem 



as , ,„ dS 

/i = _ ; ^ = ^^_ + £ Ua _ 



In particular the characteristic equation det(fj,I — A4 a (A)) = is satisfied by fi — d\S. One now 
identifies this ip with a corresponding algebro-geometric BA function via 

^ = ^exp^t^T, +22u a n Uoi j ■, n L = dn L (4.18) 

One must choose here a symplectic homology basis A/ , £?/ so that the WKB solution gives a correct 
approximation to the isospectral problem and we assume that this has been done (the problem 
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is nontrivial). The meromorphic differentials are characterized now by the following conditions: 
(41) dflxi has poles at the r points in 7r _1 (Ti) and is holomorphic outside of 7r -1 (Xi). Its singular 
behavior at these points is such that dQ^ = —dfi+ nonsingular. (4 J) d£Ljj a has poles at the r 



points in it (oo) with singular behavior dQu a — dA+ nonsingular. (4K) 



dVLrj =0. 



Matching the exponential parts in the two forms of if) gives then the modulation equations 

dS 



3T t Ti ' dU a 



= SI, 



These can be rewritten as 



ddS 



= dn 



ddS 



A=c 



dU a 



— dSljj a ; dS = /idA 



A=c 



leading to "classical" Whitham equations of the form (evaluation at A = c) 

ddfl T ddfl Tz 



ddfl 



Up 



ddflu a 



ddflu a 



(4.19) 



(4.20) 



(4.21) 



ddfl 7 
~~dU~ a 



dU a dUp ' dTi 
A similar multiscale analysis applies to the dual isomonodromic problem with WKB ansatz 

Z = ( X °(t, u, T, U, /i) + x X e + • • •) x expie-^T, U, fx)] (4.22) 

This leads to slow dynamics of the dual isospectral problem based on the dual expressions (4L) det(AI- 

= of the same spectral curve (cf. The modulation equations are obtained in the dual 
form (evaluation at fi — c) 



ddS 



dn Tt : 



MX 
dU a 



— dflun', dY> = — Adfi 



(4.23) 



(cf. 



llj). 



As before one can relate all this to SW ideas as follows: (A) One has a g dimensional period map 
h = (h ms ) — * a = (ai) where ai = § A dS for I = 1, • • • , g. One can show that the Jacobian of this 
map is nonvanishing and one may write h — h(T, U, a) (inverse period map). For a fixed this gives 
a family of deformations of the spectral curve producing a solution of the modulation equations; 
varying a one has a general solution. (B) One will have (for suitable holomorphic differentials duij) 



ddS 
dai 



= dcoi; (p dioj = 8 



(C) There is a prepotential F(T, U, a) satisfying 

— - b 
daj 1 



dS; 



13 



(4.24) 



(4.25) 



dT 
dU~ a 



d 2 T 
da j da j 



(where Hi and K a are as before). 
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5 GAUDIN MODEL AND KZ EQUATIONS 



We go to |19^ 

Zamolodchikov (KZ) equations (see e.g 



now in order to exhibit connections of isomonodromy problems and the Knizhnik- 

©0 



also jfL5|, [PR)). Here |l9|] follo ws ]130[ at first and later ]11 
3 and 4 it is appropriate to follow ]19 



114, 119, 124, 130, 134, n79], 184, 196| and cf. 



given the background in Sections 
here instead of the Russian school (to which we go later 
via [134, 135, 175, |1 76[ ) . First a generalized Gaudin mod el is based on a generalization of the XYZ 
Gaudin model to an SU(n) spin system following 13C]. Let X be a torus with modulus r, 
X = C/(Z + Zr). For integer (a, b) define 



i.e. 



°[ab] {z) 

exp[ir 



kf 



: 0,t); 

f 2i:i(m - 



(5.1) 



k)[z + k')] 



ra£Z 



Define nx n matrices (5A) J a b — g a h b and J ab — {l/n)J~ b where g — diag(l,Lu,- ■ ■ ,cj™ _1 ) with 
lu — exp(2iri / n) and h = (Si-ij); one notes that gh — uhg. The J matrices give a basis of su(n) 
over R and of sl(n, C) over C with Tr(J a bJ cd ) = S ac Sbd- As an R matrix define 



R(X) 



(o,l))6Z„xZ, 



W ab (\, V )J ab ® J ab ; W ab 



9[ab] (A + rf) 

{ab](v) 



(5.2) 



These Boltzman weights are n-periodic so the summation is as indicated. The r-matrix is the leading 
part in the r\ expansion of R at r\ = 0, i.e. 



r(A) = ^ W ab Jab ® J ah ; Wab{z) 
(a,i>)#(0,0) 



I [H(A)^oq ] (0) 
'[a6] (0)^00] (A) 



(5.3) 



where ' ~ d/dz. The sum is now over (Z„ x Z„)/{(0,0)} (alternatively one could take wqq(z) = 0) 
The r matrix satisfies the classical Yang-Baxter (YB) equation 

[r (13) (A))r (23) (/i)] 



[r( 12 )(A- M ),r (13) (A)+r( 23 )(/x)] 



(5.4) 



where e.g. r( 12 )(A) = ^w a b{X)J a b®J ab ®I- The generalized Gaudin model is a limit as 77 — > of an 
inhomogeneous SU (n) spin chain with TV lattice sites and the monodromy matrix of the spin chain 
is (5B) T(A) = Ljv(A — tjv) • • • ii(A — ti). The U are inhomogeneous parameters which eventually 
correspond to time variables in an isomonodromy problem and (5G) L^(A) = J2( a b) Wab(A)</ a & <8> 
Pi(J ab ) acting on C n ®Vi where Vi is the representation space of an irreducible representation (pi, Vi) 
of su{n). The L operators satisfy (5D) RLL = LLR and the monodromy matrix acts nontrivially 
on C n (g) V for V = V\ ® • • • §3 Vat. The leading nontrivial part in the 77-expansion of T(X) at 77 = is 
(5E) T(A) = J^i J2(ab) u 'afc(A — ii) J a b <8) pi(J ab ) and one has the fundamental commutation relation 
(5F) [T 1 ( A) , T 2 (/i)] = - [r ( A - M ) , T 1 ( A) + T 2 (/x)] where e.g. T 1 ( A) = £ 2 £ (o6) Wafc ( A - t, ) J ab ® I <8> 
Pi{J ab ) (i-e. the monodromy acts nontrivially in the first C n of C n (g> C n (g> V). Mutually commuting 
Hamiltonians Hi (i = 0, 1, • • • , N) of the generalized Gaudin model are defined via 



N 

^-Tr nxn T(X) 2 = C{P{\ + #*C(A - U) + H 
1 



JY 



(5.5) 



where V{z) and ((z) are Weierstrass functions with modulus r. The Ci are quadratic Casimir 
elements of the algebra generated by pi(J ab ), i.e. one has (5G) d — (1/2) Yl( a b) Pi(Jab)pi{J ab ) and 
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(5H) Hi = E j¥i E(ab) w ab(U - t j )p i (J abPj (J ab ) for i > with J2x Hi = 



One considers now (51) M(A) = Y^,i X](ah) w afc(A — ti)J a bAf b as a classical analogue of the T 
operator in (5E); the Af b are scalar functions of U and r (for i = 1, • • • , N) and the isomonodromy 
problem will be formulated as differential equations for these functions. Note the Af b depend on t 
but the pi(J ab ) do not; this is related to the difference between Heisenberg and Schrodinger pictures 
(as indicated below). The commutation relations for pi(J ab ) can be packed into 



{M(A) ? M(/i)} = -[r(\-ii),M(\)®I + I®M{ji)] 



(5.6) 



Here the left side is an abbreviation for (5J) J2i M W> M cd (p)}J ab ® J cd where M ab are the coef- 
ficients in M = J] M ab J ab . In terms of the residue matrix (5K) A; = Res\ =u M{\) = J2(ab) J abAf 
the Poisson structure is simply the one induced from the Kirillov-Kostant Poisson structure on 
sl(n, C); symplectic leaves are the direct product 0\ x • • • x On of coadjoint orbits Oi in sl(n, C) 
on which Aj is living. These leaves become the phase spaces of the nonautonomous Hamilton 
system described as follows: Using M(X) for T in (|]^) with Cj Casimir elements in the Pois- 
son algebra above, one thinks of Hamiltonians Hj for the classical Gaudin system in the form 
(5L) H = J2j^iJ2(ab) w ab(ti — t j) A ,ab A\ b . i?o is also a quadratic form in the Af> but is more 
complicated. The t% here are just parameters of the classical Gaudin system, which is a Hitchin 
system on the punctured torus X/{t\, • • • , £jv}. The nonautonomous Hamiltonian system will be a 
nonautonomous analogue of this Hitchin system in the form 



dAf 



= {Af, HA; = Lf, JL (h - m J2 UHA I 



(5.7) 



where rji arises via the known transformation law (5M) C,(z + 1) = Q{z) + This nonautonomous 
Hamiltonian system is an isomonodromy problem which can be put in Lax form. The key is a 
relation 

Y Y W -*,-b(» - *)Wab(p - U)J ab Af, M(/i) (5.8) 
1 (ab) 



M{X)^TrM(p) 2 



(not proved - one can use function theoretic methods as with the YB equation). Taking residues 
gives 

{M(X),Hi} = -[A(A),M(A)]; A,(A) = Y w ab (X - t,) J ab Af (5.9) 

(ab) 

and after lengthy calculations 



N 



{M(X),H a } = [4 7 r l B(A)-7 ?1 ^^A(A),M(A)]; 



(5.10) 



JV 



B W = J2J2 z - b( - x - ^) J abAf; Z ab 

1 (ab) 

The first formula can be written in the form 



Wab (x) ( o{ ab] w _ e[ ab] (oy 

4jri \9 [ab] {X) 9 [ab] {0) 



N 



{M(A), H - m Y = [iiriB(X), M(A)] 



(5.11) 
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and one obtains the Lax form of the nonautonomous Hamiltonian system as 

dM{\) dA(X) 

-^— = -[ Ai (\),M(\)]--^-; (5.12) 

^ = [2 S (A),M(A) ]+2 ^ 
One notes that t derivatives correspond to A derivatives via 

(5N) (d9 [ab] (X)/dr) = (l/im)(d 2 [ab] (X)/dX 2 ) 
These equations are integrability conditions for a linear system 

^ - M(X)\ Y(X) = 0; (5.13) 
!? +A(A)V(A) = 0; ( |- - 25(A) )y(X) = 



dti J \dr 

where the last two equations can be interpreted as isomonodromy deformations of the first equation. 



Now the elliptic KZ equations of |13C] can be written as 

+J2 J2 Wab ^ - t j)P^( J ab)Pj(J ab ) ] F(t) = 0; (5.14) 

1 (ab) 

+ Y. Z - b ^ - t j)pi(Jab)p J {J ab ) ] F(t) = 

T i,j=l (ab) 

where k = k + n with k the level of a twisted WZW model. These equation s ch aracterize N-point 
conformal blocks with irreducible representations pi sitting at Following ]179[ one adds another 
marked point at A with fundamental representation (C n ,id) to obtain 

K-gr + Wa »( X ~ t i) J abP*(J ab ) ] G(X, t) = 0; (5.15) 

1 (ab) 

k-^- +J2J2 Wab ^' - t i)pi{ J ab)pj{ jab ) + ^2wab(ti - X) Pl {J ab )J ab ] G(A, t) = 0; 

1 &i (ab) (ab) 

d N N 

K fr+1212 Zab ^ - t j )Pi(Jab)p j {J ab ) + Z "b(ti - X) Pl (J ab )J ab 

i,j = l (ab) 1 (ab) 

N 

+ 5]5I^a h (A-t J )J afc p J (J ab ) + ^Z afc (0)J ab J afc ] G(X,t) = 

1 (ab) (ab) 
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Let F — exp(S) be a fundamental solution of ( |5.14 ) (cf. ]179{ ) and consider the resulting equations 
for X = F~ 1 G of the form 



N 



w ab {\ - u)j ab A 

1 (ab) 



ab 



X(\,t) = 0; 



(5.16) 



(ab) 



5 



TV 



*^ + 2 EE^( A 

1 (afc) 



b (A - U) JabA?" X(A,t) = 0; 



i t ) JafcAf + }J Z Qb (0) J afc J afc I X(X, f) = 

(ab) 



where (50) Af = F(t)- 1 p l {J ab )F{t). This is almost the same as fl5.13| ) except that: (A) The 
last equation contains an extra term. (B) The A^ b are not scalar functions but operators on 
V = V\® ■ ■ ■ ®Vn- (C) There is an arbitrary parameter k; the previous monodromy (i.e. ( p. IS )) 
can be reproduced by setting k = — 1. Now the first descrepency (A) can b e removed by gauging it 
away via X [exp(f (t))]X- As for (B) one remarks that the system ( 5.16 ) is a quantization of the 
isomonodromy problem ( 5.13 ); this is parallel to the relation between quantized and classical Hitchin 
systems. The operators Af b depend on (t, r) but inherit the same commutativity as the J matrices 
from pi(J ab ). The passage from pi(J ab ) to Af b amounts to a change from the Schrodinger picture 
to the Heisenberg picture and classical limit now means replacing these operators by functions on a 
phase space, thus giving the isomonodromy problem. 



6 ISOMONODROMY AND HITCHIN SYSTEMS 



We will extract here from [175| which gives an excellent survey of many matters (cf. also 114, 134 
13o| , 136, |176| ). However first some remarks about RS and VB are appropriate. 



REMARK 6.1. In a first version of this survey we made up a number of remarks (seven) on 



(seve: 

M 
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The 



we will assume various 
L — > E) is a 



connections, vector bundles (VB ) , and Hitchin systems follow i ng 1 10|, O, |4§|, [l9"|, |5(j|. |51 . 63 
lOOj [iOli PI [T03l [LOT], ^ |l2(], HI [l33, m, |l3|, |l36|, [Til], |148| , |156[ |167| , U$t 

result was complete enough but lacked coherence so we will try again using \ p7\ as a springboard. 
We will first simply indicate, generally without proof, some facts about RS, VB, and sheaves from 
[ iol (for background information cf. || |7j [n], [t|, |o| ~ 
definitions are known). First a holomorphic section of a line bundle L over E (ir 
holomorphic map s : S — * L such that it o s = ids and the space of such sections i?°(S, L) is finite 
dimensional. Given a point p G E, with [/ a neighborhood (nbh) of p, V = E/{p}, and z(p) = 0, 
then, using 2 as a transition function on U n V, one produces a line bundle L p which will be useful 
later. The canonical bundle K is the bundle of holomorphic 1-forms and on E = P 1 one defines 
0(n) as the line bundle with transition function z n on U n V where e.g. € U and V = P 1 /{0}); 
the dimension of iJ^P 1 , 0(n)) is n + 1. Given a line bundle L one forms L* = L^ 1 via transition 
functions g a b(L*) = 5 a ~ 1 (i) and there are obvious relations Hom{L,L) ~ L* ® L, etc. If E is com- 
pact then genus(T,) = g = dim i?°(E, K). 

One defines sheaf cohomology in the standard manner and we recall for a line bundle L, H 1 (E, L) ^ 
i/°(E, K ® L*) (Serre duality). The isomorphism classes of holomorphic line bundles on E are given 
by elements of _ff 1 (E,C*) (called the Picard group) and for given degree d, J d ~ line bundles of 
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degree d is a complex torus isomorphic to Jac(X) (more on this below). Note also for S = O(L) = 
sheaf of holomorphic sections of L, that iJ p (E,6>) = for p > 1; if S = C or Z then H P (E,S) = 

for p > 2. From the short exact sequence (•) — » Z — > — > 0* — > 1 with ~ exp{2irif) one 
determines a long exact sequence leading to 

° |rf^| HHS.O*) - ^ 2 (S,Z) = (6.1) 



where H 2 (Y,,Z) ~ Z, tf^E,©) ~ C 9 , and ff^Z) ~ Z 2 f. The coboundary operator <5 in Q) 
yields (since if 1 (S,C*) ~ equivalence classes of line bundles) = deg(L) — c%(L) £ Z (first 

Chern class). Thus for L such that <5([L]) = d one obtains J d (E) ~ Jac(E) One shows next 
that if a section s G H (S, L) vanishes at pi £ E with multiplicities rrij then deg(L) — ^ m, so 
if deg(L) < then L has no nontrivial holomorphic sections. For vector bundles V one defines 
deg(V) — deg(det(V)) — c\{V) where det(V) ~ A m V for rank(V) = m has transition functions 
det(guw) in an obvious notation. We recall also the Riemann-Roch (RR) theorem which states 

dimH°{Yl,V) - dimH l (Y,,V) = deg(V) + rank{V)(l - g) (6.2) 

Many nice results can be obtained from this. Another classical theorem states that any rank m 
holomorphic VB V — > P 1 has the form V ~ O(ai) © • • • © C(a m ) for a, S Z. One defines di- 
rect image sheaves via f*S(U) = S(f~ 1 (U)) for / : E — » E and for 5 = 0(£) there results 
(A) #°(E, /*£>(£)) ~ H (t,O(L)) with (B) f*0(L) = 0{E) for E a rank m holomorphic VB 
(where m = deg(f) — deg(f*L p ) = #[/ _1 (p)]) and (C) If V is a holomorphic VB on E then 
f*0(L © f*V) ~ C(-E © V) (where E is given in (B)). One obtains then for L and £ as indicated 
deg(E) — deg(L) + (1 — g) — deg(f)(l — g) (note here a useful fact in proving theorems is that for V a 
VB on E if°(E, VL~ n ) = for n sufficiently large or equivalently, on P 1 , V(n) — V ®0(n) will have 
holomorphic sections for n sufficiently large). Next note that if deg(E) = above then E is trivial if 
and only if L(— 1) — L® f*0{— 1) has no nontrivial sections. Concerning Jac(E) one usually defines 
J ff_1 (E) = Jac(E) and the the image of the holomorphic map (pi, • • • ,Pg-i) — > ■ • • ^p s _i (tensor 
product) from E x • • ■ x E — > J 9_1 (E) is called the O divisor; it is of codimension 1 in Jac(E). 

Now going to Lax equations and integrable systems, take a section w G H°(T,, f*0(n)) where 
/ : E — > P 1 and for U C P 1 multiplication by w defines a linear map w : H°(f^ 1 (U),L) — > 
i/°(/- 1 (C/), J L(n)) where L(n) = X © f*Q{n). One takes here L such that L(-l) e J^- 1 / and 
/»C?(L) = O(S) with E trivial (thus deg{E) = and L(-l) = L® f*0(-l) has no nontrivial 
holomorphic sections). By definition of E this w determines a homomorphism W : H a (U,E) — > 
H°(U,E(n)) and thence 

VF : H°(P^£;)~C w, -^fr (P 1 ,£7(n))=iC TO ®H (P 1 > O(n)) (6.3) 

Thus one has an m x m matrix valued holomorphic section of O(n) so W = A(z) = Aq + Aiz + 
■ • • + A n ^ n . This gives a construction L — > A(z) and from A(z)si = wsi on sections (via local 
constructions) one is led to the algebraic curve (••) S : det(wl — A(z)) — w m + ai(z)w m ^ 1 + ■ ■ ■ + 
a m (z) = which is an m-fold covering of P 1 . One can think (rather crudely) of w embedding E 
as a one dimensional submanifold S of the total space of 0(n) (c C 2 locally) and eigenspaces of 
A(z) correspond to line bundles on E or S. This is not clearly stated and the matter should be 
clarified below. As for Lax pair equations dA/dt = [A,B] with A(z), B(z) matrix polynomials one 
sees that the spectrum of A is preserved in such evolutions so the eigenspace bundle moves in a 
complex torus and if a line bundle follows a linear motion there is a basis for which A(z) evolves as 



a Lax pair for a specific B depending on A (cf. [104] for details). Finally one comes to completely 
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integrable Hamiltonian systems (CIHS). We recall that symplectic manifolds M 2N — M have a 
nondegenerate closed 2-form uj (and we assume uj is holomorphic for complex manifolds M). This 
provides an isomorphism TM <-> T*M so that for a given function H : M — > C (Hamiltonian) 
one has a vector field Xh characterized by ui(Xh,Y) = Y{H) for vector fields Y. The Poisson 
bracked is defined via {/,<?} = Xf(g) = —X g (f) and ^7/, 9 } = One speaks of coadjoint 

orbits of a Lie group G as follows. If £ € g* then write ui^(X, Y) = £([X, Y\) where X, Y E g define 
tangent vectors to the orbit at £. Further for f,g € C°°(g*) one can express € 5 as a linear 

form on Tgi?* = 5* and write {f,g}{Q =< W(£),dg(£)] >. One writes also Ad* : g* — + g* via 
< Ad*((;),X >=< £,Ad g -iX > (and similarly < ad* x £,Y >= - < £, [X, Y] >). The momentum 
map is an equivariant map /i : M — > g* and as an example of all this let G — SL(m, C) and 
use Tr(AB) to identify g and g*. Take a product of coadjoint orbits O x •• • x of elements 
of g ~ g* with distinct eigenvalues. The moment map is (Ri,- ■ ■ ,Rk) — *• with symplectic 

quotient M = {(Ri)/Yl Ri = 0} which has dimension (Jjt) 2N — k(m 2 — m) — 2(m 2 — 1). Let 

k k 

A{z) = p(z) J2 P(*) = f[(z ~ «*); (6.4) 

det(w - A(z)) = w m + a 2 (z)w m ~ 2 H + a m (z) = p(z,w) 

The eigenvalues Afc of A{au) — Rk are fixed in advance and p(ctj, w) = det{w — Rj) = for w = A 



so the polynomial coefficients a^(z) are not arbitrary. It is shown in [ 104 1 how Lax pair equations 
arise in this example upon fixing the a, subject to constraints (cf. also |2|). In any event one defines 
a CIHS to be a symplectic manifold M 2N with N Hamiltonian functions Hj such that \H j , Hj ) = 
and A?dHj # 0. The bigger picture developed in particular in @, ^, P, [lOOj ^Olj, |Io| |o| |oj , 



134 , 13J, 136, 14S, 175, 176, 182] can be illustrated now by considering i?(z) = X^I^i^ 2 /! 2 — a j 
with Rj the residues of a matrix valued differential on P 1 . The fact Ri = corresponds to the 
vanishing of the moment map for symplectic reduction. Now replace R(z) by a 1-form on an arbitrary 
compact RS with values in End{V) for an arbitrary holomorphic VB V. The numerical miracle now 
involved consists in obtaining in suitable circumstances precisely the right number N of Hamiltonian 
functions. Thus e.g. in the example above one requires the number of independent coefficients a, in 



(3.4) to be N (cf. (♦)). In order to obtain a suitable structure in general one considers the space 1Z 
of equivalence classes of stable holomorphic rank m vector bundles V over Stable has various 
definitions (which we omit) and the facts used here are only that the only global endomorphisms 
of V should be scalars and that 1Z is a complex manifold with dim(lZ) = 1 — m 2 (l — g). Then 
T*U is a symplectic manifold and since TTZ at [V] is H 1 (Y,,End(V)) Serre duality gives T*TZ ~ 
H°(E, End(V) ® K). Thus locally a point in T*1Z is (up to equivalence) a VB V and a holomorphic 
section A of End(V) ® K. Therefore A is an m x m marix with values in K and its characteristic 
polynomial is p(z, w) — det(w — A(z)) — w rn + aiw 711 ^ 1 + ■ ■ • + a m with holomorphic coefficients a, € 
H°(E, K l ). One can show that dim i?°(E, K) = g and dim H 1 ^, K) = 1 so by RR deg(K) = 2g-2 
(g > 1) and by Serre duality again H 1 (Y 1 ,K % )* ~ ff°(S, K x ~ l ). But i^ 1 " 1 has negative degree and 
no sections for % > 1 so RR for if 1 implies dim if 8 ) = (2z — 1) (g — 1) . The number of degrees of 

freedom in choosing the characteristic polynomial p(z, w) is then g+J2^ (2i~l)(g— 1) = 1— m (1— p) 
and this is exactly the dimension of 1Z or (1/2) x dim(T*lZ) . This will then give rise to a CIHS 
based on T*1Z and p(z,w) — determines a spectral curve S in the total space of K over E; A 
is determined by a line bundle on 5* and the Hamiltonian flows are linear. Again the construction 
needs clarification (see below). 



We go now to ]17q ] and let T, g be a RS of genus g and Fi?(E, G) be the space of fiat vector 
bundles (VB), V = Vg (G — GL(N, C)), with smooth connections A. Flatness means (6A) J- a — 
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dA + (1/2) [A, A] = 0. Fixing a complex structure on S g one has A ~ (A, A) with a system of matrix 
equations 

(<9 + A)V> = 0; (d + A)^^0 (6.5) 

(linearization of (6A)). One modifies this via a parameter k € R (the level) and uses k9 instead 
of d in the first equation of ( |6.5| ). Now let /i be a Beltrami differential ^ G f2^ 1,1 ^(S £? ) (cf. [ |l()9| ). 
Then in local coordinates (6B) fi — fi{z, z)d z ® dz and one can deform the complex structure on E fl 
to produce new coordinates 

w = z — e(z, z); w — z; fi — - — — (6-6) 



(note this should be called —fi to agree with dw/dw but we retain the notation of |175|| ). Then 
dm = d + fid annihilates dw while d annihilates dz (note (d + fid)(z — e) = —Be + fi — fJ.de = 0). In 
the new coordinates (6A) has the form (6C) Ta = (d + dfi)A — ndA + [A, A] — and one arrives 
at the system 

(nd + A)^ = 0; (8 + fid + A)ip = (6.7) 

Write now fi = XL=l ^Ao where fi\, ■ ■ ■ , fi® is a basis in the tangent space to the moduli space A4 g of 
complex structures on E g (here I = 3ff — 3 fo r g > 1). Note that fid measures the deviation of d^, from 



d. Fix now a fundamental solution of (6.7) via ^(zq,zq) = /. Let 7 be a homotopically nontrivial 
cycle in E g such that (zq,zq) € 7 and set (6D) 3^(7) = VK^o, ^o)l7 = P e xp§^A (monodromy - P~ 
path ordered product as in jllj] for example where many basic ideas about connections etc. are 
spelled out). The set of matrices ^(7) generates a representation of IIi(I] 9 ,zo) m GL(N,C) and 
independence of the monodromy J- 7 to deformation of complex structure means 

d a y = (a =!,-••, d a ^^A (6.8) 



It follows that ( |6.8D is consistent with (3.7) if and only if 



<9 a A = 0; 91= -Afi° a (a = 1, ■■■,£) (6.9) 

Further this system ( |6.9[ ) is Hamiltonian where one endows FB(Y>,G) with a symplectic form 
(6E) w° = J* E < (L4, <5A > where < , >~ Trace; Hamiltonians are defined via (6F) H a = (1/2) J s < 

A,A>/x2 (o=l,---,0. 

Consider now the bundle V over .M g with fiber FB(E,G); the triple (A, A,t) can be used as a 
local coordinate and one thinks of V as an extended phase space with a closed 2-form (6G) uo = 
u>° — (1/ ft) ^ Q 5H a St. Although w is degenerate on V it produces equations of motion (6.9) since oj° 



is nondegenerate along the fibers. Gauge transformations in the deformed complex structure have 
the form 

A - f^ndf + f- l Af; A -» /- x (5 + /id)/ + /- x 4f (6.10) 

The form u> is invariant under such gauge transformations (the set of which we call Q) but not w° 
or H a independently. If one writes now (6H) A 1 = A — (l/n)fiA and uses (A, A 1 ) as a connection 
then (61) uj = / s < <L4, <5A' >. 



A gauge fixing via (6.10) plus the flatness condition (6C) is in fact a symplectic reduction from 



the space of smooth connections to the moduli space of flat connections FB{Y>, G) — SM(T,, G)/ /Q 
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where / / means (6C) plus a gauge fixing are in force. The flatness condition is called the moment 
constraint equation. Now fix the gauge so that the A component of A becomes antiholomorphic, 
i.e. (6J) 8L = via L = f (8 + /J.d)f + f~ 1 Af. This can be achieved since (6J) amounts to the 
classical equations of motion for the WZW functional Swzw (f> A) with gauge field / in the external 
field A (cf. ptpO|). Let L be the gauge transformed A, i.e. L = f^ndf + f^Af, and then (6C) 
takes the form (6K) (d + /j,d)L + [L, L] = 0. Thus the moduli space of flat connections FB(T,, G) 
is characterized by the set of solutions of (6K) along with (6J) and this space has dimension 
(6L) dim(FB) — 2{N 2 — l)(g — 1) (for g > 1). After gauge fixing the bundle V over M g becomes 
V with FB as fibers and the equations (|6.7|) - (|6.8[) become 

(nd + L)ip = 0; (a + fid + L)ip = 0; {nd a + M a )ip = (6.11) 

where we have replaced ip by f~ lr tp and M a = nd a f ]~ l ■ Note that k is not involved in the 
WZW result and thus k seems to be arbitrary in (6.11). The gauge transformations do not spoil 
the consistency of this system and from (6.11) one arrives at the Lax form of the isomonodromy 
deformation equations 

d a L- K dM a + [M a ,L] = 0; (6.12) 

nd a L - f i°L = (a + f id)M a - [M a , L] 

These equations play the role of ( |6.9| ) and the last equation in ( |6.12 ) allows one to find M a in terms 
of the dynamical variables (L,L). 

The symplectic form u on V is 

w = / < SL, SL > -- V SHJt a ; H a = \ <5L,SL> fj° (6.13) 

and we introduce local coordinates (v, u) in FB via (6M) (L, L) = (L, L)(v, u, t) with v — (vi, ■ ■ ■ , vm) 
and u = (ui, ■ ■ ■ , um) for M — (N 2 — l)(g — 1). Assume for simplicity now that this leads to the 
canonical form on FB 

uj° = / < SL(v, u, t), 5L(v, u, t) >= (5v, Su) (6.14) 
where ( , ) is induced by the trace. On the extended phase space one has then 

ui = (Sv, Su) 5Ka(v, u, t)6t a (6.15) 

a 

and variations in the K a take the form 

SK a = [ [<L,5L> [i Q a + n(<5L,d a L>- <d a L,5L>)} (6.16) 

Now because of (6K) the Hamiltonians depend explicitly on times and one considers the Poincare- 
Cartan integral invariant (cf . || ) 

6 = S^uj = (v, 5u) - - V K a (v, u, t)8t a (6.17) 

a 

Then there exist 3g — 3 = dim(M. g ) vector fields (6N) V a = nd a + {H a , ■} (a = 1, ■ • • ,£) that 
annihilate Q and one can check that (K a ~ H a here and below) 

Kd s H r - nd r H s + {H s ,H r } ua = (6.18) 
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Note from pi that V a ~ ^2(Afd Vi + Bfd Ui ) + d a and V a G fcer(0) means (d a ~ <9/<9i a ) 



so that (6NN) V Q = — {dH a Jdui)d Vi + (dH a /dvi)d Ui +d a Thereby they define the flat connection in 
V and these conditions are sometimes referred to as the Whitham hierarchy. Thus for E G M. g the 
Whitham equations determine a flat connection in FB(E, G). This gives an entirely new perspective 
for the idea of Whitham equations; it is based on geometry and deformation theoretic ideas (no 
averaging). For a given f(v, u, t) on V the corresponding equations take the form 

df(v,U,t) df(v,U,t) .jrr fX , R on> 

-^^ = K ^I^ + {HsJ} (6 - 20) 
called the hierarchy of isomonodromic deformations (HID). Both hierarchies can be derived from vari- 
ations of a prepotential F on V where (60) F{u, t) — F(u$, to) + ^2 Ju t ^ s< ^'s with ^s{d s u, u, t) — 
(v, d s u) — K s (v, u, t) the Lagrangian (d s u = SK S /Sv - note F ~ S is better notation). F then satisfies 
the Hamilton-Jacobi (HJ) equation 

Kd s F + H s (^,u,tj =0 (6.21) 
and F ~ log{r) where r corresponds to a tau function of HID (the notation is clumsy however). 



Note that equations such as ( |6.2l| ) arise in |137| with period integrals aj as moduli (cf. also Section 
8 - recall that the Toda type Seiberg-Witten curves for massless SU(N) have N — 1 = g moduli 
with a 1 — 1 map to the a,). Singular curves are important for low genera situations but we omit 



this here (cf. 175 



Consider next t he modul i spac e 1Z = TZ 9 ^n of stable holomorphic GL(N,C) ve ctor bund l es V 
over E = E 9 (cf. fl0|, pcj, [l78| for stable and semistable VB and note that in UJ, 111= HH 



one develops the theory for E ~ E„ n with n marked points but this will be omitted here. 1Z 
is a smooth variety of dimension g = N 2 (g — 1) + 1. Let T*1Z be the cotangent bundle to 1Z 
with the standard symplectic form. Then Hitchin defined a completely integrable system on T*1Z 
(cf. |{100[ ), The space T*1Z can be obtained by a symplectic reduction from T*1Z s g N = {(cj>,A)} 
where A is a smooth connection of the stable bundle corresponding to d + A and is a Higgs 
field 4> G fi°(E, End(V) ® K) where K is the canonical bundle of E (K ~ holomorphic cotangent 
bundle - cf. pa]). One has a symplectic form (6P) w° = J s < 5(f), 8 A > which is invariant under 

the gauge group Q = C°°(E, GL(N, C)) where (6Q) <j> -► /"V/ and i -► + J" 1 !/ with 

ft gAr = n s g N /G. Let now Pskd^' 1 8) dz be (-fc + 1, 1) differentials (p afe G fl" 1 ^, Tk ~ 1 ® #) where 
r' " 1 ~ fc — 1 times diffcrentiable sections, and s enumerates the basis in iJ 1 (E,r fc_1 K) so 
Ps,2 ~ A*s)- By Riemann-Roch dim(H 1 (Yj : T k ~ 1 if)) = (2fc — l)(g — 1) and one can define gauge 
invariant Hamiltonians 

H s<k = i jT < fe > p s , fc (A = 1, • • • , N; s = 1, • • • , (2k - l)(g - 1)) (6.22) 
where the Hamiltonian equations are 



>%<> ■■ [d a = A ; a =(«,*;)] : =</ V'.s./. (0.23) 
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The gauge action produces a moment map \i : T*1Z S — > gl (N,C) and from (6P)-(6Q) one has 
(6R) fi = 5<fr + L4, </>]. The reduced phase space is the cotangent bundle T*lZ g ^ N ~ T*TZ S //Q = 
(0) /Q (gauge fixing and flatness) and finally the Hitchin hierarchy (HH) is the set of Hamiltonians 
( p^22| ) on T*K g , N . Note also that the number of H s , k is £f (2fc ~ l ){9 ~ l ) = N *(9 ~ !) + 
1 = 5 = dim{T*TZ g ^). Since they are independent and Poisson commuting (HH) is a set of 
completely integrable Hamiltonian systems on T*TZ gt N- To obtain equations of motion for (HH) fix 
the gauge of A via A = fdf^ 1 + fLf^ 1 so that L = f~ 1 4>f is a solution of the moment constraint 
equation (6S) dL + [L,L] = (flatness condition). The space of solutions of (6S) is isomorphic to 
-ff°(S, End(V) <E) K) ~ cotangent space to the moduli space H g ,N- The gauge term / defines the 
element M a = d a f / _1 E gl(N, C) while the equations 

(X + L)Y = 0; (d a + M a )Y = 0; (6.24) 
{d + ^2x k -H Sik p 3tk +L)Y = 

are consistent and give the equations of motion for (HH) . To prove this one checks that consistency 



follows from (6S) and in terms of L the equations (6.12) take the form 

a L + [M a , L] = (Lax); d a L - dM a + [M a , L] = L k - X p s , k (6.25) 



where a — (s, k). The Lax equation provides consistency of the first two equations in (6.24) while 



the second equation in (S.25) plays the same role for the second and third equations in (6.24); this 



equation also allows one to determine M a from L and L. 

Due to the Liouville theorem the phase flows of (HH) are restricted to the Abelian varieties 
corresponding to a level set of the Hamiltonians H s k = c s . k - This becomes simple in terms of 
action-angle coordinates defined so that the angle type coordinates are angular coordinates on the 
Abelian variety and the Hamiltonians depend only on the actions. To describe this consider 

P(A, z) = det(X - L) = \ N + hX"- 1 + ■■■ + bjX N ~i + --- + b N (6.26) 

where bj — ^Mirij where Mirij ~ principal minors of order j, and b\ = Tr(L) with b^ — det(L). 
The spectral curve C C T*E is defined via (6T) C = {P(X,z) = 0} and this is well defined since 
the bj are gauge invariant. Since L 6 H°(Y,, End(V) (g> K) the coefficients bj £ H Q (Y, g , K^) and 
one has a map (6U) p : T*7Z 9i n — ► B = 0^ v iJ°(S, K^). The space B can be considered as the 
moduli space of the family of spectral curves parametrized by the Hamiltonians H Syk ; the fibers 
of p are Lagrangian subvarieties of T*TZ 3y N and the spectral curve C is the A^-fold covering of the 
base curve S: ir : C — > E. One can say that the genus of C is the dimension g of lZ g .N (recall 
g = N 2 (g — 1) + 1). There is a line bundle C with an eigenspace of L(z) corresponding to the 
eigenvalue A as a fiber over a generic point (A, z); thus (6X) C C ker(X + L) c tt*(V). It defines a 
point of the Jacobian Jac(C), the Liouville variety of dimension g — g(C). Conversely if z S S s is 
not a branch point one can reconstruct V for a given line bundle on C as (6Y) V z = ® v ^-\i z \L v . 
Let now u)j where j = 1, • • ■ ,g, be the canonical holomorphic differentials on C such that for the 
cycles ai, • • • , a g , • • • , (3 g with oti ■ a>j — Pi ■ (3j — one has § a ujj — Sij. Then the symplectic 
form (6P) can be written in the form 

r N 

uj° = < 5L, SL >= V SXjS^j (6.27) 
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Here £j are diagonal elements of SLS where SLS^ 1 = diag(Xi, • • • , Ajv)- There results (6Z) u° = 
J c SXSt;. Since A is a holomorphic 1-form on C it can be decomposed as A = ajUj and conse- 
quently (I) u!° — J2i^ a j I Wj<5£. The action variables can be identified with (II) a,j — $ a .X. To 
define the angle variables one puts locally £ = dlog(ij}); then if (p m ) is a divisor of tp 



Pn 



u}jS£_ = ^2 w j Zo sW = (6-28) 

C m Jpo 

The 4>j are linear coordinates on Jac(C) and (III) lo° = 5a,j5(f>j. 

REMARK 6.2. For completeness let us give here some comments on the above following 
p8| , pp| , pl[ . As before take 1Z for the moduli space of stable holomorphic GL{N, C) vector bundles 
V over £ = £ ff with dim(7Z) = g = N 2 (g - 1) + 1 (note r<mfc(V) = iV and write d = deg(V). 
The cotangent space to TZ at a point V is Ty7\L = H°(End{V) <g) X) (note H°(End{V) (g> K) ~ 
H°(E, End(V)® K)) and Hitchin's theory says T*7?. is an algebraically completely integrable system 
(ACIS). This means there is a map /i : T*7vL — > £?, for a <7 dimensional vector space £?, that 
is Lagrangian with respect to the natural symplectic structure on T*1Z (i.e. the tangent space 
to a general fiber h~ 1 (a) for a £ J] is a maximal isotropic subspace relative to the symplectic 
form). Then by contraction with the symplectic form one obtains a trivialization of the tangent 
bundle T^-i^ ~ Oft-i(o) ® r*_B. This gives a family of Hamiltonian vector fields on h^ 1 (a), 
parametrized by T*B, and the flows generated by these fields on h^ 1 (a) all commute. Algebraic 
complete integrability means in addition that the fibers /i _1 (a) are Zariski open subsets of Abelian 
varieties on which the flows are linear (i.e. the vector fields are constant). In a slightly more general 
framework let K be the canonical bundle as before with total space K = T*T, g and think of a K 
valued Higgs pair (V, <fi : V — ► V <g> K) where V is a VB on S g and 4> is a K valued endomorphism. 
Imposing a stability condition this leads to moduli spaces TZk (resp. 7Z B K ) parametrizing equivalence 
classes of semistable bundles (resp. isomorphism classes of stable bundles). Let B = Bk be the 
vector space parametrizing polynomial maps p a ■ K — > J£ N where p a = x N + aix N ^ x + • • ■ + 
with ai e H°(K & ) (i.e. B = B K = ®^H°{K' si )). The assignment (V,<j>) -> det(xl - cj>) gives a 
morphism hx ■ 71 k ^ Bk (to the coefficients of det(xl — (/)))■ Then the map h is the restriction of 
hx to T*1Z which is an open subset of 7V K and dim(B) = g (miraculously - see Remark 6.1). The 
spectral curve £ = S a defined by a G Bk is the inverse image in K = T*£ of the zero section of 
Tt® N under p a : K — > K. N . It is finite over £ of degree N and for £ a nonsingular the general fiber 
of Iik is the Abelian variety Jac(£). 



REMARK 6.3. We omit discussion from 
complicated and refer also to fh' 




since the notation becomes 
1481 h78|, pM |93 for more on all of 



this ([156] is an excellent reference for symplectic matters). 



REMARK 6.4. There are two other versions of the Hitchin approach related to physics, namely 
those of Mq, f49[ and pO, f79| . Both serve as a short cut to some kind of understanding but both are 



flawed as to details (cf . [pl|, 148 for more detail but with a lack of concern regarding connection 
ideas). Let us look at |70j [79| as the most revealing, especially |79[ . Thus let £ be a compace RS 
of genus g and let G be a complex Lie group which can be assumed simple, connected, and simply 
connected. Let A be the space of g valued (0, 1) gauge fields A = Agdz on £ and take for T*A 
the symplectic manifold of pairs (A, <j>) where <f> = $ z dz = <\>dz is a g valued (1, 0) Higgs field. The 
holomorphic symplectic form on T* A is (XXIV) J„ Tr6(j)5A where Tr stands for the Killing form 
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on g suitably normalized. The local gauge transformations h £ Q = Map(E, G) act on T*A via 

A -> A h = hAh- 1 + hdh- 1 ; <f> -> <p h = h^h" 1 (6.29) 

and preserve the symplectic form. The corresponding moment map fi : T*A — > g* ~ A 2 (E)(g)(7 takes 
the form (XXV) fj,(A, <j>) = dip + A<j> + (pa (note here that Adz A (j>dz + <fidz A ^4c?z = [A, tj>]dz A 
so one has adjoint action). The symplectic reduction gives the reduced phase space (XXVI) V = 
fx~ 1 {{0})/G with the symplectic structure induced from that of T*A. As before V can be identified 
with the complex cotangent bundle T*Af to the orbit space N = A/Q where N is the moduli 
space of holomorphic G bundles on E (of course the identification here really should be restricted 
to gauge fields A leading to stable or perhaps semi-stable G bundles but is rather cavalier 
about such matters and we are delighted to follow suit). Now the Hitchin system will have V as 
its phase space and the Hamiltonians are obtained via (XXVII) h p (A,cf>) — p{4>) — p(Q z )(dz) dp 
where p is a homogeneous Ad invariant polynomial on g of degree d p . Since h p is constant on the 
orbits of Q it descends to the reduced phase space (XXVIII) h p : V — ► H°(K dp ). Again K is 
the canonical bundle of covectors proportional to dz and H°(K dp ) is the finite dimensional vector 
space of holomorphic d p differentials on E. The components of h p Poisson-commute (they Poisson- 
commute already as functions on T*A since they depend only on the "momenta" <ft). The point 
of Hitchin's construction is that by taking a complete system of polynomials p one obtains on V a 
complete system of Hamiltonians in involution. For the matrix groups the values of h p at a point 
of V can be encoded in the spectral curve C defined by (XXIX) det(<f> — £) = where £ £ K. This 
spectral curve of eigenvalues £ is a ramified cover of E; the corresponding eigenspaces of <f> form a 
holomorphic line bundle over C belonging to a subspace of Jac(C) on which the h p induce linear 
flows. For example for the quadratic polynomial pi — (l/2)Tr the map h P2 takes values in the space 
of holomorphic quadratic differentials H°(K 2 ). This is the space cotangent to the moduli space M. 
of complex curves E. Variations of the complex structure of E arc described by Beltrami differentials 
5 fx = b\x\d z dz such that z 1 = z + 5 z with d s Sz = Sfif gives new complex coordinates (see below for 
connections to the notation in (|6.6| )). The Beltrami differentials 5fi may be paired with holomorphic 
quadratic differentials j3 ~ (3dz 2 via (XXX) (/3, 8p) = J s f35/i. The differentials <5/i = d(8£) for a 
vector field on E describe variations of the complex structure due to diffeomorphisms of E and they 
pair to zero with f3. The quotient space i? 1 (_ft' _1 ) of differentials S/j, modulo d(6£) is the tangent 
space to the moduli space M. and H°(K 2 ) is its dual. The pairing (XXX) defines then for each 
[5(i] £ if 1 (if -1 ) a Hamiltonian (XXXI) hsu = h p2 5jjL and these commute for different Sfj,. 



REMARK 6.5. In (|6.6|) one has ji = fi(z,z)d z ® dz with w — z — e(z, z), w = z, and 



fi = de/(l — de) (which should be \i = —de/(l — de) — dw/dw) while from |79 the notation is 
z' = z + Sz and Sfx = Sfi^d z dz with d(Sz) = 6 fig. Thus Sz ~ — e and one surely must think of 
8fi z s ~ fx{z, z) in which case — de ~ fi (adequate for small de). Thus for small de the notations 



of (6.6)~ |175[ and |79j ca n be com pared at least. To co mpare ideas with the Koda ira-Spencer 



theory of deformations (cf. |122[ [167|) one can extract from |109| (cf. also |6fj, |9S|, |103| ). Thus the 
space of infinitesimal deformations of E is determined by i? 1 (E, G) where is the sheaf of germs of 
holomorphic vector fields on E. This in turn can be identified with the tangent space To(T(E)) of the 
Teichmiiller space T(E) (or T g ) at the base point (E, id.). We further recall that the moduli space M g 
corresponds to T(E)/Mod(E) where Mod(E) is the set of homotopy classes of orientation preserving 
diffeomorphisms E — ► E (modular group). Recall M. g is the set of biholomorphic equivalence classes 



of compact RS of genus g with dimension 3g — 3. Now, more precisely, going to [ 109 1 for notation 
etc., for a given RS E consider pairs (i?, /) with orientation preserving diffeomorphism / : E — > R. 
Set (R, f) = (S,g) if g o f^ 1 : R — > S is homotopic to a biholomorphic map h : R —> S. Then 
write [(R, /)] for the equivalence class and the set of such [(R, /)] is T(E). For [(R, /)] £ T(E) work 
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locally: (U,z) -> (V,w) with f(U) C V, and write F — wo foz^ 1 with fi f = F- Z /F z = dF/dF. This 
is the Beltrami coefficient and one can write fi = df/df in a standard notation. There are transition 
functions for coordinate changes leading to an expression fif — fi(dz/dz) where jif is a (— 1, 1) form. 
In terms of the space M(E) of Riemannian metrics on E one has T(E) ~ M(E)/.Di//o(E) and 
.M s ~ M(E) / Dif f+(E) where Diff + ~ orientation preserving diffeomorphisms of E and Diffo ~ 
elements in Dif f + homotopic to the identity. Next one writes ^(E) for the space of holomorphic 
quadratic differentials 4> = 4>{z)dz 2 and </> € _42(E)i if ||0||i = 2 J s |0| 2 <ix<iy < 1 (cf. below for the 
factor of 2). Then T(E) is homeomorphic to -42(E) (and hence to R 6 ™~ 6 ). One notes again the 
natural pairing (XXX) of A2 with Beltrami forms via (<j),jif) — J^fj,(f>dzdz formally (recall for 
z = x + iy and z = x — iy one has dx f\dy = (i/2)dz A dz so strictly 2 J g dxdy = i J g dz A dz). Now 
it is easy to show that there is an isomorphism 

and 9 = 0(K- V ) with -4 2 (E) = #°(S, lfi (K )) (recall if - holomorphic cotangent space of E). 
We note that fl^E^ ' 1 ^- 1 )) - Beltrami differentials idzj/dzj} while w S H ^^ ^ ^- 1 )) 
corresponds to a C°° vector field v ~ {«j(9/9zj} so 9w = {/j,j(dzj /dzj} where /Xj = dvj/dzj. 
Consequently there is a canonical isomorphism (XXXII) Ao (i5*) _1 : iJ 1 (E,6) — > ^(E) where 
A[/x](</>) = J^fitpdxdy for /i e 7?°(E, f 0,1 (Ar -1 )) and S ^(E). This gives the isomorphism 
_ff 1 (E,0) ~ To(T(E)) mentioned above. These facts will help clarify some remarks already made 
above (note H 1 ^- 1 ) ~ F^E.O) and its dual H°(K 2 ) - _4 2 )- 



We return now to |l75| again and consider HID in the (scaling) limit k —> (critical value). 
One can prove that on the critical level HID coincides with the part of HH relating to the quadratic 



Hamiltonians in ( 3.22 ). Note first that in this limit the A connection is transformed into the Higgs 
field (A -> 4> as k -> - cf. ( |6T0| ) and recall L = /"V/) and therefore FB{Z,G) -> r*^ g ,Ar 



(perhaps modulo questions of stability) . But the form w on the extended phase space V appears to 
be singular (cf. (6G) and ( |6. 13 ) ) and to get around this one rescales the times (IV) t = T + nt H 



where t H are the fast (Hitchin) times and T the slow times. Assume that only the fast times are 
dynamical, which means (V) 5fi(t) = K^ s ^"<5tf r where one writes = dn s . After this rescaling 
the forms (6G) and ( |6. 13 ) become regular. The rescaling procedure means that we blow up a vicinity 
of the fixed point /i° in A4 g , n and the whole dynamics is developed in this vicinity. The fixed point is 
defined by the comple x coo rdinates (VI) wq = z — ^ s T s t s {z, z ) wit h wq — z. Now compare the BA 
function ip of HID in ( 6.11 ) with the BA function Y of HH in ( |6.24 ). Using a WKB approximation 



we assume (VII) t/j — Qexp[(S Q / k) + S 1 ] where $ is a group valued function and S°, S 1 are diagonal 
matrices. Put this in the linear system ( |6.11 ) and if (VIII) dS°/dwo = = dS°/dt^ then there 



are no terms of order n 1 . It follows from the definition of the fixed point in the moduli of complex 
structures (VI) that a slow time dependent Sq emerges in the form 



S a = S° ^T li '-',T l \z-^2T a e.(z,z)^ 



.31) 



In the quasiclassical limit put (IX) dS° = A. so in the zero order approximation we arrive at 
the linear system of HH (namely ( |6.1l| )) defined by the Hamiltonians H s ^ (k = 1,2) and the BA 
function Y takes the form (X) Y — <&exp{J2 s tf {dS° /dT s )\. The goal now is the inverse problem, 
i.e. to construct the dependence on the slow times T starting from solutions of HH. Since T is 
a vector in the tangent space to the moduli space of curves M g it defines a deformation of the 
spectral curve in the space B (cf. (6U)). Solutions Y of the linear system ( |6.24 ) take the form 
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Y = Qexp[^2 s ty Cl s ] where the fl s are diagonal matrices. Their entries are primitive functions of 
meromorphic differentials with singularities matching the corresponding poles of L. Then in accord 
with (X) we can assume that (XI) ddS/dT s — d£l s so the T s correspond to Whitham times (along 
with the t s from ( 6.18) ), ( ]6.20| ), and ( |6.21 ) - note t s = T s + ntf so t s and T s are comparable in this 
spirit). These equations define the approximation to the phase of tp in the linear problems fl6.1l| ) of 
HID along with (XII) ddS/daj = dujj. The differential dS plays the role of the SW differential and 
an important point here is that only a portion of the spectral moduli, namely those connected with 
H s .k for k — 1,2, are deformed. As a result there is no matching between t he ac tion parameters of 
the spectral curve a,j (j ' = 1, • • • ,g) in (II) and deformed Hamiltonians (cf. |l94| and Section 3). 

Next the KZB equations (Knizhnik, Zamolodchikov, and Bernard) are the system of differential 
equations having the form of non-stationary Schrodinger equations with the times coming from M 9i n 
(cf. ]114| ). They arise in the geometric quantization of the moduli of flat bundles F -B(S, G). Thus let 

V = V\®- ■ -®V n be associated with the marked points and theHilbert space of the quantum system is 
a space of sections of the bundle £ v ^ ant (S ffi „) with fibers FB(E, G) depending on a number K i uant - 
it is the space of conformal blocks of the WZW theory on S g n . The Hitchin systems are the classical 
limit of the KZB equations on the critical level where classical limit means that one replaces operators 
by their symbols and generators of finite dimensional representations in the vertex operators acting 
in the spaces Vj by the corresponding elements of coadjoint orbits. To pass to the classical limit in 
the KZB equations (XIII) (K quant d s + H S )F = (note these are kind of flat connection equations) 
one replaces the conformal block by its quasi-classical expression (XIV) F = exp(!F '/%) where 
h = (where k 



..quant 



/Ti as in [ 134 1 ) and consider the classical limit K 



quant 



which 

leads to HID as in (6.21) so (XIII) is a quantum counterpart of the Whitham equations). One 
assumes that the limiting values involving Casimirs C\ (i = 1, • • • , rank(G)) and a — 1, ■■■,n), 
corresponding to the irreducible representations defining the vertex operators, remain finite and 
this allows one to fix the coadjoint orbits at the marked points. In this classical limit (XIII) 
then is transformed into the HJ equation for the action T — log(j) of HID, namely ( 3.21 ) (note 
K quant d s F — > K,d s T etc.). The integral rep rese ntations of conformal blocks are known for WZW 
theories over rational and elliptic curves (cf. 112 , 70| , (73)) so (XIV) then determines the prepotential 
T of HID. The KZB operators (XIII 
■Mg,n w hh fibers £ v tK mant (£ 9i „) (cf. [72 



) lay t he role of flat connections in the bundle r pi uant over 
"In)) with 



[ K quant d s + H s , K quant d r + H r ] = 



(6.32) 



In fact these equations are the quantum counterpart of the Whitham hierarchy ( 5. IS ) 



7 WHITHAM AND SEIBERG-WITTEN 



We mention first that a general abstract theory of Whitham equations has been developed in 
HI 105 , 106 , 127 , 12£, 12£] and this was summarized in part and reviewed in 

(a 



|26|, |33|. One 

, N) and a collection 



deals with algebraic curves S s ,jv having punctures at points P a (a = 1, 
of Whitham times Ta and corresponding differen tials dilA is constructed. The Whitham equa- 
tions (XXXVIII) dAdVlB = dsdVlA persist as in ( 2.38| ). The construction of differentials becomes 
complicated however and we will not deal with this here. Rather we will follow the approach of 

© H M " 



112 




(summarized as in 



accurate development follows |83J, 18E] (as recorded in |26L E8I]) and we will extract here from 




The most revealing and 
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7.1 Background from |83 



We take a SW situation following |§ H EH EH EH ^ H II H HI EOH HI HI HH^ HH I 



and recall the SW curves S s (<? = N — 1) for a pure SU(N) susy YM theory 



tfetjvxiv[£H - A] = 0; P(A) = A JV [w + 



w H 



(7.1) 



AT 



P(A) = A Ar -^ Ufc A JV - fc =n(A-A J ); u fc = (-l) fc £ V"Ai, 
2 l ii<---<4 

Here the Uk are Schur polynomials of hk — y i a the formula (A) log(X~ N P(X)) = 

— ^ fc (/ifc/A' £ ); there are g = N — 1 moduli Uk and we refer to [ 1 13 1 for the Lax operator L. Thus 
i*o = 1, u\ =0, U2 — hi, U3 = hs, U4 = h,4 — (\/2)h\, etc. (hi — for SU(N)). One also has the 
representation 

1 



y 



P 2 (A)-4A^ V ; y = A 



,iV 



w 



(7.2) 



giving a two fold covering of the punctured Riemann sphere with paramet er A . Such Toda chain 
curves are characterized by a function (B) 2A N w = (P + y). Note also from (71) - (7.2) one obtains 



5P + P'SX = NP6log(A) + y—; SP = -V \ N - k 6u k ; P' 

in — ' 



<9P 
9A 



On a given curve (fixed it/c and A) 

dw P'dX 



dw XdP 

; dbsw = A — = 

y w y 



ddS, 



sw 



duk 



X N ~ k dw X N ~ k d\ 



P> 



dv k : k 



.N 



where the dv k are g = N — 1 holomorphic one forms with 



a, 



dSsw] & 



ik 



dv k 



dog 
du k 



(7.3) 

(7.4) 
(7.5) 

(7.6) 



* Ai J Ai 

and du>i — ((Ti k )~ 1 dv k are the canonical holomorphic differentials with (C) § A dujj — 6ij; § B . dujj = 
P>ij. Note in (7.5) it is necessary to assume w is constant when the moduli u k are varied in order to 
have ddSsw/duk = dv k holomorphic. Now the periods = § A . dSsw define the as functions of 
Uk (i.e. hk) and A, or inversely, the Uk as functions of a, and A. One proves for example (see (7.3£) 
below) 

du k du k , , 

and generally A and T\ can be identified after suitable scaling (cf. 17, [H |H |64|, H). This is 
in keeping with the idea in (64j| that Whitham times are used to restore the homogeneity of the 
prepotential when it is disturbed by renormalization (cf. also 



For the prepotential one goes to 113, |168[ for example and defines differentials (D) dfl n 
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(£ " 1 + 0(l))d£ for n > 1 with § A <Kl n — (pick one puncture momentarily). This leads to the 
generating function of differentials 



W(t C) = E nC^dCdn^) = d(d c E(£, C) 



(7.8 



where J5 is the prime form (cf. |7]J| ) and one has 



(HF 



0(D = E^^rC"- 1 rfC + 0(i) 



(7.9) 



(cf. here |2(| [27j] for comparison to the kernel K (fi, A) = l/(P{p) - P{\)) of Section 2.4). One can 
also impose a condition of the form (E) ddQ / d (moduli) = holomorphic on differentials dfl n and 
use a generating functional (note we distinguish dS and dSsw) 



a 

dS = 2_. T n dCl n — oii<ko 



oo 

E 



(7.10) 



(we have added a Todflg term here even though To = in the pure SU(N) theory - cf. |6J, P-72| ). 
Note that there is a possible confusion in notation with dfl n since we will choose the dCl n below 



to have poles at oo± and the poles must balance in (7.1C). The matter is clarified by noting that 
(D) determines singularities at £ = and for our curve £ = ~ oo± via £ = ur^ 1 ^ which in 
turn corresponds to P(X) 1 ^ N for A = 1 (see ( 7.19] ), (7.28), and remarks before (7.16) - cf. also 
( [7.30 ) for dfl^). Thus in a certain sense d£l n here must correspond to dfl^ + dfl~ in a hyperelliptic 
parametrization (cf. ( 7.30| )) and T n ~ T+ = T~ (after adjustment for the singular coefficient at 
oo±). The periods 014 



dS can be considered as coordinates on the moduli space (note these 
are not the a, of |23|, 168 ). They are not just the same as the but are defined as functions of hk 



Ai 



and T n (or alternatively hk can be defined as functions of a, and T n so that derivatives dhk/dT n 
for example are nontrivial. One will consider the variables oij and T n — —(l/n)Res£=o£ n dS(£) as 
independent so that 

8dS , ddS JO 

oai dT n 

Next one can introduce the prepotential F((Xi, T n ) via an analogue of af = d!F/dai, namely 



OF 
da* 



Then one notes the formulas 
d 2 F 1 



dT m dT n 2irin 



Resot; 



dF 



ddS 



1 



2nin 



Res C n dS 



dT„ 



2irin 



ResoS, n dQ m = 



1 



27TITO 



Res £- m dn n 



(7.12) 



(7.13) 



(the choice of £ is restricted to w ±x / N in the situation of (7.2)) and factors like n _1 arise since 
£~ n ~ l d£ = — d(£ _ ™/n). Below we use also a sl ightl y differ ent n ormalization dQ n ~ ±w ±n / N (dw/w) = 
(N/n)dw ±n / N near oo± so that residues in ( 7.12 ) and ( 7.13 ) will be multiplied by N/n instead of 
1/n. Note also that in accord with the remarks above Res^ = o will correspond to the sum of residues 
at 00 ± involving £ = m^ 1 /" which in turn corresponds to P(X) 1/N . By definition F is a homogeneous 
function of at and T n of degree two, so that 



2F = a, 



dF 

dai 



dF _ ^ d 2 F 
n &T n 1 3 daidaj 



2a.,T„ 



d 2 F 

da t dT n 



d 2 F 

dT n dT„ 



(7.14) 
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Note however that F is not just a quadratic function of on and T n ; a nontrivial dependence on these 
variables arises through the dependence of du>i and d£l n on the moduli (such as Uk or hk) which 
in turn depend on on and T n . The dependence is described by a version of Whitham equations, 



obtained for example by substituting (7.1C) into (7.11). Thus 



ddS 



(Kin + T r 



ddQ m dui 
du t dT n 



G?fL 



^ m dT n 



ddtl„ 
due 



dfl r 



(7.15) 



(since § A dfl n = 0). Now (cf. ]23|, |113fl ) to achieve (E) one can specify dd£l m /due = Y^Pijdwj so 
§ A (dd£l m /due) = § A . Y^PTjd^j = 0$ and since dtl n = dSl n + ^Cjdujj of necessity there results 
(F) J2(due/dT n ) Tmftel = ^{due/ 'dT n )au = — c™ which furnishes the Whitham dynamics for u p 
in the form du p /dT n — — c"a lp (cf. [^3], |113| - the formulas in Q and copied in an earlier version 
of |2(| were too rushed). 



Now the SW spectral curves (7.2) are related to Toda hierarchies with two punctures. We recall 
(jfjj) - Q and note therefrom that (G) to* 1 = (1/2A N )(P ± y) ~ (1/A JV )F(A)(1 + 0(\- 2N )) 
near A = oo± since from y 2 = P 2 - 4A 2JV we have (y/P) = [1 - (4A 2Ar /P 2 )] 1 / 2 = (1 + 0(\- 2N )) 
so w ±J = (P/2A W )[1 ± (y/P)} = {P/2A N )[2 + 0(\- 2N )}. One writes w(X = oo+) = oo and 
w(X = oo_) = with £ - ly^i/JV (i.e. £ = ^-i/JV _ A -i at QO+ and £ = _ A -i at QO _). 

Near A = ±oo by (G) one can write then w^ 1 ^ ~ P(A) 1 / A ' (for A = 1) in calculations involving 
W ±n/N w ^ n n <- 27V. The to parametrization is of course not hyperelliptic and we note that (D) 
applies for £ = w ±1 / JV with all d£l n and later only for £ = A -1 in certain differentials dfl n . It would 
now be possible to envision differentials dVL^ and dQ„ but dtl n = dti+ + dtl~ is then clearly the 
only admi ssib le object (i.e. dCl n must have poles at both punctures); this is suggested by the form 
dw/w in (7.4) and the coeffic ients of w n / N at oo + and of w~ n / N at oo_ must be equal (see also 
remarks above about d£l n in ( 7.10| )). This corresponds to the Toda chain situ atio n with the same 
dependence on plus and minus times. Moreover one takes differentials dCl n for (7.2) (A = 1 here for 
awhile to simplify formulas) 



dfl r 



R n (X) 



dw 



P 



.n/N 



(A) 



dw 



(7.16) 



These differentials satisfy (E) provided the moduli derivatives are taken at constant w (not A) and 
we can use the formalism developed above for £ = w Tl / N . Note the poles of dtl n balance those of 
dfl n = dfl n — (Xa(Ju dfl n )du>j and we have 



(7.17) 



where we keep n < 2N for technical reasons (cf. [ 1 89 1 ) . The SW differential dSsw is then simply 
dSsw = dfli, i.e. 



dS\ Tn —g n 1 = dSsw', a i\T n =6 n i = a i> a 



i lT„=5„ i 



(7.18) 



With this preparation we can now write for n < 2N, where ~ P(A) 1 / JV , with dti n as indicated 

in Q , 

f)f? AT / \ 

dT = 2^i ( Res °°+ wn/NdS + ReSoo_w- n ' N dS) = (7.19) 



N 
2mn 



[Re Soo+ w n ' N +Re Soo _w- n / N ) ^£ T m P? /N (X^J ^ = 
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A " ^T m i? eSco (p™ /Ar (A)dP"/ JV (A)) =-J^Y, T ™ Res °o(P n/N WdPf N (\)) 

ra rn 

One can introduce Hamiltonians here of great importance in the general theory (cf. ||(| Q |3[ |113| , 
137 1) but we only indicate a few relations since our present concerns lie elsewhere. Then evaluating 

ResooP n ^ N (\)d\ = —H n+1 (7.20) 



at T n = 5 n> i (7.19) becomes 

dF N 2 „ .. N 



dT n iim 2 iirn 
where 

H n+1 = -^Res^dX = £ tlfl ( » ) ^ ^ K . . . hi 

k>l ii H Hk=n+1 



hn+1 ~ih ^ h^ + Oih 3 ) (7.21) 



This can be rephrased as 



H = ^ E - + • • (7.22) 

where 



H m+ i,n+i = -—Resoo (P n/N dP™ /N ) = -H n+1 . m+1 : (7.23) 



W n+1 = W n+ i ;2 = -—Re Soo P n ' N d\ = h n+ i + 0(h 2 ) 
n 



Now for the mixed derivatives one writes 

d 2 F f _ I 



dQ n = — — Res C n dLJ t = (7.24) 



dctidTn J B . 2nin 
A fijesoo.w"/^^ + J Res 00 _w-" /Ar dw l N ) = —Re Soo P n/N du 



. I UUJ. (i I l^J LI 

Imn V / tirn 

Next set (H) P"/JV = E^l&A* so that (I) Re Soo P n / N dtu l = E-ooPnk Res ooX k d^i. Then e.g. 

d.M) - ^dv\x) = = (1 + ° {X ~ 2N)) = 

= -^ d J^m dX(l + 0{ X-^)) (7.25) 
From (A) and aj£ — du^/daj one obtains then 

^(A> (1 + (A-,) - = E (7 - 26 » 

n>2 ' n>l 

so for k < 27V, (J) Res oa X k duJi = dhk+i/dai. Further analysis yields (K) ReSooW n / N diOi 
= ResooP n t N duji = dHn+i/dcn leading to 

a — = -. — i?es oc P(A) i du l = (7.27) 

oa.iOl n iirn inn oa.i 
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For the second T derivatives one uses the general formula (7.13) written as 



d 2 F 



1 



-Res C n dttr 



N 



dT n dT m 2itin ' N "' 2-nin 
while for the second a derivatives one has evidently 

o 2 f _ r 

da.ida 3 J B . 



Resoo + w n + ReSoo_w n/N d£l m ^ (7.28) 



dujj = Bij 



(7.29) 



Note that one can also use differentials dfl defined by (D) with £ = A 1 (not £ = w Tl ^ N ); recall 
oo ± ~ (±, A — > oo ) in the hyperelliptic parametrization. This leads to 



dttl 



± 



(w ±n ' N + 0(1)) — = -dw ±n ' N + 
V / w n 



n 



-^kp^dX 
i 

Putting (H) and ( |7.30|) into (|7.28|) gives then 

d 2 F 



n 1 



dT m dT n 



N 2 m 

Trunin < * 



(7.30) 



(7.31) 



where dflf = dfl~l 



Further analysis in 0] involves theta functions and the Szego kernel (cf. [fnl |83|). Thus let E be 
the even theta characteristic associated with the distinguished separation of ramification points into 
two equal sets P(A) ± 2K N = J]f (A - r± ). This allows one to write the square of the corresponding 
Szego kernel as 



*!(a,aO 



P(X)P(fi) - 4A 2JV + y(X)y(jj) dXdfx 



2y(A)y(^) 



(A-m) s 



We can write (cf. @) 



n>l " 



2.y 



(l + 0(A- 2JV ))dA near oo ± 
0{X- 2N dX near oo T 



and utilize the formula 



(7.32) 



(7.33) 



dzidzj 



log6 E {d\B) 



(cf. |71|, |8J). Here one uses ([7.261) and (|7j) to get (1 < n < 27V, ( ~ 

= E 5^ Ht^) ! = A- 1 *^) - p^(A); 



(7.34) 



dO n (A) = A"- 1 (*|+(A) + *|-(A)) - 2 P ^ i (A) 



(7.35) 
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where (L) p % n = (l/n)(dh n+ i/daj)dfjlog 9e(0\B). From this one can deduce with some calculation 



d 2 T 



dT m dT r ' 



N 
nin 



H 



Next one notes 



Then there results 



E 



duk 



dlog(A) 



2N dHm+l dHn+l Q 2 

mn da,i da lJ 
ddSsw 



8?Aogd E (P\B] 



duk 



ddS s 



ir 



Ai dlog(A) 



E 



du k daj 
dlog(A) duk 



ddS s 



M dlog(A) 



-N 



P dw 
~P 1 ~ 



-N 



PdX 

Ai y 



(7.36) 



(7.37) 



(7.38) 



P + V -dX = -2NA N I ^ = -2NA N 4 wdv N 



= -N 

Ui y Ja z y ja % 

Here one is taking dSsw = Xdw/w and using ( |7.3|) in the form (§P = Sw = 0) SdSsw 1 'Slog(A) = 
5X(dw/w)/5log)A) = (NP/P')(dw/w). Then (|7.4|) gives NPdX/y and the next step involves 
§ A d\ = 0. Next (B) is used along with @. Note also from (M) XdP = X[NX N ^ - J2( N ~ 



k)u k X N - k - 1 ]dX = NPdX + £ kukX dX one obtains via (7.4), (UM, and (7.38) (middle term) 



Eduk da, 
— , ; — — ~ = N 



dlog(A) duk 



PdX 



Ai y 



(7.39) 



[XdP 



^ ku k X 



N-k 



dX 



/a 4 V y " y 

which evidently implies (7.7). Now from ( |7.38 ) there results 



duk dai 
dlog(A) duk 



= 2NA 



N 



wdX 



= N 



a, y 



E ku k 



P + y 

At y 



daj 
du k 



dX = 



(7.40) 



dh 



2N i> *|(A) = 2Np\ = 2N^d?Aog0 E {d\B 



from which 



At 

du k 



dlog(A) 



-2N 



da 

du k du 2 „ 2 , o m 
da-da- d - l ° 99Eim 



(7.41) 



Here one can re place Uk b y any function of Uk alone such as hk or r H n +\ (note u 2 = h 2 ). Note also 
(cf. J2(| Q 106, 155 1) that identifying A and T\ (after appropriate rescaling hk — * I^/ifc an( l 
Uk -> TfWfe) one has (/3 = 2JV) 



3kj(A) ~ 2 7 ri^ i '" z/ (7 - 42) 

(this equation for dF/dlog(A) also follows directly from (|7.20| ) - fl7.22|) when T„ = for n > 2 since 
H 2 = h 2 ). 

Finally consider O) in the form (N) P'SX - J2k XN ~ kSu k = NP5log(A). There results for 



(5a,; = (cf. ( fQsl )) 

(5a, = 



. , v-^ c f X N k dw 

A 5X -^ 5Uk t A —- 



N5log{A) 



P dw 



(7.43) 
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On the other hand for on = Tkh + 0{T 2 , T 3) • • •) 



F d 

5a % =a l 5log{T 1 )+T 1 d> 5X—+0(T 2 ,T 3 , 
JAi w 



so for constant A with T„ — for n > 2 (while m and T n are independent) Sai — implies 



E 



du k 



dlog(Ti 



Since XdP = NPdX + J2 k ku k X N - k d\ it follows that (cf. (7.7)) 



dm. 



9io5(Ti) 



a MA) 



kui 



XdP 



duk 
'' da,. 



(7.44) 



(7.45) 



(7.46) 



(cf. ( |7.7| ) - note the evaluation points are different and ai = T\a,i + 0(T2,T 3 , •••)). This relation is 
true for any homogeneous algebraic combination of the u k (e.g. for h k and TL k ). We will return to 
the log(A) derivatives later. For further relations involving Whitham theory and A derivatives see 

0, gj n, n, n ~ 



105, 106] , 113] and references there. 



THEOREM 7.1. Given the RSjO)- ( |7~2j) on e can determine dF/dT n fr om (^20| ), d 2 F/daidT n 
from ( ffT27| ), d 2 F/dT n dT rn from (|7.28|) orj7j|), and d 2 F/daidcu from (|7.29|) (also for T„ = 
<5n,i, Qii = a.; as in (7.18) with A = 1). Finally the equation (7.42) for dF / 'dlog(K) corresponds 
to an identification T\ ~ A and follows from ( 7.20| ) - ( 7.22| ) when T n = for n > 2. Therefore, 
since F ~ -Fsw for Ti = 1 we see that all derivatives of Fsw are determined by the RS alone 
so up to a normalization the prepotential is completely determined by the RS. We emphasize that 
Fsw involves basically T n = 8 n _\ only, with no higher T n , and ai — ai, whereas F involving ai 
and T n is defined for all T n (cf. here pi|, |8j|). In fact it is really essential to distinguish between 
Fsw = F sw and general F — F w = Fw — Fwhit and this distinction is developed further in p8| . 
The identification of A and T\ is rather cavalier in and it would be better not to set A = 1 
in various calculations (the rescaling idea is decpetive although correct). Thus since F sw arises 
from F w by setting T n = 6i, n it is impossible to compare AOaF with T\d\F w directly. Indeed 
a statement like equation (7.42) is impossible as such since T\ doesn't appear in F sw . In m 



took a simple example (elliptic curve) and computed everything explicitly; the correct statement 
was then shown to be T\d\F w = A8aF w . In addition it turns out in this simple example that the 
Whitham dynamics lead directly to homogeneity equations for various moduli and the prepotential. 



7.2 Connections to ]|168| 

The formulation in Section 7.1, based on |83|, d iffer s from [|64| , 113| , 168 1 in certain respects and we 
want to clarify the connections here (for [105, 106, 127, 128, 129 we refer to the origi nal papers 
and to 1 26, |3j|). Thus first we sketch very briefly some of the development in [ 168 1 (cf. also 
[fl6| , |23|). Toda wave functions with a discrete parameter n lead via T k = et k , T k — et k , To = —en, 
and dj = ieOj {= § A dS) to a quasiclassical (or averaged) situation where (note T does not mean 
complex conjugate) 



dS 



9 



aidiOi 



n>0 



n>l 



(7.47) 
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1 (« dF dF dF \ 

F = 2 £ a ^ + £^ + E T «Md 



(7.48) 



where (ffi n ~ dO+ cMl„ ~ <iil n , T n ~ T_„, and near P + 



while near P_ 



dfl+ = 



dSl- = 



-5 n0 2 1 - E r ™ 
i 

oo 



dz (n > 1); 
cb (n > 0) 

dz (n > 0); 
dz (n > 1) 



(7.49) 



(7.50) 



Here dQo has simple poles at P± with residues ±1 and is holomorphic elswhcre; further dQ^ = 
dil,Q = dflo is stipulated. In addition the Abelian differentials dSl^ for n > are normalized to have 
zero Aj periods and for the holomorphic differentials diOj we write at P± respectively 



dujj = — <jj m z m 1 dz; dujj = — a 



Jin" 



v dz 



m > 1 m > 1 

where z is a local coordinate at P±. Further for the SW situation where (g = N — 1) 
d5 = AdP = APU\. y2 = p2 _ A 2iv : p(A) = x n + y UN _ kX " 

v v V 



(cf. (7.1) where the notation is slightly different) one can write near P± respectively 



dS 



OF 

dfn 



(7.51) 



(7.52) 



(7.53) 



leading to 



dS = ( - E nTnZ-"- 1 - Toz- 1 - ^f^ 1 ] dz 

n>l n>l 



F=-\Y.^ i fjS-Y J T n Res +Z - n dS- 



(7.54) 



T n Res-z n dS — To[Res+log(z)dS — Res-log(z)dS] 



n>l 
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(the 2iri is awkward but let's keep it - note one defines a® = § B dS). In the notation of [ 168 1 one 



can write now (•) h = y + P, h = —y + P, and hh = A with h 1 ~ z at P+ and z ~ h at 
P- (evidently /i~wof Section 2). Note also (••) h + (A 2N /h) = IP and 2y = h- (A 2N /h) yielding 
y 2 = p 2 _ A 2N and ca i cu i at i ons j n [jiTjgj g i ve (. m m ) dS = XdP/y = Xdy/P = Xdh/h (so h ^ w 
in Section 7.1). Further the holomorphic du>i can be written as linear combinations of holomorphic 
differentials (g = N — 1) 



dv k = 



X^dX 

y 



2g+2 



(k = l,---,g); y 2 = P 2 -1= Y[(X-X a ) 



(7.55) 



(cf. (7.5) where the notation differs slightly). Note also from fl7.55| ) that 2ydy = Y^i N Yia^pi^ ~ 
X a )dX so dX = corresponds to y = 0. From the theory of |168| (cf. also ^3|) one has then Whitham 
equations 

ddu>j ddtUi ddaj.j. ddflA ddfln ddflA 

(7.56) 



ddiOi ddiOj ddflA ddfli 



den da,j ' 8Ta 
along with structural equations 



dai 



ddQ, 



dT A dT B 



ddS , ddS m , ddS , n _ ddS 

aw,; ttt^ — = dil„; = aS2„ ; -rr^ = alio 



dai 



dT n 



dT 



(7.57) 



Finally we note that in 1 168 one presents a case (cf. also |26J ) for identifying N = 2 susy Yang-Mills 
(SYM) with a coupled system of two topological string models based on the An-i string. This 
seems to be related to the idea of ti fusion (cf. J59|). 



8 SOFT SUSY BREAKING AND WHITHAM 

The idea here is to describe briefly some work of Edelstein, Gomez-Reino, Marino, and Mas about 
the promotion of Whitham times to spurion superfields and subsequent soft susy breaking J\f = 2 — > 

8.1 Remarks on susy 

We begin by sketching some ideas from (lj, [u| where a nice discussion of susy and gauge field 
theory can be found. In particular these books are a good source where all of the relevant notation 
is exhibited in a coherent manner. One says that a quantum field theory (QFT) is renormalizablc 
if it is rendered finite by the renormalization of only the parameters and fields appearing in the 
bare Lagrangian (for renormalization we refer also e.g. to |lj, [l3| and the bibliography of pflj). We 
denote by A the renormalization scale parameter. One notes that renormalization of bare parameters 
occurs as a quantum effect of interaction and the shifts thus generated are infinite, which means 
that the bare parameters were also infinite, in order to produce a finite measured value. Further 
in order to implement gauge invariance in weak interactions for example one must find a method 
of generating gauge vector boson masses without destroying renormalizability. Any such mass term 
breaks the gauge symmetry and the only known way of doing this in a renormalizable manner is 
called spontaneous symmetry breaking. This arises when e.g. when there are nonzero ground states 
or vacua which are not invariant under the same symmetries as the Lagrangian or Hamiltonian. Once 
such a vacuum is chosen (perhaps spontaneously by the system "settling down") the symmetry is 
broken. One can then define new fields centered around the vacuum which have zero vev (vacuum 
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expectation value) and the Lagrangian expressed in the new fields will no longer have the same 
symmetry as before. Such new fields (with nonzero vev) have to be scalar (not vector or spinor) 
and are called Higgs fields; this kind of spontaneous symmetry breaking is nonperturbative (the 
vevs are zero in all orders of perturbation theory). In the case of continuous global symmetry 
in the Lagrangian there can be a subspace of degenerate ground states and massless modes called 
Goldstone bosons arise. In any event the idea now is to break a local gauge invariance spontaneously 
in the hope that the break will induce gauge boson masses while the (now hidden) symmetry will 
protect renormalizability. This is referred to as the Higgs mechanism and what happens is e.g. 
that the Goldstone bosons are "eaten" by the gauge transformed massive boson field and a scalar 
Higgs field remains (we recall that massive terms are quadratic in the Lagrangian). In the case of 
nonabelian gauge theories one includes Yukawa couplings of fermions to the scalar fields in order to 
have fermion masses emerge under spontaneous symmetry breaking. Further magnetic monopoles 
may arise in spontaneously broken nonabelian gauge theories and instantons arise in general (which 
are classical gauge configurations not necessarily related to spontaneously broken symmetry); we 
omit any discussion of these here. 

Now, turning to susy (following (lSj) one introduces a spinor geometry to supplement the bosonic 
generators of the Poincare group. This leads to a natural description of fermions and in a susy 
theory the vanishing of the vacuum energy is a necessary and sufficient condition for the existence 
of a unique vacuum. Further, every representation has an equal number of equal mass bosonic and 
fermionic states. Generally the nonzero masses of observed particles are generated by susy breaking 
effects so one looks first at representations (and their TCP conjugate representations) of the N = 1 
susy algebra that can be realized by massless "one" particle states. This leads to supermultiplets 
for N = 1 involving (A ~ helicity) (A) chiral: quarks, leptons, Higgsinos (A = 1/2) with squarks, 
sleptons, Higgs particles for A = (scalar particles) (B) vector: gauge bosons (A = 1) with gauginos 
(A = 1/2) (C) gravity: graviton (A = 2) with gravitino (A = 3/2) along with TCP conjugate 
representations (A)' A = -1/2 and A = (B') A = -1 and A = -1/2 (C) A = -2 and A = -3/2. 
For N = 2 susy one has supermultiplets involving 4-D real representations of U(2) (in the absence 
of central charges) (D) vector: A = 1, double A = 1/2, and A = (E) hypermultiplet: A = 1/2, 
double A = 0, and A = -1/2 (TCP self-conjugate) (F) gravity: A = 2, double A = 3/2, and A = 1 
along with TCP conjugations for (D) and (F). 

One goes then to superfields S(x,0,0) with Grassman variables 0, (for which a nice discussion 
is given in (l|). There will be expansions 

$(x M , 0,G)=<t>+ V29tp + 00F + id^Oa^e - -^=008^^0 - -8^^0000; (8.1) 

V 2 4 

$t = 0t + V20~i) + 00F^ - id^Oa^O + -^=000^8^ - -8^^0000 

v2 4 

for chiral superfields where ip ~ left handed Weyl spinor, (/>, F ~ complex scalar fields, and <7° = I2 
with a 1 ~ Pauli matrices (i — 1,2,3). It is useful to note also that 5F is a total divergence under 
susy transformations. Similarly vector superfields can be written 

V(x, 0, 0) = C{x) + i0 X (x) - i§x(x) + l — [M(x) + iN(x)} - % — [M{x) - iN(x)] + (8.2) 



+9a fl 0V fl (x) +U 



D - Idnd^C 
2 M 
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and 5D will be a total divergence (along with 8(d^d^C)). Here %, A are Weyl spinor fields, is a 
real vector field, and C, M, N, D are real scalar fields. Generally one refers to the coefficient of 88 
in product expansions j or $i$j&k for example as an F term and the coefficient of 8888 as a D 
term. Then susy Lagrangians involving chiral superfields will have the form 



C = Yfei*i\D + (XW($)]F+HC)~ J d 4 6j2®t®i+ (J d 2 8W(<S>) + HC 



(8.3) 



(HC ~ Hermitian conjugate) where W is called a superpotential and involves powers of <&i only 
up to order three for renormalizability. The latter equations arises since the superspace integration 
projects out D and F terms (recall / d88 = 1, J dB = 0, / d 2 6>6><9 = 1, (d/d8)f(8) = J d6f(6), etc.). 
Note that, apart from a possible tadpole term linear in the <&i, one will have 

W($) = -rriij^j + ^A,,,'l>,<l>,<h (8.4) 

There is then a theorem which states that the superpotential for TV = 1 susy is not renormaliz- 
able, except by finite amounts, in any order of perturbation theory, other than by wave function 
renormalization. Regarding susy breaking one must evidently have this since we do not see scalar 
particles accompanied by their susy associated fermions. To recognize when susy is spontaneously 
broken we need a vacuum |0 > which is not invariant under susy (or alternatively |0 > should not be 
annihilated by all the susy generators). As a consequence whenever a susy vacuum exists as a local 
minimum of the effective potential it is also a global minimum. For the global minimum of the effec- 
tive potential (physical vacuum) to be non susy it is therefore necessary for the effective potential to 
possess no susy minimum. In theories of chiral superfields one needs < 0|Fj|0 >^ for spontaneous 
susy breaking where the tree le vel e ffective potential is V — F} F, = \Fi\ 2 and F} = — dW (4>) / dfa 
(here Fi is an F term arising in ( |3.3| )). Once spontaneous susy breaking occurs a massless Goldstone 
fermion appears which for Fi will be the spinor ipi in the supcrmultiplet to which Fi belongs. When 
global susy becomes local susy in supergravity (sugra) theories the Goldstone fermion is eaten by 
the gravitino to give the gravitino a mass. For theories involving vector superfields there is also 
another possibility when there is a D term with < 0|D(x)|0 >= <fr ^ 0. Finally one notes that the 
renormalized coupling constants necessarily depend on the mass scale A but the physics described 
by the bare Lagrangian is independent of A so the coupling constants must "run" with A. The RG 
equations specify how these coupling constants vary. 

Now to couple with sugra one recalls first the Nocther procedure for deriving an action with a local 
symmetry from an action with a global symmetry. For example given So = i J d x^^d^ip invariant 
under the global symmetry ip — ► exp(—ie)ip one lets e = e(x) and considers (Jit) ip — > exp(—ie(x))ip. 
Then SSq = J d 4 x-tpj li ipd ll e — J d 4 xj fi d ll e where — •tpjVip is the Noether current. To restore 
invariance a gauge field is introduced transforming under (Jit) as (4) — > + d^e and a 
coupling term is added to Sq to obtain S = So — J d A xj^A^ = J d i xi4>^ il (d ll + iA^)ip. This S is 
invariant under (Jit) and (Jb). One can use this technique to construct a locally susy action from 
the global susy action for the sugra multiplet. Next one extends the pure sugra Lagrangian of the 
graviton and gravitino to include couplings with matter fields. Recall that whereas in global susy 
theories susy breaking manifests itself in the appearance of a massless Goldstone fermion, in locally 
susy theories the corresponding effect is the appearance of a mass for the gravitino which is the 
gauge particle of local susy. The most general global susy Lagrangian for chiral superfields is 



-Glob 



J d i 8K(<p\<i>) + J d 2 8(W(<S>)+HC) (8.5) 
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where K is a general function (since nonrenormalizable kinetic terms cannot be excluded in the 
presence of gravity). Similarly the superpotential may contain arbitrary powers of the The 
sugra Lagrangian turns out to depend only on a single function of 4>* and <pi, namely (•) G(<p* , <ft) — 
J{4>*,4>) + log\W\ 2 ; J = -3log(-K/3), where G (or J) is called the Kahler potential. Note G 
is invariant under (••) J — > J + h((j>) + h* ((/>*); W — > exp(—h)W . The sugra Lagrangian C may 
be written as C = Cb + £f.r" + C-f where Cb contains only bosonic fields, Cfk contains fermionic 
fields and covariant derivatives (supplying the fermionic kinetic energy terms), and Cf has fermionic 
fields but no covariant derivatives (see p| for details). Now regarding spontaneous susy breaking, for 
theories with local susy breaking the vacuum energy is no longer positive semidefmite. There are in 
particular the following possibilities. First at least one of the fields in the theory must have a vev not 
invariant under susy. Under certain assumptions one has then an F term generalization involving 
dW/dcj) + 4>*W 7^ or a D term generalization involving G l (T )j.^j ^ where G 1 = ((ft 1 )* + 
(l/W)dW/d<t>i (here (T )jj ~ generators of the gauge group in the appropriate representation). 
There are also other possibilities (cf. fig]]). 

For string theory both IIA and IIB are unsuitable to describe the real world but the heterotic 
string is perhaps tenable, with the extra 16 left mover dimensions providing the gauge group for 
the resulting 10-D theory (upon compactifying on a 16-D torus or variations of this). Further 
compactification of 6 dimensions is then still necessary leading at first to an N — 4 susy theory 
if toroidal compactification is used (which is unsuitable since N > 2 susy models are nonchiral). 
However orbifold compactifications will yield a 4-D theory having N = 1 susy and then some 
symmetry breaking must be induced. One can also compactify on Calabi-Yau (CY) manifolds and 
modular invariance can be achieved. We leave ||l2| , [l3| in what follows in order to concentrate on 
material more directly related to the projects at hand. 



d A xTrF mn F mn + (8.6) 



In 1143] one takes now an TV = 1 Yang-Mills (YM) action 



1 

8tt 



t d 4 x d 2 8TrW a W a 



Qyai 
'32tt 2 



• / <i'.rj r - ' ~ ~ ~ ' ' 
where (440 T = (^ym/^) + (47ri/<? 2 ). Here the W a are chiral spinor superfields of the form 



W a = -i\ a (y) + 9 a D(y) - -(a m a n 9) a (d m v n - d n v m ){y) + (69)a™ $ d m \P (y) (8.7) 

W a = -i\ a (y) + a D(y) + 9^af na F nn (y) - {ee)a m ^ a V m X^y) 
where F mn = (l/2)e mn P c 'F pq (dual field strength) and 

F mn = d m v n - d n v m + i[v m , v n \, m X p = d m X p + i[v m , X 13 ] (8.8) 

We omit the background considerations involving the Wess-Zumino (WZ) gauge etc. (cf. also [fl3|| ). 
For N = 2 susy one thinks of 8 and with D a — * D ai D a and / d 2 9 — * J d 2 6d 2 9, etc. Then an 
N = 2 chiral superfield is an N = 2 scalar superfield which is a singlet under global SU(2) and 
satisfies 

D & V(x, 6, 9, 9, 9) = 0; = (8.9) 

where e.g. 

D a = d a + 2ia^d m ; D & = -B & (8.10) 
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Set also ( XXX) y m = x m + iOa m + i9u 9. Then expanding an TV = 2 chiral superfield in powers 
of 9 the components are N = 1 chiral superfields. Thus 



.11) 



(W a is an A^ = 1 chiral spinor superfield). For N — 2 YM, if one forgets about renormalizability, 
there arises 



1 

47T 



d 2 9d 2 9TrT(^) 



.12) 



with T = (1/2)t* 2 and constraints on * of the form (XXX) (D aa D b a )^ = {D%D b6t )^ where a, b 
are SU(2) indices. Writing J C ' Q ($) = dF/d$ a and _F o6 = d 2 T /d$ a d<$> b the Lagrangian flS.12| ) can 
be written in terms of N — 1 superfields as 



47T 



i /" d 2 9T ab (<P)W aa W b a + J d 4 9(^e 2V ) a T a (<P) 



.13) 



where V will be clarified below. 



8.2 Soft susy breaking and spurion fields 



We go now to |61j, |62|, |63|, [147| and one works upon the foundation of []83| sketched in Section 7.1. The 
idea of soft susy breaking goes back to |pl[ for example and was developed in a form relevant here in 
j|, |5|, ||, 146 1 . From |^l], ^32|, 147] one extracts the following philosophical comments: Softly broken 



susy models offer the best phenomenological candidates to solve the hierarchy problem in grand 



unified theories. The spurion formalism of |81 provides a tool to generate soft susy breaking in a 
neat and controlled manner (i.e. no uncontrolled divergences arise). To illustrate the method, start 
from a susy Lagrangian L(& , $!,-••) with some set of chiral superfields, and single out a particular 
one, say $o. If you let this superfield acquire a constant vev along a given direction in superspace, 
such as e.g. < $o >= co + 9 2 Fo, it will induce soft susy breaking terms and a vacuum energy of 
order |F | 2 . Tur ning the argument around, you could promote any parameter in your Lagrangian 
to a chiral superfield, and then freeze it along a susy breaking direction in superspace giving a vev 
to its highest component (the F n terms below). In the embedding of the SW solution within the 
Toda-Whitham framework we have obtained an analytic dependence of the prepotential on some 
new parameters T n (the Whitham or slow times). Then these slow times can be interpreted as 
parameters of a non-supersymmetric family of theories by promoting them to be spurion superfields. 
In |5|, |(| |146[ this program was initiated with the scale parameter A and the masses of additional 
hypermultiplcts rrii as the the only sources for spurions. To deal with the times T n now one writes 

T n = T n T{ n ] u k = Tfu k ; a,{u k ,T n ) = T^^Uk, A = 1) + 0{T n>1 ); (8.14) 

hi = a i (u fc ,Ti,r„ > i = 0) = Tiai(u kl K = 1) = ai(u kl K = 7\) 
Then the Whitham times T n (or T n ) are promoted to spurion superfields via 

si = -ilog(A); s n = -if n ; Si = si + 8 2 F 1 - (8.15) 

Vi = \d x 9 2 ~9 2 - S n = s n + 9 2 F n - V n = \D n 9 2 ~9 2 

A = exp («Si); Si = — ; r = — + — K 2N - exp{2mr) (8.16 
N 2tt g z 



53 



(note the 9 in r has a different meaning from the Grassman 9 in the superfields). One also writes 
Hm+i,n+i = T™ TL m +i,n+i and in the manifold T„>i = the T n are dual to the Tln+i via 



/T ^ TTAT , dT iV- a? 7 
(DUAL) — — — — = — H 2 ] — 

dlog(A) m &T n 



T m =S„ 



- N -H 

— — — n-n+l 

mn 



(these are special cases of more general formulas below and in Section 7.1 - cf. ( 7.20| )). Further the 
T~L n +i are homogeneous combinations of the Casimir operators of the group; this means that one 
can parametrize soft susy breaking terms induced by all the Casimirs of the group and not just the 
quadratic one (associated to A) . In this way one extends to JV = the family of AT = 1 susy breaking 
terms first considered in ||. Note the Uk in the SW curve (7.1) can be written as Uk =< Ok > 
for Ok = (l/k)Tr<f) k + lower order terms (basic observables for a complex scalar field cf) as in [147] 
for example - see Section 10); the 'H m ,n are certain homogeneous polynomials in the Uk and the 



(Casimir) moduli hk of (7.1) refer basically to a background Toda dynamics (cf. [ 1 13 1 ) . Thus there 



are at least two points of view regarding the role of the T n : (1) The philosophy of (8.15) implies that 



T n or T n correspond to coupling constants while (2) The philosophy of looking at moduli dynamics 
of Uk or h n depending on T n puts them in the role of deformation parameters. We should probably 
always treat T n and ay in parallel, either as coupling constants or deformation parameters (this does 
not preclude treating a,j = aj(T n ) however). 



We indicate now some formulas arising from [ pl| , p3| which arc discussed further in 18E] in con- 
nection with [137, 145, 162]. Thus, in the notation of ]6l|, |62| ], one defines Hamiltonians H. m +i,n+i = 
'Hm+i.n+i with Hrn+i = 'Hm+1,2 (homogeneous polynomials in the Uk (or hk) via 



N 



Ti-m+l.n+l — -KeS 

mn 

and setting — dF/ds n write 



Re Soo (p m ' N dPl /N ) = H n +i, m+ v, Hn+i = ^Re Soo (p m ' N d\) (8.17) 



d 2 F 
daiCtj 



d 2 F 

ddidSn ' 



d 2 F 
ds m ds 

n 



There results 



n 



2tt 



TYl S m$n,7~(' m +1, n +1 

m > 2 m,n>2 



(8.18) 



(8.19) 



JL 

2im 



'Hn+i + i ^2 ms m H m+ i <n+ x 



A = JL 
1 2tt 



da 



2 , . dHn+i 

n>2 1 



n & dTin+1 ii 



27m da; 



An 



i = -2TlTjd Tij lo 9 e B (0\E); 



2rlTfd Tij logQ E {Q\r); r mn = - 2^Tj l d Tt] logQ E {Q\r) 



Here &e designates 

@b(0|t) = 0[a, /3](£V|t) = expliTTTij-niUj + itViTii — in ^""^ n l \ 

where V l = du 2 /da i and a = (0, • • • , 0) with (3 = (1/2, ■ • • , 1/2). 



.20) 
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ALERT! Henceforth we assume all Ti p , dfc, Ui, etc. have hats but we remove them for muta- 
tional convenience. 

Now the spurion superfield Si appears in the classical prepotential as (D) F = (N / n)SiH.2 and 
one obtains the microscopic Lagrangian by turning on the scalar and auxiliary components of S\. 
Next the remaining S n are included and one expands the prepotential around S2 = • • • = sw-i = 0; 
the D n and F n will be the soft susy breaking parameters (more on this below). The microscopic 
Lagrangian is then determined by 



yy^i N 
F = — > — o n n n +i + - — ; 
ir i — d n 2m 

(with s n = 0) and one is primarily interested in 



m,n>2 



SnSrnHn+1 



m+1 



.21) 



rred _ iv \ " 1 c nj 

J~ — > — ^niLn+l 

n * — ' n 

l 

which in fact is the relevant prepotential for Donaldson- Witten (DW) theory. Note 



q2 jrred 



dS m dS n 



2N 2 dH m+ i dH n+1 1 



irimn den 



daj 



d Tij loge E (0\T) 



.22) 



.23) 



are essentially the contact terms of [137] (cf. Section 9.1). Expanding ( 8.22 ) in superspace one has 
a microscopic Lagrangian and this gives an exact effective potential at leading order for the J\f = 
theory, allowing one to determine the vacuum structure. Detailed calculations for SU(3) theory are 
given in 61 62]| . 



So are we dealing with coupling constants T n or deformation parameters? One answer is "both" 
and we refer to Section 9 for more on this. The spurion variables parametrize deformations of the SW 
differential and the T mn and r™ of Q8.19D have nice transformation properties under Sp(2(N — 1), Z); 
the s n behave like the ay in many ways. Promotion of T n — > S n for susy breaking should however 
correspond to a coupling constant role; the T n are parameters of a non-susy family of theories; thus 
role of s n is parallel to a?. The prepotential determines the Lagrangian via Ta-, Fab, Fabc, the- 
D n and F n , A, ip gluinos, <j) — scalar component of M = 2 superfield, etc. Thus C — Lki n + Li n t 
with 

1 



4-7T 



3 [(V^KV^r + i(V^)ia^ b FZ- 



.24) 



47T 



iF ah X a ^{V^X) b - - A F ah (F a F h ^ + iF«F b n 



F AB F A {F*) B - l -F abC ((VV)0n C + {X a \ h )F c + tV2(^ a X b ) 



D 



c 



1 



-F AB D A D l 



W (KfbcD b F c + V2{(^X) a F^ b - (j~\) a F a } 



Here A and ip are t ne gluinos and <f> is the scalar component of the J\f — 2 vector superfield. 
The f? c are structure constants of the Lie algebra. The indices a, b, c, ■■■ belong to the adjoint 
representation of SU(N) and are raised and lowered with the invariant metric. Indices A, B, ■ ■■ run 
over both indices in the adjoint and over the slow times (to, n, ■ ■ •). Since all spurions corresponding 
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to higher Casimirs are purely auxiliary superfields the Lagrangian can be simplified as follows. Set 
° a = -(KLsY lC (( bclaSS )? D rn + M(g0tfLf a )) with F a = - (b^ ass ) ac {b class )™F m where the classical 
matrix of couplings 65df s is defined via b class = (\/An)T class where 



T(! abi 'clas 



0; (T clas 



N 



nim dcf> a nim 



-Tr 



(8.25) 



where the dots denote the derivative with respect to 4> a of lower order Casimir operators. This leads 
to 

1 



B 



,,„.... , F m F r : + -D m D„ ) + f b e c (b class r(b class )aeD, 



.26) 



i p){ —class \ m 
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9 RENORMALIZATION 



We extract here from p6| where a number of add i tiona l topics concerning renormalization also 
appear. In particular we omit here the work in 1 105 , |l06| (sketched in |26[| and partially subsumed 
in the formulation of Section 7.1) and other work of various authors on the Zamolodchikov C theorem. 
The formulas in Section 7.1 involving derivatives of the prepotential all have some connection to 
renormalization of course and to indicate this briefly we refer to the original elliptic curve situation 
(cf. ||, H HI HI HH for example). Thus consider (CC) y 2 = (A - A 2 )(A + A 2 )(A - u) 
for example with a = (V2/vr) ff 2 [(A - u)/(A 2 - A 4 )] 1 / 2 ^. One will have then for F ~ F sw 



2F = aF, 



2iu 2iu 

; A_Fa = = Snibiu 

7T 7T 



(9.1) 



where b\ = 1/Air 2 is the coefficient of the 1-loop beta function. This is the only renormalization term 
here and (in the more enlightened notation of Section 7.1) KF^ — TiF^ shows that an important 
role of the Whitham times is to restore the homogeneity of F w which can be disturbed by renor- 
malization (cf. p6l |2^, Q for more on this - and see below for more details about renormalization). 
We remark also that beta functions for this situation are often defined via 



13(t) = Ad A r| u=c ; (3 a (r) = Ad A r| a=c 
where r is the curve modulus (e.g. t = F™ - cf. O, ||, ||). 



(9.2) 



In any event renormalization is a venerable subject and we make no attempt to survey it here (for 

m m n n n ~ 



105, 106 ; 133, 155, 18C, 182 



In 
the 



renormalization in susy gauge theories see e.g 

particular there are various geometrical ideas which can be introduced in the space of theories 
space of coupling constants (cf. here || || ||, || ||, IM HI- [HI HI PI)' We extract here now 
mainly from fl42| where it is argued that RG (= renormalization group) flow can be interpreted as 
a Hamiltonian vector flow on a phase space which consists of the couplings of the theory and their 
conjugate "momenta" , which are the vacuum expectation values of the corresponding composite 
operators. For theories with massive couplings the identity operator plays a central role and its 
associated coupling gives rise to a potential in the flow equations. The evolution of any quantity 
under RG flow can be obtained from its Poisson bracket with the Hamiltonian. Ward identities can 
be represented as constants of motion which act as symmetry generators on the phase space via 
the Poisson bracket structure. For moduli g regarded as coupling constants one could obtain beta 
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functions via K n (dg/dK n ) = d n g for T n — log(n n ) and d/dT n = K n (d/dK n ). Whitham dynamics 
on moduli spaces corresponds then to RG flows and one obtains e.g. dhk/dT n or diik/dT n which 
could represent beta functions j3^. Thus following |^6|, ^2) (revised for |3(]]) we may consider the 
moduli space of Uk or h n as a coupling constant space M. of elements g a (space of theories) and 
interpret RG flows as a Hamilonian vector flow on a phase space T*A4 (in this context a 3 should 
also be regarded as a deformation parameter but we will ignore it here). Take T = log(n) and set 
P a {g) = nd K g a (note, corresponding to A = exp{is\) in ( |8.15 ) with is\ — log(A) one would obtain 



(i a ~ Ad\g a and recall that (3 — Ad\r is a standard beta function). We have now a tangent bundle 
TM. ~ (g a ,(3 a ) with T*M ~ (g a ,4> a ) where cj) a = dw(g : t)/dg a for some free energy W = —log(Z) 
with W ~ J wd D x. Here Z could correspond to J T>cj)exp(—S((f))) and 1 = J V(j)exp[—S(<p) + W] 
implies dW —< dS >= j V(j)dS{(j))exp[~S + W] (and S ~ / Cd D x for some Lagrangian so an 
underlying D-dimensional space is envisioned). The approach here of |42[ is field theoretic and 
modifications are perhaps indicated for the N = 2 susy YM theory; thus the constructions are 
heuristic. Now one can display a Hamiltonian (HH) H(g,<j>) = Y.^is)^ + /3 r (5,r)</> r which 
governs the RG evolution of (g a , <j> a ) via 



dg a _ dH 
~dT~W n , 



d4a _ dH 
dT dg a 



(9.3) 



Here Y (cosmological constant) is a coupling associated with the identity I whose conjugate momen- 
tum is the expectation value of the identity - in fact one can take heuristically 

/=T; p r (g,Y) = ^ = -DY + U r (g); <f> r = *° (9.4) 

where T = log(n) could refer to any T n (~ n n ). One has then a symplectic structure and a Hamilton- 
Jacobi (HJ) equation 

g + «( 9 ,^)=0=g + yj /3 " W . + /S r *- (9-5) 

For j3 a = f3 a (g,T) one could take T as an additional coupling and work on M. = (g a ,Y,T) with 
f3 T = 1 and 4> T = d T w = -H(g, <j>, T) (T - log(n) implies nd K T = f3 T = 1). 



Now identify w ~ F + Yn so <Pt = wt = Ft + YDk so that (9.5) says e.g. (with g a ~ h a ) 



dF v-^ dhk dF r n , 

(note /3 r 4>r ~ (-DY + U r )n D DYn D + U r n D and the DYn D term cancels). We know by 

Whitham dynamics that hk — hk(ctj,T n ) satisfies some homogeneity equations (A) Y] a Adhk/ dai)- 



YT n d n h k = so a typical equation like (EE) 2F = J2a 3 (dF/da 3 ) +J2 T nd n F (cf. (|7.14|) ) plus 
(9.6) for T ~ T n implies (we think now of w n = F + Yk® for fixed F while U T is also held fixed for 



different T n 



= T ndnF -J2T n (-d n F - U t k») = 2 T ndnF + (V T n n») U T (9.7) 
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This seems to say (for T n exp(DT n ) = T) 

F -J2T n d n F = \tU 1 



(9.8) 



where F — F(T n ,aj) and h k = h k ( T n , ct j) with T n and ctj independent allow us to set ctj = ctj(hk) 
for T n fixed. In particular from (7.14) and ( |9.6| ) with say T n — <5„ ; i and T% ~ A (there is an 
implicit switch here to the notation of ( |3.14 )) we see that ( [7.14[) sa ys 2F — Ad\F + ^ aj(dF/da,j) 
(and F = F sw now) while (JD| says (for T - Zot?(A) as in fl8.15|)) and no other T n (C) Ad A F + 
J2Hdh k /dA){dF/dh k ) + U^JP = 0. This means (cf. (O) with more variables inserted) 



2F 



' <9a 7 - 



A9 A F 



2iu2 



-U L A 



(9.9) 



v 8F dh k _ T D v . dF 

Note ?/2 ~ /i2 ~ ^2 (cf. below for more such notation) is a Hamiltonian and we see that it has the 
form (HH) (after shifting W ~ F + TA^ 1 ). In addition t/ r is determined as 



DttT 



2i 



-u 2 



dF 
dh k 



One notes that HJ equations reminiscent of (9.5) appear in [137 in the form 



~dt 



-H [ a,—,t 
oa 



(9.10) 



(9.11) 



in connection with Lagrangian submanifolds (with H somewhat unclear) and this is surely connected 
to formulas such as (6.21) referring to Whitham flows. 



Thus our heuristic picture of renormalization based on j2(| [42| seems suggestive at least, modulo 
various concerns over terms like T, W — J wdPx with w ~ F + Tk d for t ~ log(n), etc. We remark 
that in one has a dictionary of correspondence between QFT or statistical mechanics with RG 
and classical mechanics. This has the form 



QFT ~ Statistical Mechanics 
Couplings ~ g a (t) 
Beta functions ~ (3 a (t) 
vev's ~ ^ a (t) 
Bare couplings ~ > <%) 
Generating functional ~ w(g(t), go,t) 
H(g,^t)=p a (g,t)cj> a 

3A 



/3 a = 
(7(. 9 ,<)=/3 A + ^A 



Classical Mechanics 
Coordinates ~ q a (i) 
Velocities ~ q a (i) 
Momenta ~ p a (t) 
Initial values ~ (9o>Pa) 
Action ~ S(q(t) ,qo,t) 
H = ^g ab (q)PaPb + U(q,t) 

■a _ dH_. ■ _ _ 

tf(g,f) 
^ = -<fl7 



9H 



(9.12) 



<9w j- \ 







w, +H(g(t),g^. 

No explicit k dependence in f3 
Anomalous dimensions 
RG invariant {9, H} = 



5, 



)>§>*) 



= 



Conservative system 
Pseudo — forces (Coriolis) 
Constant of motion {9,H} = 
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In [177 1 one also finds a dictionary of comparisons between QFT and magnetic systems based on 
the effective action and we expand on this idea of effective action as follows. 



REMARK 9.1. In fact the effective action is isolated by A. Morozov in [ 165 1 as the source of 



integrability (cf. also 



113; 166 1 ) . Indeed effective action of the form 



exp[S eff {t\<f>)} = Z{t\4>) = / T></>exp[S(t\<t>)] 



(9.13) 



in matrix models for example will correspond to a tau function of integrable systems such as KP 
and the "time" variables t can be thought of as coupling constants. Further investigation of related 
generalized Kontsevich models (GKM) leads to effective action involving both KP times t n along 
with "Whitham" times T k which have a rather different origin. Thus for D — n 2 , modulo a few 
details (C is defined below) 

Zgkm (L\V P+1 ) =C J d D Xexp[Tr(-V p+1 (X)+XL)} 

where L is an n x n Hermitian matrix, V v +i is a polynomial of degree p + 1, and W l 



(9.14) 



is a polynomial of degree p. Further C in ( 3.14 ) is a prefactor used to cancel the quasiclassical 
contribution to the integral around the saddle point X = A, namely, modulo inessential details, 



C = exp[Tr{V p+1 {K)) - Tr(A^ +1 (A))]de^ 2 [d 2 ^ +1 (A)] 



(9.15) 



The time variables are introduced in order to parametrize the L dependence and the shape of Vp+i 
via 

T k = ^Tr(A- k y 7 f k = \Tr{k- k ); 
k k 

P 



(9.16) 



L = W P (A) = A"; t k = fc(p _ k) Res p W^ k ^(fi)d^ 



Then one can write 



Zgkm = exp[-J r p (T' fe |t n )]Tp(T fc + t k ) 
where r p is a p-reduced KP tau function and 



J D 



\Y, A < 



[t)(fi + t t )(f 3 +tj); A l3 = Re Soo W^(X)dW{ /p {\) 



(9.17) 



(9.18) 



is a quasiclassical (Whitham) tau function (or the logarithm thereof). In addition such a GKM 
prepotential is a natural genus zero component in the SW prepotential (cf. Q and compare ( 9.18) ) 
with (7.23) for example). 



REMARK 9.2. It is possible now to apprehend how the Whitham times can play two appar- 
ently different roles, namely (1) coupling constants, and (2) deformation parameters. We follow 
the Russian school of Gorsky, Marshakov, Mironov, and Morozo v (fj2| |83| , |88|, ^9[ [M)|, pT6]|l65] [l66| , 
15 4 154 1) and especially the latter as in ]165] , |166| (cf. also [113, 137Q. The comments to follow are 
mainly physical and heuristic and the matter seems to be summarized in a statement from |165| , 
namely: The time variables (for SW theory), associated with the low energy correlators (i.e. the 
renormalized coupling constants) are Whitham times (the deformations of symplectic structure). 
Thus both roles (1) and (2) appear but further clarification seems required and is possible. One 
relevant theme here is described in [113, 165 1 roughly as follows. Given a classical dynamical system 
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one can think of two ways to proceed after exact action-angle variables are somehow found. One can 
quantize the system or alternatively one can average over fast fluctuations of angle variables and get 
some effective slow dynamics on the space of integrals of motion (Whitham dynamics). Although 
seemingly different, these are exactly the same problems, at least in the first approximation (non- 
linear WKB). Basically the reason is that quantum wave functions appear from averaging along the 
classical trajectories - very much in the spirit of ergodicity theorems. In string theory, the classical 
system in question arises after some first quantized problem is exactly solved with its effective ac- 
tion (generating function of all the correlators in the given background field) being a tau function 
of some underlying loop-group symmetry for example. The two above mentioned problems concern 
deformation of classical into quantum symmetry and renormalization group flow to the low energy 
(topological) field theory. The effective action arising after averaging over fast fluctuations (at the 
end point of RG flow) is somewhat different from the original one (which is a generating functional 
of all the matrix elements of some group) ; the "quasiclassical" tau function at the present moment 
does not have any nice group theoretical interpretation. The general principle is in any case that 
the Whitham method is essentially the same as quantization, but with a considerable change in the 
nature of the variables; the quantized model lives on the moduli space (the one of zero-modes or 
collective coordinates), and not on the original configuration space. 

Note that in higher dimensional field theories the functional integrals depend on the normalization 
point [i (IR cutoff) and effective actions describe the effective dynamics of excitations with wave- 
lengths exceeding [iT 1 . The low energy effective action arises when \i — > and only a finite number 
of excitations (zero modes of massless fields) remain relevant. Such low energy effective actions are 
pertinent for universality classes and one could say that AT = 2 susy YM models belong to the uni- 
versality class of which the simplest examples are (0 + f )-dimcnsional integrable systems. We could 
supplement the above comments with further remarks based on the fundamental paper |8^| (some 
such extractions appear in |2(J). Let us rather try to summarize matters in the following manner, 
based on |B2, 



1 65 1 . First think of a field theory with partition function Z(t\4>) = f, T><j) exp[S (t\(f>)] 



( 9.13| ) . The dynamics in space time is replaced by the effective dynamics in the space of coupling 



constants g , 8, and tj via Ward identities which generate Virasoro conditons and imply that Z is 
the tau function of some KP-Toda type theory. The parameter space is therefore a spectral curve 
for such a theory and the family of vacua (~ <fio) is associated with the family of spectral curves 
(i.e. with the moduli space). Now the averaging process or passage to Whitham level corresponds to 
(A) Quantization of the effective dynamics corresponding to some sort of renormalization process 
creating renormalized coupling constants Ti for example, and (B) Creating RG type slow dynamics 
on the Casimirs hk of the KP-Toda theory; then via a 1 — 1 map ihk) — > (u m ) to moduli one deals 
with moduli as coupling constants while the T n are deformation parameters. In this spirit it seems 



that in |6lj, 62, |63| both features are used. Promoting T n to spurion superfields corresponds to 
a coupling constant role for the T„ (as does treating the prepotential as a generating function of 
correlators) while dealing with the derivative of the prepotential with respect to the T n in the spirit 
of renormalization involves a deformation parameter aspect. 

REMARK 9.3. We remark that there is also a brane picture involving RG flows and Whitham 



times (cf. |85|, |86j, |87], |39| |92j, |9CJ, 116] for example) in which Whitham dynamics arises from the 
motion of certain D — 4 branes. This dynamics generates conditions for the approximate invariance 
of the spectral curve under RG flow and provides the validity of the classical equations of motion 
in the unperturbed theory if suitable first order perturbation is allowed. It is also conjectured that 
both Hitchin spin chains and Whitham theories can be regarded as RG equations (Hitchin times 
are to be identified with the space RG scale - t ~ log(r) - and the fast Toda-Calogero system - 
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motion of certain D — branes - corresponds to RG flows on a hidden Higgs branch, providing susy 
invariant renormalization for the nonperturbative effects, while the Whitham system is operating on 
the Coulomb branch). In this spirit the spectral curve is the RG invariant and the very meaning of 
integrability is to provide the regularization of nonperturbative effects consistent with the RG flows. 



9.1 Contact terms 

Another arena where Whitham times play an important role involves the structure of contact terms 
in the topological twisted N = 2 susy gauge theory on a 4-manifold X (cf. |6l[ ffi ^37j pSl, |l47| , 
|162| , |l89fl ). We will not try to cover the background here. An earlier version of this section (written 
in 1998) made some attempt at this but it would have to be considerably enlarged and revised to 
be instructive now. Since enlarged versions already exist in the references above there seems to be 



no point in an inadequate sketch. Thus we follow mainly [145, 147, 189 1 here for the SU(N) theory 
(as in Sections 7 and 8) and go directly to the u-plane integral (for b 2 + = 1) 



Mcoulomb 



[da da]A(u k ) x B(uk) a exp (^2 PkUk + ^ 



, fkfrnTk,' 



(9.19) 



Here S S H 2 {X, Z) and the T^ m are so called contact terms. We omit discussion of the other terms. 
The contact terms can be derived via blowup procedures and this introduces the tau function of a 
periodic Toda lattice. In fact one can show that (notation as in Sections 7 and 8) 



Ti 



km ( d 2 T r 



fe+l,m+l 



4/V 2 \dT k dT„ 



(9.20) 



T n>2 =0 



(cf. (3.23)) where T r is defined in ( |8.22 ). The key idea is that the slow times T k are dual to the 
Ttk+i hi t he su bspace T„> 2 = (cf. (DUAL) after ( B.16 )); note also Hk+i — U k+i + Oju k)- Further 



the fk in ( 9.19 ) are proportional to the fast Toda times t^- The development in 145 , 147 , 18£] shows 
how this is all a very natural development. 



10 WHITHAM, WDVV, and PICARD-FUCHS 

We add a few comments now relating Whitham times to the flat times of Frobenius manifold theory 
and the WDVV equations. This involves connections to TFT, Landau- Ginzburg (LG) models and 
Hurwitz spaces, and Picard-Fuchs (PF) equations. 



10.1 ADE and LG approach 

Connections of TFT, ADE, and LG models abound 
for N = 2 susy YM we go to fll| (cf. also @, | 
Thus one evaluates integrals ai — § A Xsw and af = § B Xsw using Picard-Fuchs (PF) equations. 
One considers Pn(u,Xi) = det(x — Qr) where R ~ an irreducible representation of G and is a 
representation matrix. Let u,; (1 < i < r) be Casimirs built from of degree e 2 ; -I- 1 where is the 



% H 11 H Hi' jgj . Pl : Hi) and 

First we extract from [112| as in pfl. 



exponent of G (see below). In particular 



Ml 



h where h is the dual Coxeter number of G (h 



quadratic Casimir and u r ~ top Casimir of degree 
r + 1 for A r ). The quantum SW curve is then 



P R (x, z, Ui) = P R x, m + Si 



z + 



(10.1) 



Gl 



where fi 2 — A 2h with A ~ the dynamical scale and the Ui are considered as gauge invariant moduli 
parameters in the Coulomb branch. This curve is viewed as a multisheeted foliation x(z) over CP 1 
and the SW differential is Asw = x{dz/z). The physics of N = 2 YM is described generally by 
a complex rank(G) dimensional subvariety of the Jacobian which is a special Prym variety (cf. 
|H |l~54] l). Now one writes ( |10.1[ ) in the form 



z H V u r = Wq(x, u\ , 







(10.2) 



For the fundamental representation of A r for example one has (we have organized the indexing to 



conform more closely to |4q, 83L 154, 168 



W r A +1 = x r+1 



UiX T 1 — • • • — u r ^ix; 



(10.3) 



and setting (SP) Wq (x, u\, ■ ■ ■ , u r ) — Wq(x, Ui, ■ ■ ■ , u 



r-l) 



; r it follows that W^ 1 is the funda- 
mental LG superpotentials for A r type topological minimal models (cf. also |59| |53|, [35|). The 
can be thought of as coordinates on the space of TFT. We will concentrate on A r but D r and other 
groups are discussed in 1 112]. For comparison to we recall (cf. [^6|) that for a pure SU(N) susy 
YM theory 

det NxN [L(w) — A] = 0; P(A) = A N ' 



N 



(10.4) 



P(A) = A JV -^^A w - fc = n(A-A J ); v k = (-l) k E A n ---A !t 

2 1 ii<-<i h 

Here the Vk are Schur polynomials of hk — (1 /fc) X3i ^ ne formula (A) log{A~ N P{\)) = 

— ^ /c (ft,fc/A' £ ); there are g = N — 1 moduli Uk and we refer to [ 1 1 3 1 for the Lax operator L. Thus 
v = 1, vi =0, v 2 = h 2 , v 3 = h 3 , V4 — I14 — (l/2)h 2 , etc. (hi — for SU(Nj). One also has the 
representation 

' P 2 (A)-4A 2Ar ; y = A N 1 



V 



w 



(10.5) 



giving a two fold covering of the punctured Riemann sphere with parameter A. Thus r + 1 ~ N ~ 
h, g = N—l, A ~ x, and ~ Vk+i with w r ~ ujv while P ~ W'^ 1 = where j4at_i ~ SU(N) 

or SL(N). To clarify the curve correspondence write 



w = z + /i 2 z ; it) = P(x); y = z — prz 1 ; y 2 = P 2 — 4/i 2 



(10.6) 



with /i = A leading to the equivalence of curves y = P — 4A . Recall also from Section 
7 (7.4), etc., that dSsw — X(dw/w) = X(dP/y) = \(dy/P). For the moment we ignore relations 
between w and w and refer to 154 for discussion of the algebraic geometry picture for the w 
parametrization. 

Now in 2-D TFT of LG type A r with superpotential (SP) the flat time coordinates for the 
moduli space are given via 



T l = c l fdxW§(x,u) e '/ h (i = l,.-.,r) 



(10.7) 



(the &i will be discussed later and recall h — r + 1). Note that no Riemann surface is involved 
here since the integral corresponds to a residue calculation for x ~ p at 00 (cf. Remark 10.2). 
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Thus the TFT here depends only on Wq (x, u) with x a formal parameter. These are residue 
calculations for times Tj (with no reference to Whitham theory) which will be polynomials in the uj 
(the normalization constants c,; are specified below. One defines primary fields 



kRl ^ _ dW*(x,u) 



(10.8) 



where 0^ = 1 is the identity ~ puncture operator. The one point functions of the gravitational 
descendents a n ((j)f) (cf. [fl E39l p3|) are evaluated via 



< a n (4>?) >= b n ,i ^2 VH f W§(x, up/V+n+i ( n = o, 1, . . •) 



for certain constants b n> i (cf. |112{ for details). The topological metric rjij is given by 

d 2 



r/ij =< <j>?<j>fP >= b 



'o.r 



dTidT; 



W§(x,u) 1+ W 



(10.9) 



(10.10) 



and rjij = <J e4+e . ^ can be obtained by adjustment of c; and b n ^. The primary fields generate the 
closed operator algebra 



where 



<j>f{x)ct>f{x) = J2 C^(T)^(x) + Q^x)d x W§(x) 
i 

d 2 Wg{x) 



dTjdTi 



(10.11) 



(10.12) 



(for details on ( |10.11[ ) - ( |10.12[ ) we refer to [fl [39] pO] |66|, |67|, |94j, p0l|). Note again x ~ p is a formal 
parameter and X ~ T r would be a natural identification. At this point it is not evident (nor is it 
true) that the times Tj of ( |10.7[ ) correspond to the Whitham times of Section 7 based on the Toda 
curve y 2 = P 2 - 4A 2JV . Indeed we note from ( Og) a nd e.g. ( fT5l| ), that T k = -(l/k)Res £=0 £ k dS is 
more or less tautological whereas a prescription ( 10.7 ) determines the Tk as functions of Uj. Now the 
structure constants are independent of R since cy^ = cf.-r^fc is given via Cijk(T) =< 



> 



which are topologically invariant physical observables. In 2-D TFT one then has a free energy F 
such that djk = d^F/dTidTjdTu (cf. § g || ||). 



Now for the Picard-Fuchs (PF) equations, the SW differential A 



sw 



A = (xdz/z) can be written 



as (cf. p6l ||) X SW = [xd x W/y/W 2 - Afi 2 ]dx (for W ~ Wq) and one has then, for fi fixed 



9A 



sir 



y/W 2 - A/i 2 



dx + d 



dW 



^JW 2 - 4^ 2 



(10.13) 



(total derivative terms will then be suppressed). Suppose W is quasihomogeneous leading to 



dW 



qiT,— = hW (i.e. W{tx,t^Ti) = t h W(x,Ti)) 



(ft 



1 is the degree of Tj). Then 



hWdx 



y/W 2 - 4/i 2 



(10.14) 



(10.15) 



G3 



and applying the Euler derivative E QjTj (d / dTj) to both sides yields (after some calculation) 



sw 



d 2 Xsw 



dT? 



The calculation goes as follows; first 

52 'I- 1 '< } ' X - (52 <l, l ,'>j : 52 '/' / '''' A = hdx 52 '<' ' ''' 



= 



w 



52 'l-'h ' - 'j*:; + 52 ~ 1 i7 .; A ; " 52 ' - X ' ' A 



^W 2 - Ajj? 
hWdx 



(10.16) 



(10.17) 



^W 2 - 4^ 2 



-hdx 



52 



^Jw 2 - v " V) 3/2 



= X <i,(i,'l ,1 ,\,, + £ ?j(gj - 2)TjAj + A = focte 
Next note that 



W 



hd 



oW 



^JW 2 - Ay? 



= hdx 



W 



A\x 2 (hW - xW x ) 

v / M/ 2_ 4/i 2 + ( ^2_ V)3 /3 

Ay?xW x 



^JW 2 - Ay 2 (l^ 2 - V) 3 / 2 



(10.18) 



so the right side of ( |10.17| ) can be written as AW h 2 dx / (W 2 - 4/i 2 ) 3 / 2 + M [xW/^W 2 - 4/i 2 
We note also that by degree counting (see below for more on this) 



d 
dT r 



52 



du.i d 
dT r dm 



d 

8u r 



dW 



= 1 



(10.19) 



since dW/du r = -1 (note here du k /8T r = -S kr => ( |l0.19| )) and thus from ( |10.13|) {W r = 1, W 
0) 

<9A cfe , / x \ 



dT r ^w 2 - An 2 
Wdx + d( dr 



^W 2 - An 2 



(10.20) 



5 2 A 



dT 2 (W 2 - 4/i 2 ) 3 / 2 



Wdx 



^JW 2 - Ay 2 J (W 2 -Afi 2 r/2 



Consequently the right side of ( 10.17 ) is equivalent to 4/i 2 h 2 d 2 A / dT 2 and ( 10.17 ) becomes 



2,2<9 2 A d 2 X 



OX 



dTidTj 



£<&(<& -2)^— + A = 



(10.21) 



which is equivalent to ( 10.16| ). We note that the second term in ( 10.16] ) represents the scaling 
violation due to /i 2 = A 2/l since ( 10.16[ ) reduces to the scaling relation for Xsw m the classical 
limit fi 2 — > 0. Note that Xsw(Ti, /i) is of degree one (equal to the mass dimension) which implies 
(444) (Ei qtT t {d/dT t ) + h(j,(d/dfi) - 1) Xsw = (see Remark 10.1 below) from which ( pU6| ) can 
also be obtained. In this respect we note that 



AxfiW x dx 

f = ( W 2 _ ^2)3/2 



(52<liTidi-l)X = -h^d fi X = - 



Ah/i 2 xW x dx 
(W 2 - Ap?fl 2 



(10.22) 
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Then, using (|10.18[) 



(^q i Tid i -l)X = -hfid tl X 
xW \ 



hd 



4hfi 2 xW x dx 

(w 2 - V) 3 / 2 

Whdx \ 



^JW 1 - Afi 2 J \ \/W 2 - Aii 2 



and this implies, via ( 10.2C ) 

{^q&idi-lf X = -hyd la, 



xW 



\ 4h 2 fi 2 Wdx 2 2 d 2 X 

r M /W^V J + (W^ 2 - V)3/2 ~ ^ df 2 



REMARK 10.1. In connection with scaing we note from 

W(tx, t n+1 u n ) = t N x N ~ Ul t 2 (xt) N - 2 t N u N -i = t N W{x, tin) 



that (corresponding to (10.14)) 



xd x W + + l)u n I —J = NW 



(10.23) 



(10.24) 



(10.25) 



(10.26) 



Alternatively we can write W(x,u n ) = x N — U2X N 2 — ■ ■ ■ — un with W(tx,t n u n ) — t N W(x,u n ) 
leading to [xd x + J2 nu n (d/du n )]W = NW. Similarly we see that W x = Nx N ~* —(N— 2)uix 



„N-3 



■ ■ ■ — uat_2 implies W x (tx, t n+1 u n ) = t N 1 W x (x, u n ) and consequently for A = xW x dxj \JW 2 — Ay? 
one gets 

/ ti 1 1 n \ txt ^V x tdx , , 

X(tx,t ^ u n ,t fi) = U7N „ n ., 9N ? , 1/9 = tX{x,u n ,fi) 



(fN W 2_ U 2N lJ 2y/2 

?d x + J^i 71 + 1 ) u ™q^~ + N ^ ) A = A 



Next from (10.7) we have (x — > t£) 



Ti(t n+1 u n ) = Cl jdxW{t n+1 u n ,xT' N = 

/dT 
d£W(u n ,0 = t qi Ti(u n ) => Y,(n + 1 ) u "q^- = H T i 



(10.27) 



(10.28) 



Now t o confirm Q10.14Q wri te W( tx,t^T if i) ~ W(tx, t n+1 u n , fi) = t N W(x,Ti, fi) (via ( |10.25| ) and 
( |l0~28l )). Further note (via ( gOjg ) 



X(tx,t 9i T h r fi) ~ X(tx,t n+L u n ,t N fi) =tX(x,Ti,fi) 



(xd x + <li T i d i + Nfidf. 



A = A 



(10.29) 



Let us try now to de rive t he relation (444), namely q{Tidi + Nfid^ — 1)A = 0. Thus, (444) = 
( |10.22| ) and thence ( |10.23| ) which means 



1 A 



^/W 2 - 4^ 2 



(10.30) 
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and this corresponds to (10.15) which we know to be true. This shows only however that an 
integrated form of (444) is valid (i.e. (444) is valid modulo hd(xW/ y/W 2 — 4/i 2 ). 



Another set of differential equations for Xsw is obtained using ( 10.11 ), to wit 



o 2 



dT t dT, 



■Xsw = Y. c U T )^K&T r Xsw 



(10.31) 



To see how this arises one writes from (10.11), (10.12), and (10.13) (using (10 



Wi 3 dx 



,Wdx 



d 2 x 

dT t dT 3 ~ ^J W 2 - 4 M 2 ' {W 2 - 4 M 2 )3/2 - 
dQjj (EC&k + QijWs) Wdx 



y/W 2 - Afi 2 



(W 2 - 4^ 2 ) 3 /2 



C*W k Wdx 



Then observe that from ( 10.13| ) 




(W 2 - V) 3 / 2 dT k dT r 



(10.32) 



(10.33) 



resulting in Ql0.31fl . Then the PF equations (based on J10.24Q and ( |10.31[ )) for the SW period 
integrals Id = ~§ Xsw are nothing but the Gauss-Manin differential equations for period integrals ex- 
pressed in the flat coordinates of topological LG models. These can be converted into u k parameters 
(where du k /dT r = —8 kr ) as 



qiUi 



d 

dui 



i n-4^ 2 /i 2 



d 2 u 

du 2 



(10.34) 



where 



d 2 u ^ A . . d 2 u ^ _ . . on 

+ ^ Aijk (u) + ^ B ijk (u) — 

J 1 ' I 



dT m dT n du k ^ 



i dui duj &Ti 



du k du r 

B ijk {u) 



£ 



d 2 T n du k 
duiduj dT n 



(10.35) 



which are all polynomials in tij. One can emphasize that the PF equations in 4-D N — 2 YM are 
then essentially governed by the data in 2-D topological LG models. 



10.2 Frobenius algebras and manifolds 

We sketch here very briefly and somewhat incompletely some basic material on Frobenius algebras 
(FA) and Frobenius manifolds (FM) following |5^, [55| (cf. also [125, 143| , 144 ) with special emphasis 
on LG models and TFT. This is not meant to be complete in any sense but will lead to a better 
understanding of Section 10.1 and serve simultaneously as a prelude to WDVV. Generally one is 
looking for a function F(ti, • • • , t n ) such that c Q/ 3 7 = d 3 F / dt a dt@ dt 1 satisfying (C) rj a p — ciaplt) is a 
constant nondegenerate matrix with rf 1 ^ = {r} a p)~ 1 , (D) c^g = "q ie c ta (){t) determines a structure of 



GO 



associative algebra A t : e a ■ ep — c^e 7 where e\, ■ ■ ■ , e n is a basis of R" with e\ ~ unity via c^ Q = ^ , 
and (E) F(c d H l , ■ ■ • ,c d "t") = c^i^t 1 , • • • , i") which corresponds to £ B F = E a d a F = d F ■ F 
for £7 = E a d a with i? Q = d a i Q here (we use t k ~ tfc and later, for a certain LG model as in 



Remark 10.2, t k ~ T k - cf. ( |10.41[ )). In @ one looks at di = 1 and physics notation involves 
d a = 1 — q a , d,F = 3 — d, q n = d and g Q + q n -a+i = d. The associativity condition in (D) reads as 
(WDVV equations) 

d 3 F Am d 3 F 8 3 F Am 8 3 F 

dt a dtPdt xV dtidt s d» ~ dPdt^dt^ dt a dt s dt u ( ' ' 

Next one defines a (commutative) FA with identity e via a multiplication (F) (a, b) — >< a, 6 > 
with < ab, c >—< a, 6c >. Here if u> € A* is defined by w(a) =< e, a > then we have < a, b >= tv(ab). 
The algebra ^4 is called semisimple (ss) if it contains no nilpotent a (a m = 0). Given a family A t 
of FA one often identifies A t with the tangent space T t M at t to a manifold M (t S Af). M is 
a Frobenius manifold (FM) if there is a FA structure on T t M such that (G) < , > determines a 
flat metric on M, (H) e is covariantly constant for the Levi-Civita (LC) connection V based on 
< , >, i.e. Ve = 0, (I) If c(it, v, w) =< u ■ v,w > then (V 2 c)(it, w, w) should be symmetric in the 
vector fields (u,v,w, z), and (J) There is an Euler vector field E such that V(V-E) = and the 
corresponding one parameter group of diffcomorphisms acts by conformal transformation on < , > 
and by rescaling on T t M . One shows that solutions of WDVV with d\ ^ are characterized by the 
FM structure (d a ~ d/dt a , £ E as in (E)) 

d a -dp = c^(t)9 7 ; < d a , dp >= riap; e x = d\ (10.37) 

and CeF — dpF + A a pt a t^ + B a t a + c (where the extra terms in CeF can be killed when dp ^ 
0, dF — d a ^ 0, and dF — d a — dp ^ 0). 

REMARK 10.2. The case of interest here is based on M = {W(j>) = p n+1 + a„p n_1 + 
•■■ + ai; at G C} where TwM ~ all polynomials of degree less than n and Aw on TwM is 
A w = C/W'(p) where ' - d/dp and < f,g > w = Res 00 [f{p)g{p) /W {p)\. Then e - d/dai and 
E = (l/(n + 1)) J2( n ~ * + l)ai(9/dai). This should be compared to Section 11.1 where W — 
x r+1 — uix 1 "^ 1 M r so at ~ — u r -k+i and e ~ —d/du r with n = r. Note the indexing leads to 

W(tj>, t"- fe+2 c fe ) = f i+1 I^(p, ojt) =► (pd p + £(n -k + 2 W^r) W = (n + l)W (10.38) 

and (n — k + 2)ft£, ~ — (n — k + 2)«„_/ £+1 ~ — (to + l)tt m (cf. ( |10. 26j ) ) . One checks the vanishing of 
the curvature for the metric < f,g >w— R e Soo[fg/W] as follows. Consider p — p(W) inverse to 
W = W(p) obtained via Puiseaux series 

1 I 't n t n ~ x t 1 \ ( 1 \ 

p = p ^ = k + —{T + l^ + --- + k;) +0 {k^) ^ 

where k n+1 — W; this determines the coefficients t l (a^) where 

p(k) n+1 + anpik) 11 - 1 + ■ ■ ■ + ai = k n+1 = W (10.40) 



will determine the expansion ( 10.39| ). There is then a triangular change of coordinates (L) 



-V- + fi(t l+1 , ■ ■ ■ ,t n ) for i = 1, • • • , n. Evidently 



f" = U + 1 flesoo VF 2 ^ 1 (p)^ (10.41) 

?i — a + 1 



G7 



(cf. (10.7)) which suggests that a = n — i + 1 and q = — (n + l)/(n — i + 1) in Section 10.1 and 

we identify t k and T k ). One can verify ( |10.41| ) by looking at dp = dk - (l/(n + l))[{t n /k 2 ) H h 

(nt 1 /fc n+1 )]dfc + 0(dk/k n+2 ) with = To prove that the t a are fiat coordinates one uses the 

thermodynamic identity 

d a {Wd P )\ p=c = -d a (pdW)\ w = c (10.42) 

which follows from ^(^(W, t), t) — W via d a W\ p=c + {d p W)d a p\w=c = which says d a (Wdp)\ p=c + 
d a (pdW) | w=c = since d p W ~ dlV/dp. Next we have for 1 < a < n 



d a (Wdp)\ p= 



[k a ~ 1 dk] + 



(10.43) 



where [fdk] + = [f (dk/dp] + dp. To see this note k = W 1 ^ 11 ^ = p + 0(l/p) via ( |10.4C| ) and from 
(|10.40D and (|10.42|) we have 



- d a (Wdp)\ p=c = d a (pdW)\ w=c = (8 a p)dk n+1 
1 1 „ ( 1 



(10.44) 



1 k n - a+1 



O 



(Ik 



n+1 



L dk + 



The left side is polynomial inp and [0(l/fc)dfc]+ = since dk = dp+0(l/p 2 )dp while (1/fc) = 0(l/p) 
so (10.43) is proved. Now use (10.44) in a standard formula (cf. 153]) 



< d a , dp >= Res 



d a (W(p)dp)dp(W(p)d P ) 



Re 



dW{p) 

d a (Wdp)dp(Wdp)d 7 (Wdp) 



(10.45) 



which yields 



< d a , dp >= Res p=c 



dpdW{p) 
k a - x kP- x dk 



1 



'Ja+P,n+l 



,p=oc ~n + l 

so the t a are flat coordinates since r\ a p is constant. The corresponding F arises via 



(10.46) 



d a F 



{a + l)(n + a + 2) 



(10.47) 



(based on (39[|3)). Thus one will have some WDVV equations based on the LG model corresponding 
to those specified by Wq in Section 10.1. Strictly speaking there is no Whitham theory here - only 
a LG model based perhaps on a dispersionless KdV hierarchy (see Remark 10.3). 

10.3 Witten-Dijkgraaf-Verlinde-Verlinde (WDVV) equations 



We go first to pi2| and continue the context of Section 10.1. The main point is to show that the 
WDVV equations for SW theory (involving variables ai) arise from those of the corresponding LG 
model (this is quite cute). Thus one takes z + (fi 2 / 1 z) = Wg{x,T\, ■ ■ ■ ,T r ) as in (10.2) - ( 10. 3j ) 
but Wg is now expressed in terms of flat c oordin ates a s in (10.7) (again only A n is considered). 



We tent atively identify T a and t a in ([L0.7D and (10.4l|) , where 



in ( 10.8 ) with <j) r = 1 and flatness of rjij in ( |l0.10 ) (where P ~ r = 1) implies ( 10.12 ), namel y 



One writes 



a.s 



dxQij(x) = didjW(x) with rjij —< fa^j^r >— S ei + ej ,h- Associativity of the chiral ring (cf. ( 10.25] )) 
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MMk) = (<f>i<f>j)<f>k implies that C£C™ = C e jk C% or (M) [C u Cj] = where (C,)J = C*. From 
fyfc = CfjT]£k one obtains then WDVV in the form 



i)jk 



F, 



(10.48) 



This is all based on TFT for the LG model with W ~ P in Ql0.2fl - fll0.4j ) or in Remark 10.2. 

Now look at the SW theory based on W with (N) /x 2 = A 2JV /4 where TV ~ h ~ r + 1. We 
used ^ 2 = A 27V in Section 10.1 (as in [jm) (9712018)) but switch now to (N) plus A = X sw = 
(l/2wi)(xdz/z) (instead of A = xdz/z in Sectio n 10.1) in order to conform to the notation of |112| 



(9803126). The PF equati ons (|10.34| ) - (J10.35| ) can be written in flat coordinates as (cf. |112 
simply integrate § Xsw in ( 10-24 ) and ( 10.31 )) 



C U = [J2 liTidi -in- 4 M 2 /i 2 



d 2 u 

dT} 



(10.49) 



£. y -n = didjU - c-jdkdrU = o {d t 



_d_ 

dT, 



) 



Now one makes a change of variables Tj — ► a« = §. A (philosophy below) so that (di — d/dTi 

\ a 2 n 



a. J 



darda. 



Pf.™ =0 
13 daj 



(10.50) 



where = didjaj — Y^i Cijdk®r a i ■ Since aj satisfies Cijaj = one knows that II = aj satisfies 
(10.50) and = 0. Next take H — af — dT '/dai to get the third order equation for T ~ F 
(Pli = 0) 



sw 



■Fijk — CjjFirk', ^ijk — d i a I djajdkaKJ~ijK\ 3~ijk 



d 3 T(a) 
daidajdax 



(10.51) 



Defining a metric by Qij = Tij r one has (O) Ti — CiQ for Ti = {T)jk = Tijk and from commuta- 



tivity of the C, there results 



(10.52) 



Hence the Q~ x Ti commute and consequently the matrices (P) T k x Ti = (G~ 1 Tk)~ 1 G~ 1 Ti also 
commute for fixed k. Therefore we obtain (Q) T%T^ Tj = TjT k x Ti and removing the Jacobians 
dai/dTi from (Q) implies the general WDVV equations 



TiT^Tj = TjT k x Ti 



(10.53) 



as m (T| [EH, f[4{| gH, |160[ [l6§ where an association a/ ~ dwj ~ holomorphic differential is used 



in the constructions and proof (cf. also []26] |3§ |5J, |ll] |L27j |128|, |129|, |143|, |144| for WDVV). Here 
jF ~ pSW an( j ^.jjg corr esponding WDVV equations are a direct consequence of the associativity of 
the chiral ring in the A n LG model (hence of WDVV equations of the form ( 10.36| )). 

Regarding philosophy note we have assumed no a priori connection between F and T. F comes 
from the TFT for the LG model while T is defined via af = dT/dai = § B Xsw and (general) 
WDVV for T follows from the PF equations. In particular F is not related to F w ~ Whitham 
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prepotential for the SW curve. There is also a WDVV theory for a Whitham hierarchy on a RS 
involving the LG type for 1 < k < n plus other variables including the ai (see Remark 10.3 
below) and we note that the metric is quite different from Q insofar as the a,{ are concerned. 



REMARK 10.3. In this direction one can make the following comments. In (10.7) or (10.41) 
we have a formula (AS) T rn — — [(n + l)/(n+l — m)]ResW 1 ~^ m ^ n+1 ^(p)dp for LG times based on a 

superpotential (AT) W = p n+1 — qip 11 ^ 1 ■ — q n [qu ~ uu and a/. = — q n -k+i)- This expresses the 

T m as functions of the q^ and one requires m = 1, • • • , n so there are n primary times and n coefficients 
<7fc. Such objects W as above often arise from the dispersionless form of n-KdV situations for 
example and one can refer to them as dispersionless times for a TFT of LG type. Note that diagonal 
coordinates (Riemann invariants) are obtained via u, = W(pi) where W'(pi) = (i = 1, • • • , n). 
Further in |33| for example one shows how this extends, in the context of Hurwitz spaces, to a new 



class of Whitham times as follows (note this differs from Section 7 and see also [127, 128, 129 1). 
Consider g-g&p solutions of KdV problems based on L = d n+1 — q\d n ~ 1 — ■ ■ ■ — q n and let M ffjn+ i 
be the space of such g-gap solutions. This moduli space M. = M gn+ i is the moduli space of 
algebraic curves T, g (with fixed homology) whose ramification is determined by the meromorphic 
function (AT). We recall by Riemann- Roch that [# (zeros) — #(poles)](dW) — 2g — 2 and for 
p ~ 1/z a pole p n dp corresponds to —z~ n ~ 2 dz so this has order n + 2 (not n). Hence dW will have 
N = 2g — 2 + (n+2) = 2g+n zeros pi where dW(pi) — and one defines Ui = W(pi) (i = 1, • • • , N) as 
local coordinates at the branch points pi. The averaged KdV hierarchy (or Whitham hierarchy) then 
has the form (AU) d m d£li = dxdil m (or more generally d m dft s = d s dil m ) where the dVt m corresond 
to dW m K n+1 ^ + regular terms as p — > oo with § A rfO m = (see below - note for z ~ l/p near oo 
we want to deal with dil m ~ ~mz~ rn ~ 1 dz + regular terms. Further one identifies p here with Q±, 
i.e. dp = and defines a flat metric via (AV) gu — Res Pi (dp 2 / dW) as ds 2 = ~^2i gu(u)du 2 and 
the corresponding flat times U (1 < i < N — 2g + n) for ds 2 have the form 

W r v \i-i/(n+i) 

t t = -In + l)i?esoo — — dp (i = 1, • • • ,n); (10.54) 

n + 1 — i 

t n +a = tt— f pdW \ t g+n+a = (I dp (a = 1, • • • , g) 
2tt« J Aa J Ba 



(as indicated earlier t l ~ ti ~ Ti for 1 < i < n). Then ( 10.46| ) holds, < d n + a , d g + n +f3 >= 3a, p 



and otherwise < , > is zero. Thus the set of TFT times (AS) is enlarged to co ntain 2g additi onal 



primary times and this is a basic sort of Whitham theory (cf. also |2E], |53j, 127, 128, 129 for 
further enlargments). The associated Whitham-LG theory for A n now involves a coupling space 
M = M gt . n+ i with flat coordinates Ti, • • • , Tjy as above and primary fields 4>\ ~ —dp, fo, • • • , 4>n 
of the form (AW) (pi = — (n + l)dQi (i = 1 • • • ,?i) with 4> n + a = duj a and </> ff + rl + Q = da a (a = 
1, • • • , g) where the duj a are holomorphic differentials with § A dio^ — 2iriSjk and there are multivalued 

holomorphic differentials da a = da^ satisfying AB a da a = —dW and § A da a = (where AB a da a — 

da a (P + B a ) — do- a (P)). Further the LG potential is W — W(p) where p ~ Jp g dp ~ Jp g dQi and 



d(Wdp) 



dt a 



d{pdW) 



p—c 

with correlation functions (note the minus sign in (fn above) 

\ ^ 



dt a 



(10.55) 

W=c 



< 4> a (j)p >= r) a p = Res dw=0 — ; (10.56) 
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< 



>= X! ReS dW=Q-J } 



dWdp 



The chiral algebra c^(i) has the form (AX) a <^/3 = c^^dp modulo dW divisible differentials and 

one will have an expansion (AY) pdW = (n + l)df2 n +2 + ta4>a for the multivalued differential 
pdW . Writing k n+1 = W one has also 



pdW = kdW - + Oik- 1 ) 



(10.57) 



along with (s = 1, • • • , g) 



pdW 



2mt n+s ; A As (pdW) 



0; A Bs (pdW)=t 



n+g+s 



dW 



(10.58) 



In keeping with dn m ~ dW m ^ n+1 ^> as indicated and (AW) one has <f> m = (-k™- 1 + 0(k- 2 ))dk 
with f. m = 0. One can also determine a partition function or prepotential F such that c Q ^ 7 = 
d^F/dtadpd-y. This development is an obvious ancestor to the SW theories of more recent vintage. 
The analogous SW theory is best developed in a two puncture Toda framework with two points 
oo± corresponding to p 



31 105, 10e, 127, 128, 129, h6§); the analogous W term in 



oo (cf. I 

not polynomial but logarithmic which severly limits the number of primary times. We will not 
develop this further here. 
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